Next Article in Journal
Development of a Non-Hydroxamate Dual Matrix Metalloproteinase (MMP)-7/-13 Inhibitor
Next Article in Special Issue
The Oligomeric Form of the Escherichia coli Dps Protein Depends on the Availability of Iron Ions
Previous Article in Journal
Valorization of By-Products from Commercial Fish Species: Extraction and Chemical Properties of Skin Gelatins
Previous Article in Special Issue
The Role of Metal Binding in the Amyotrophic Lateral Sclerosis-Related Aggregation of Copper-Zinc Superoxide Dismutase
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 22
Issue 9
10.3390/molecules22091549
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Understanding of MYB Transcription Factors Involved in Glucosinolate Biosynthesis in Brassicaceae

by
Mi-Suk Seo
and
Jung Sun Kim
*
Genomics Division, Department of Agricultural Bio-Resources, National Institute of Agricultural Sciences, Rural Development Administration, Wansan-gu, Jeonju 54874, Korea
*
Author to whom correspondence should be addressed.
Molecules 2017, 22(9), 1549; https://doi.org/10.3390/molecules22091549
Submission received: 12 August 2017 / Revised: 8 September 2017 / Accepted: 12 September 2017 / Published: 14 September 2017
(This article belongs to the Special Issue Metallopeptides)
Browse Figures
Versions Notes

Abstract

:
Glucosinolates (GSLs) are widely known secondary metabolites that have anticarcinogenic and antioxidative activities in humans and defense roles in plants of the Brassicaceae family. Some R2R3-type MYB (myeloblastosis) transcription factors (TFs) control GSL biosynthesis in Arabidopsis. However, studies on the MYB TFs involved in GSL biosynthesis in Brassica species are limited because of the complexity of the genome, which includes an increased number of paralog genes as a result of genome duplication. The recent completion of the genome sequencing of the Brassica species permits the identification of MYB TFs involved in GSL biosynthesis by comparative genome analysis with A. thaliana. In this review, we describe various findings on the regulation of GSL biosynthesis in Brassicaceae. Furthermore, we identify 63 orthologous copies corresponding to five MYB TFs from Arabidopsis, except MYB76 in Brassica species. Fifty-five MYB TFs from the Brassica species possess a conserved amino acid sequence in their R2R3 MYB DNA-binding domain, and share close evolutionary relationships. Our analysis will provide useful information on the 55 MYB TFs involved in the regulation of GSL biosynthesis in Brassica species, which have a polyploid genome.

1. Introduction

Plants produce various secondary metabolites that are involved in traits such as taste, color, and scent, and have roles in plant defense against environmental changes or stress unrelated to the primary functions of plants, such as development, reproduction, and photosynthesis [1]. More than 200,000 secondary metabolites are known in plants, and humans have utilized plants for their various benefits (natural medicines, flavors, insecticides, and industrial materials) obtained through the wide chemical diversity of metabolites [2]. Recent studies have reported that cruciferous (Brassicaceae) vegetables are rich in secondary metabolites, including carotenoids, flavonoids, anthocyanins, and glucosinolates (GSLs) [3]. GSLs are derived from amino acids and sugars and are one of the largest known groups of secondary metabolites in the Brassicaceae family. GSLs and their breakdown products have recently attracted attention owing to their various beneficial roles, such as anticarcinogenic and antioxidative activities in humans, and defense against pests and pathogens in plants [2,4]. Therefore, understanding of the GSL biosynthesis pathway will increase their nutritional value and provide agriculturally useful information related to the defense mechanisms in Brassicaceae plants.
The Brassicaceae family consists of 338 genera and about 3700 species, including Arabidopsis thaliana, which is widely studied as a model plant with a small genome, and many plants with agronomic importance such as vegetables, fodder, oil crops [5]. The agricultural and nutritional properties of Brassicaceae plants have resulted in their extensive cultivation. The genus Brassica, which is a member of Brassicaceae, contains the vegetable species, Brassica rapa (A genome, Chinese cabbage, pak choi), Brassica oleracea (C genome, kale, broccoli, cauliflower), Brassica nigra (B genome, black mustard), and Brassica napus (A and C genome, canola). These Brassica species are known as the ’triangle of U’ regarding the genome evolution and the relationships between members of the Brassica genus, including three diploid and three allopolyploid genomes. The genomes of these Brassica species contain triplicated whole genomes of A. thaliana [6]. Thus, a comprehensive approach including A. thaliana and Brassica species will provide useful information on polyploidy-related genome evolution. Recently, the genome sequences of B. rapa and B. oleracea, with diploid genomes, and B. juncea and B. napus, with allopolyploid genomes, have been published and used to explain the evolutionary relationships between Brassica species [7,8,9,10]. The completion of the polyploid genome sequence of the Brassica species permits the identification of various gene families of secondary metabolites by comparative genome analysis with A. thaliana.
Recent studies have shown that genome polyploidy leads to the amplification or redundancy of genes, and duplicated genes with functional diversity distinct from their ancestral function [11]. One of the representative functional genes in late flowering, FLC (FLOWERING LOCUS C) have five copies in the B. rapa genome; these five copies are the result of genome triplication (BrFLC1, 2, 3) and additive recent duplication (BrFLC3a, b) as well as gene introgression (BrFLC5) [12]. In our previous study, BrFLCs showed different expression levels in tissue (leaf and root) and in unvernalized and vernalized B. rapa [13]. Such variation in the copy number of genes could play a major role in regulatory processes in the Brassica genus. A genome-wide comparison of genes related to GSL biosynthesis suggested that genome duplication resulted in the expansion of GSLs-related genes in the Brassica genome [10,14,15]. Furthermore, functional analysis of these genes has demonstrated that GSL biosynthesis in Brassica species is controlled by more regulators than in Arabidopsis, because of their polyploid genomes [16,17]. This indicates that GSL biosynthesis of Brassica species has a regulatory network distinct from Arabidopsis. Recently, some MYB TFs were described as the regulatory factors directly or indirectly activating GSL pathway genes [18,19,20]. Studies have used various molecular biological approaches to investigate GSL biosynthesis in A. thaliana and Brassica species [16,17,18]. Some R2R3 MYB transcription factors (TFs) with two MYB repeats (R2R3) control GSL biosynthesis. In particularly, MYB28, MYB29, and MYB76 are positive regulators of aliphatic GSL synthesis in Arabidopsis [19]. Furthermore, MYB34, MYB51, and MYB122 have been shown to regulate indolic GSL synthesis [21].
In this review, we describe these various findings to provide an integrative understanding of divergence and a functional diversification of the MYB TFs involved in GSL biosynthesis following genome duplication in A. thaliana and Brassica species.

2. Metabolic and Regulatory Pathways of GSL Biosynthesis

2.1. Metabolic Pathway of GSLs

GSLs are a large group of sulfur-containing secondary metabolites, with more than 120 different GSLs identified in Brassicaceae [22]. GSLs are classified into three groups derived from various amino acids: (1) the aliphatic group, derived from Met, Ala, Leu, Ile or Val; (2) the indolic group, derived from Trp; and (3) the aromatic group, derived from Phe or Tyr [23]. The GSLs of the three groups are synthesized in three stages (Figure 1). The precursor amino acids of the aliphatic and aromatic groups are chain-elongated in the early phase, and then oxime is converted into the core GSL structure in the second stage, respectively. Members of the indolic group are transformed into the core GSL structure from the oxime at the first stage, without undergoing the chain-elongation stage. Finally, the side chains of the core structures are modified by oxidation, elimination, alkylation or esterification in the aliphatic and indolic groups [14]. Aliphatic GSLs are the most abundant, at about 57–97% of the total GSL content, while aromatic GSLs are minor components in Brassica species [24,25,26]. The GSLs of the aliphatic and aromatic groups have unique properties for human health and livestock and have roles in oncogenesis, disease, and nutrition [27]. Glucoraphanin, a precursor of sulforaphane, can act as an anticarcinogenic agent in human cells [28]. Sinigrin has been reported in high amounts in black mustard (B. nigra), broccoli (B. oleracea), and Indian mustard (B. juncea), and has been shown to have antioxidant, anticancer, and antifungal effects [29]. Furthermore, phenethyl isothiocyanate (PEITC), which is formed from gluconasturtiin of the aromatic GSL group, has been shown to have a significant chemo-preventative effect against human prostate cancer [30]. Conversely, progoitrin causes goiter disease in mammals, and impedes the use of Brassica crops for cattle feed [31]. GSLs of the indolic groups, such as glucobrassicin, 4-hydroxy glucobrassicin, neoglucobrassicin, and 4-methoxy glucobrassicin, have demonstrated that these breakdown products contribute to defense against biological stresses, including pests and pathogens [32,33]. Tryptophan-derived indole-3-acetaldoxime in the indolic group—also known as a precursor for auxin related to plant growth and development—and camalexin biosynthesis controlling deter bacterial and fungal pathogens [34,35]. Recently, some aliphatic and indolic GSLs have been reported that are involved in primary metabolism in plants, such as nitrogen and sulfur sources, and abiotic stresses, such as salinity, light, and elevated CO2 [36,37,38,39,40].

2.2. Regulators of the GSL Biosynthesis Pathway in the Model Plant Arabidopsis

Most of the structural genes and transcription factors involved in GSL biosynthesis have been identified from diverse molecular and genetic studies in Arabidopsis [18,19,33]. To date, more than 40 structural genes and eight TFs involved in the GSL biosynthesis pathway have been identified in Arabidopsis. Dof1.1 for DNA binding with one finger has been reported as the regulator of networks that positively control indolic GSLs in A. thaliana [41]. IQD (IQ-domain) 1.1 encodes a novel nuclear-localized protein that positively regulates GSL accumulation, and controls the expression of several GSL pathway genes [42]. A further six TFs belong to the group of MYB TFs that contain a R2R3 DNA-binding domain. MYB28, MYB29, and MYB76 are positive regulators of aliphatic GSLs in A. thaliana, and high levels of MYB28 transcription result in the production of a large amounts of aliphatic GSLs. MYB34, MYB51, and MYB122 have been shown to regulate indolic GSL biosynthesis [18,19,43,44]. Additionally, recent studies have shown that the MYB34, MYB51, and MYB122 MYB TFs are involved in the biosynthesis of jasmonate acid (JA), abscisic acid (ABA), ethylene (ET), and salicylic acid (SA), which are involved in plant defense [21]. Additionally, these MYB TFs regulate GSL biosynthesis in cooperation with MYC-bHLH (MYC-like basic helix-loop-helix) TFs, known as signaling components of the jasmonic acid pathway [45]. The R2R3 MYB TFs influence the expression of biosynthetic genes related to GSL biosynthesis in A. thaliana [19,46]. AtMYB28 has been reported to regulate the expression of aliphatic GSL biosynthetic genes, such as AtBCAT-3, AtLeuC1 and MAM1 [19]. In particular, biosynthetic gene AOP2 (2-oxoglutarate-dependent dioxygenases) expression was increased the transcript levels of the MYB28 and MYB29 of the GSL pathway, controlling GSL biosynthesis [47]. The TFs MYB34, MYB51, and MYB122 positively regulate the accumulation of indolic GSLs and have the potential to upregulate GSL biosynthetic genes, such as CYP79Bs [48]. Moreover, all six MYB TFs from A. thaliana was shown to control genes associated with primary sulfate metabolism and are closely related to the GSL biosynthesis network [49].

3. MYB TFs Involved in GSL Biosynthesis in Genome-Sequenced Brassica Species

3.1. Characterization of MYB TFs in Brassica Species

MYB TFs comprise one of the largest gene families of plant TFs, and play significant roles in the regulation of multiple biological processes, such as developmental and environmental responses, and metabolic pathways. These MYB TFs are classified into three subfamilies depending on the number of DNA-binding domain repeats and are known to regulate various biological processes by modulating the rate of transcription initiation of target genes [50]. The R2R3 MYB subfamily, containing the DNA-binding domain of the helix–loop–helix repeats R2 and R3, is the largest group of the MYB family and regulates plant-specific processes including primary and secondary metabolism [50]. Interestingly, the DNA binding domain of the R2R3 MYB family contains conserved amino acid sequence motifs, despite the divergence of the amino acid sequence downstream and the conserved motifs contributing to functional conservation [51]. Therefore, an understanding of R2R3 MYB TFs in Brassica species may help to elucidate the regulation of secondary metabolism by polyploidy-genome evolution in Brassicaceae.
Six MYB TFs belonging to the R2R3 MYB family have been identified as regulators of GSL biosynthesis in Arabidopsis [43,52]. These MYB TFs were also identified using the NCBI (National Center for Biotechnology Information) and Brassica databases (http://brassicadb.org/brad/) of genome-sequenced Brassica species [9,14,20,53]. In total, we identified 63 orthologous copies corresponding to five MYB TFs, except MYB76. To date, MYB76 TFs have not been defined in the genomes of Brassica species, although this has been identified as the positive regulator of aliphatic GSLs in A. thaliana. A summary of these 63 orthotonus gene sequences of MYB TFs related to GSL biosynthesis in Brassicaceae is shown in Table 1. More than two copies corresponding to each MYB TF were found in the Brassica species. These results indicated that genome duplication events have contributed to the expansion of the R2R3 MYB gene family in the Brassica species. The Brassica genus, initiating with B. rapa [7] and ending finally with B. juncea [10] has had its complete genome sequence reported. Only 68% of B. nigra (Accession YZ12151) and 85% of B. oleracea (var capitata line 02-12) of the estimated genome size have been sequenced, and there is some lack of completed gene annotation in their genomes [54]. There are some insufficient gene sequence fragments of orthologous Brassica genes for which we can only confirm 52 full-length genes.
Although 11 partial sequences of MYB TFs exist in Brassica, we were only able to identify the R2R3 MYB DNA-binding domain sequence in BnMYB28.2, BnMYB51.3, and BnMYB51.8. In total, 55 orthologous genes contain a complete R2R3 MYB DNA-binding domain sequence, and these genes showed high conservation at orthologous and paralogous levels (Figure 2). These conserved domains consist of 102 amino acids (AA) and exhibited more than 90% sequence identity with each A. thaliana. However, only three motifs (N[R/K/H]VA) were conserved in the 52 full-length MYB TFs of Brassica species, despite the 12 motifs conserved in the C-terminal regions of six AtMYB TFs. Consequently, variations in gene length (in base pair) and nonsynonymous amino acid sequences are caused by polymorphisms in the downstream C-terminal region. We have identified and reported the conserved R2R3 MYB DNA-binding domains of 13 MYB TFs related to GSL biosynthesis in B. rapa [16,24]. A phylogenetic tree was constructed using the amino acid sequences of 61 R2R3 MYB DNA-binding domains containing A. thaliana to elucidate evolutionary relationships (Figure 3). Fifty-five MYB TFs were grouped based on six AtMYB TFs. Three MYB TFs of Brassica related to aliphatic GSLs, including MYB28, 29 and 76, were subgroups, with 76 located between subgroups 28 and 29. MYB34, MY51, and MYB122, involved in indolic GSLs, clustered in a large group, and the MYB34 group formed a subgroup with MYB51. These results indicated that the six MYB TFs related to GSL biosynthesis were evolutionarily conserved in Brassica species and exhibit functional conservation. The AtMYB76 showed a high sequence identity with all MYB29 TFs in the conserved R2R3 MYB DNA-binding domain, but showed a low sequence identity of 69.3% and 62–66% with AtMYB29 and MYB29 TFs of Brassica species in the full length of the MYB29 TFs (no data). Polyploidy can result in chromosomal rearrangements and gene loss, due to unequal rates of sequence evolution of duplicated genes and changes in DNA methylation [55,56,57]. Although the genome assembly of Brassica species is not complete, the loss of the MYB76 TF in Brassica species may have been caused by genome duplication in evolutionary time. Furthermore, the R2R3 MYB TFs of B. juncea and B. napus, which have allopolyploid genomes, are closely related to those of B. rapa, B. oleracea, and B. nigra, which have diploid genomes. The amino acid sequence identity of the analyzed R2R3 MYB DNA-binding domains with diploid plants showed the evolutionary origin of the allopolyploid plants in Table 2. Four BjMYB28 TFs from B. juncea (AABB) revealed a high level of sequence similarity to B. rapa (AA) and B. nigra (BB). BjMYB28.1 and BjMYB28.3 were found to possess 100% sequence identity with BniMYB28.1 and BrMYB28.3, respectively, while BjMYB28.2, BrMYB28.2, and BniMYB28.2 were fully conserved among the three species. Although BjMYB28.4 was found to possess a high level of sequence similarity (99.02%) to BrMYB28.1, it is possible that an unidentified MYB28 TF of B. nigra exists with higher sequence similarity. Similarly, BnMYB28.1 and BnMYB28.4 from B. napus (AACC) showed 100% sequence similarity with BolMYB28.3 of B. oleracea (CC), and BnMYB28.2 was found to be 100% conserved in three species. BnMYB29s, BnMYB34s, BnMYB51s, and BnMYB122s of B. napus were also found to possess high sequence identity with either BrMYBs or BolMYBs. These results revealed that the BjMYB and BnMYB TFs originated from the genomes of BrMYBs, BniMYBs, or BolMYBs. Therefore, MYB TFs related to GSL biosynthesis are evolutionary closed and conserved in their R2R3 DNA-binding domains despite duplication and hybridization of two diploid Brassica genomes. Furthermore, the number of MYB TFs has increased during evolution, which may have allowed functional diversification and the development of complex networks for the regulation of GSL biosynthesis in polyploidy Brassica species.

3.2. Functional Description of MYB TFs Related to GSL Biosynthesis in Brassica Species

Most studies on MYB TFs related to GSL biosynthesis have been performed using A. thaliana of the Brassicaceae family as a model plant with a small genome size [21,43,44]. The recent completion of the genome sequencing of Brassica species permits the identification of various gene families in its genome. B. rapa is a model dicot plant for use in studies of polyploidy-related genome structure and evolution because of the small size of its genome (529 or 485 kb) in the Brassica genus [54]. Many putative biosynthetic and regulatory genes related to GSL biosynthesis have been identified in the genome sequence of B. rapa [7]. The 13 BrMYB TFs that possess a complete coding sequence indicate that paralogous genes arising through gene duplication have led to functional diversity and changes in expression patterns, reflected by genotype-specific variation in B. rapa subspecies [24]. The expression of some BrMYB TFs, such as BrMYB28, BrMYB34, and BrMYB51 also increased under abiotic and biotic stress conditions. Furthermore, functional analysis of three BrMYB28 TFs has been performed using Agrobacterium-mediated transformation in B. rapa [16]. The three BrMYB28 TFs are involved in the regulation of aliphatic, indolic, and aromatic GSL biosynthesis, and in the expression of biosynthetic genes, such as BrAOP2 and BrGSL-OH in transgenic B. rapa. These results suggested that the regulation of GSL biosynthesis involves a GSL pathway that is more complex than that in Arabidopsis due to the complexity of the polyploidy genome in B. rapa. Four paralogs of BjMYB28 in B. juncea have been reported as regulators of aliphatic GSL accumulation in transgenic A. thaliana and the gene silencing lines of B. juncea with a low GSL content [20,59]. Studies on MYB TFs in B. oleracea are still not sufficient for the mechanism of BolMYB TFs in the regulation of GSLs biosynthesis to be elucidated, although some studies have revealed various expression patterns of biosynthetic genes and BolMYB TFs [51,60]. Genetic studies using association mapping have identified multiple loci controlling GSL biosynthesis in B. napus and B. juncea [61,62]. Recently, the MYB28 TF was identified as the regulator of aliphatic GSL biosynthesis by associative transcriptomics using transcriptome sequencing in B. napus [63]. Additionally, association mapping of B. napus and B. juncea confirmed that MYB28 TFs are associated with GSL content [64,65].
Previous studies have demonstrated the role of MYB 28, 29, and 76 as aliphatic GSLs, and MYB 51, 34 and 122 as indolic GSLs, in the regulation of GSL biosynthesis in A. thaliana. In the case of Brassica species, only a few recent studies have shown that MYB28 TFs positively regulate the accumulation of GSLs. The complexity of the Brassica genome, due to an increased number of paralog genes through genome duplication, makes it difficult for a molecular approach to be used to determine the regulation of GSL biosynthesis. In this review, we discussed the role of MYB TFs as important regulators associated with the GSL biosynthesis pathway in Brassica species, and provided useful information on the 55 MYB TFs for improved understanding of the regulatory mechanism of GSL biosynthesis in Brassica species.

4. Conclusions

Recently, six R2R3 MYB TFs controlling the accumulation of various GSLs were reported as regulators of different stress responses and hormones, such as ABA, ethylene, and jasmonate, in A. thaliana. Although Brassica crops have commercial and scientific value, our understanding of the roles of most MYB TFs is lacking, with the exception of a few MYB TFs related to GSL biosynthesis in Brassica species. The 55 R2R3 MYB TFs identified as ortholog genes with A. thaliana in Brassica species shared a close evolutionary relationship, with a highly conserved DNA-binding amino acid sequence. This will provide valuable information on the mechanisms of MYB TF regulation on unique properties, such as stress responses and various metabolites containing GSL biosynthesis in Brassica species with polyploid genomes. Further extensive functional studies of the 55 MYB TFs will help to elucidate the functional diversity of genes via genome duplication in polyploidy plants.

Acknowledgments

This work was supported by a grant from the National Academy of Agricultural Science (PJ010157), Rural Development of Administration, Korea.

Author Contributions

M.-S. Seo and J.S. Kim wrote this manuscript.

Conflicts of Interest

The authors declare no conflicts of interest

References

  1. Verpoorte, R.; Memelink, J. Engineering secondary metabolite production in plants. Curr. Opin. Biotechnol. 2002, 13, 181–187. [Google Scholar] [CrossRef]
  2. Yonekura-Sakakibara, K.; Saito, K. Functional genomics for plant natural product biosynthesis. Nat. Prod. Rep. 2009, 26, 1466–1487. [Google Scholar] [CrossRef] [PubMed]
  3. Cartea, M.E.; Francisco, M.; Soengas, P.; Velasco, P. Phenolic compounds in brassica vegetables. Molecules 2010, 16, 251–280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Hayes, J.D.; Kelleher, M.O.; Eggleston, I.M. The cancer chemopreventive actions of phytochemicals derived from glucosinolates. Eur. J. Nutr. 2008, 47, 73–88. [Google Scholar] [CrossRef] [PubMed]
  5. Beilstein, M.A.; Al-Shehbaz, I.A.; Kellogg, E.A. Brassicaceae phylogeny and trichome evolution. Am. J. Bot. 2006, 93, 607–619. [Google Scholar] [CrossRef] [PubMed]
  6. Lysak, M.A.; Koch, M.A.; Pecinka, A.; Schubert, I. Chromosome triplication found across the tribe brassiceae. Genome Res. 2005, 15, 516–525. [Google Scholar] [CrossRef] [PubMed]
  7. Wang, X.; Wang, H.; Wang, J.; Sun, R.; Wu, J.; Liu, S.; Bai, Y.; Mun, J.H.; Bancroft, I.; Cheng, F.; et al. The genome of the mesopolyploid crop species Brassica rapa. Nat. Genet. 2011, 43, 1035–1039. [Google Scholar] [CrossRef] [PubMed]
  8. Liu, S.; Liu, Y.; Yang, X.; Tong, C.; Edwards, D.; Parkin, I.A.; Zhao, M.; Ma, J.; Yu, J.; Huang, S.; et al. The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat. Commun. 2014, 5, 3930. [Google Scholar] [CrossRef] [PubMed]
  9. Chalhoub, B.; Denoeud, F.; Liu, S.; Parkin, I.A.; Tang, H.; Wang, X.; Chiquet, J.; Belcram, H.; Tong, C.; Samans, B.; et al. Early allopolyploid evolution in the post-neolithic Brassica napus oilseed genome. Science 2014, 345, 950–953. [Google Scholar] [CrossRef] [PubMed]
  10. Yang, J.; Liu, D.; Wang, X.; Ji, C.; Cheng, F.; Liu, B.; Hu, Z.; Chen, S.; Pental, D.; Ju, Y.; et al. The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection. Nat. Genet. 2016, 48, 1225–1232. [Google Scholar] [CrossRef] [PubMed]
  11. Adams, K.L.; Wendel, J.F. Polyploidy and genome evolution in plants. Curr. Opin. Plant Biol. 2005, 8, 135–141. [Google Scholar] [CrossRef] [PubMed]
  12. Kim, J.S.; Chung, T.Y.; King, G.J.; Jin, M.; Yang, T.J.; Jin, Y.M.; Kim, H.I.; Park, B.S. A sequence-tagged linkage map of Brassica rapa. Genetics 2006, 174, 29–39. [Google Scholar] [CrossRef] [PubMed]
  13. Kim, S.Y.; Park, B.S.; Kwon, S.J.; Kim, J.S.; Lim, M.H.; Park, Y.D.; Kim, D.Y.; Suh, S.C.; Jin, Y.M.; Ahn, J.H.; et al. Delayed flowering time in Arabidopsis and Brassica rapa by the overexpression of FLOWERING LOCUS C (FLC) homologs isolated from Chinese cabbage (Brassica rapa L. ssp. pekinensis). Plant Cell Rep. 2007, 26, 327–336. [Google Scholar] [CrossRef] [PubMed]
  14. Zang, Y.X.; Kim, H.U.; Kim, J.A.; Lim, M.H.; Jin, M.; Lee, S.C.; Kwon, S.J.; Lee, S.I.; Hong, J.K.; Park, T.H.; et al. Genome-wide identification of glucosinolate synthesis genes in Brassica rapa. FEBS J. 2009, 276, 3559–3574. [Google Scholar] [CrossRef] [PubMed]
  15. Guo, R.; Huang, Z.; Deng, Y.; Chen, X.; XuHan, X.; Lai, Z. Comparative transcriptome analyses reveal a special glucosinolate metabolism mechanism in Brassica alboglabra sprouts. Front. Plant Sci. 2016, 7, 1497. [Google Scholar] [CrossRef] [PubMed]
  16. Seo, M.S.; Jin, M.; Chun, J.H.; Kim, S.J.; Park, B.S.; Shon, S.H.; Kim, J.S. Functional analysis of three BrMYB28 transcription factors controlling the biosynthesis of glucosinolates in Brassica rapa. Plant Mol. Biol. 2016, 90, 503–516. [Google Scholar] [CrossRef] [PubMed]
  17. Borpatragohain, P.; Rose, T.J.; King, G.J. Fire and brimstone: Molecular interactions between sulfur and glucosinolate biosynthesis in model and crop brassicaceae. Front. Plant Sci. 2016, 7, 1735. [Google Scholar] [CrossRef] [PubMed]
  18. Sonderby, I.E.; Burow, M.; Rowe, H.C.; Kliebenstein, D.J.; Halkier, B.A. A complex interplay of three R2R3 MYB transcription factors determines the profile of aliphatic glucosinolates in Arabidopsis. Plant Physiol. 2010, 153, 348–363. [Google Scholar] [CrossRef] [PubMed]
  19. Hirai, M.Y.; Sugiyama, K.; Sawada, Y.; Tohge, T.; Obayashi, T.; Suzuki, A.; Araki, R.; Sakurai, N.; Suzuki, H.; Aoki, K.; et al. Omics-based identification of Arabidopsis MYB transcription factors regulating aliphatic glucosinolate biosynthesis. Proc. Natl. Acad. Sci. USA 2007, 104, 6478–6483. [Google Scholar] [CrossRef] [PubMed]
  20. Augustine, R.; Majee, M.; Gershenzon, J.; Bisht, N.C. Four genes encoding MYB28, a major transcriptional regulator of the aliphatic glucosinolate pathway, are differentially expressed in the allopolyploid Brassica juncea. J. Exp. Bot. 2013, 64, 4907–4921. [Google Scholar] [CrossRef] [PubMed]
  21. Frerigmann, H.; Gigolashvili, T. MYB34, MYB51, and MYB122 distinctly regulate indolic glucosinolate biosynthesis in Arabidopsis thaliana. Mol. Plant 2014, 7, 814–828. [Google Scholar] [CrossRef] [PubMed]
  22. Hopkins, R.J.; van Dam, N.M.; van Loon, J.J. Role of glucosinolates in insect-plant relationships and multitrophic interactions. Annu. Rev. Entomol. 2009, 54, 57–83. [Google Scholar] [CrossRef] [PubMed]
  23. Kliebenstein, D.J.; Kroymann, J.; Brown, P.; Figuth, A.; Pedersen, D.; Gershenzon, J.; Mitchell-Olds, T. Genetic control of natural variation in Arabidopsis glucosinolate accumulation. Plant Physiol. 2001, 126, 811–825. [Google Scholar] [CrossRef] [PubMed]
  24. Seo, M.S.; Jin, M.; Sohn, S.H.; Kim, J.S. Expression profiles of BrMYB transcription factors related to glucosinolate biosynthesis and stress response in eight subspecies of Brassica rapa. FEBS Open Biol. 2017. [Google Scholar] [CrossRef]
  25. Padilla, G.; Cartea, M.E.; Velasco, P.; de Haro, A.; Ordas, A. Variation of glucosinolates in vegetable crops of Brassica rapa. Phytochemistry 2007, 68, 536–545. [Google Scholar] [CrossRef] [PubMed]
  26. Cartea, M.E.; Velasco, P.; Obregon, S.; Padilla, G.; de Haro, A. Seasonal variation in glucosinolate content in Brassica oleracea crops grown in northwestern spain. Phytochemistry 2008, 69, 403–410. [Google Scholar] [CrossRef] [PubMed]
  27. Verkerk, R.; Schreiner, M.; Krumbein, A.; Ciska, E.; Holst, B.; Rowland, I.; De Schrijver, R.; Hansen, M.; Gerhauser, C.; Mithen, R.; et al. Glucosinolates in Brassica vegetables: The influence of the food supply chain on intake, bioavailability and human health. Mol. Nutr. Food Res. 2009, 53. [Google Scholar] [CrossRef] [PubMed]
  28. Brooks, J.D.; Paton, V.G.; Vidanes, G. Potent induction of phase 2 enzymes in human prostate cells by sulforaphane. Cancer Epidemiol. Biomark. Prev. 2001, 10, 949–954. [Google Scholar]
  29. Mazumder, A.; Dwivedi, A.; du Plessis, J. Sinigrin and its therapeutic benefits. Molecules 2016, 21, 416. [Google Scholar] [CrossRef] [PubMed]
  30. Powolny, A.; Takahashi, K.; Hopkins, R.G.; Loo, G. Induction of gadd gene expression by phenethylisothiocyanate in human colon adenocarcinoma cells. J. Cell. Biochem. 2003, 90, 1128–1139. [Google Scholar] [CrossRef] [PubMed]
  31. Pearson, A.W.; Greenwood, N.M.; Butler, E.J.; Fenwick, G.R. Biochemical changes in layer and broiler chickens when fed on a high-glucosinolate rapeseed meal. Br. Poult. Sci. 1983, 24, 417–427. [Google Scholar] [CrossRef] [PubMed]
  32. Kim, J.H.; Lee, B.W.; Schroeder, F.C.; Jander, G. Identification of indole glucosinolate breakdown products with antifeedant effects on Myzus persicae (green peach aphid). Plant J. 2008, 54, 1015–1026. [Google Scholar] [CrossRef] [PubMed]
  33. Pfalz, M.; Vogel, H.; Kroymann, J. The gene controlling the indole glucosinolate modifier1 quantitative trait locus alters indole glucosinolate structures and aphid resistance in Arabidopsis. Plant Cell 2009, 21, 985–999. [Google Scholar] [CrossRef] [PubMed]
  34. Zhao, Y.; Hull, A.K.; Gupta, N.R.; Goss, K.A.; Alonso, J.; Ecker, J.R.; Normanly, J.; Chory, J.; Celenza, J.L. Trp-dependent auxin biosynthesis in Arabidopsis: Involvement of cytochrome P450s CYP79b3 and CYP79b3. Gene. Dev. 2002, 16, 3100–3112. [Google Scholar] [CrossRef] [PubMed]
  35. Glawischnig, E.; Hansen, B.G.; Olsen, C.E.; Halkier, B.A. Camalexin is synthesized from indole-3-acetaldoxime, a key branching point between primary and secondary metabolism in Arabidopsis. Proc. Natl. Acad. Sci. USA 2004, 101, 8245–8250. [Google Scholar] [CrossRef] [PubMed]
  36. Marino, D.; Ariz, I.; Lasa, B.; Santamaria, E.; Fernandez-Irigoyen, J.; Gonzalez-Murua, C.; Aparicio Tejo, P.M. Quantitative proteomics reveals the importance of nitrogen source to control glucosinolate metabolism in Arabidopsis thaliana and Brassica oleracea. J. Exp. Bot. 2016, 67, 3313–3323. [Google Scholar] [CrossRef] [PubMed]
  37. Huseby, S.; Koprivova, A.; Lee, B.R.; Saha, S.; Mithen, R.; Wold, A.B.; Bengtsson, G.B.; Kopriva, S. Diurnal and light regulation of sulphur assimilation and glucosinolate biosynthesis in Arabidopsis. J. Exper. Bot. 2013, 64, 1039–1048. [Google Scholar] [CrossRef] [PubMed]
  38. Schonhof, I.; Klaring, H.P.; Krumbein, A.; Schreiner, M. Interaction between atmospheric CO2 and glucosinolates in broccoli. J. Chem. Ecol. 2007, 33, 105–114. [Google Scholar] [CrossRef] [PubMed]
  39. Gigolashvili, T.; Kopriva, S. Transporters in plant sulfur metabolism. Front. Plant Sci. 2014, 5, 442. [Google Scholar] [CrossRef] [PubMed]
  40. Koroleva, O.A.; Cramer, R. Single-cell proteomic analysis of glucosinolate-rich S-cells in Arabidopsis thaliana. Methods 2011, 54, 413–423. [Google Scholar] [CrossRef] [PubMed]
  41. Skirycz, A.; Reichelt, M.; Burow, M.; Birkemeyer, C.; Rolcik, J.; Kopka, J.; Zanor, M.I.; Gershenzon, J.; Strnad, M.; Szopa, J.; et al. Dof transcription factor Atdof1.1 (OBP2) is part of a regulatory network controlling glucosinolate biosynthesis in Arabidopsis. Plant J. 2006, 47, 10–24. [Google Scholar] [CrossRef] [PubMed]
  42. Levy, M.; Wang, Q.; Kaspi, R.; Parrella, M.P.; Abel, S. Arabidopsis IQD1, a novel calmodulin-binding nuclear protein, stimulates glucosinolate accumulation and plant defense. Plant J. 2005, 43, 79–96. [Google Scholar] [CrossRef] [PubMed]
  43. Gigolashvili, T.; Yatusevich, R.; Berger, B.; Muller, C.; Flugge, U.I. The R2R3-MYB transcription factor HAG1/MYB28 is a regulator of methionine-derived glucosinolate biosynthesis in Arabidopsis thaliana. Plant J. 2007, 51, 247–261. [Google Scholar] [CrossRef] [PubMed]
  44. Gigolashvili, T.; Engqvist, M.; Yatusevich, R.; Muller, C.; Flugge, U.I. HAG2/MYB76 and HAG3/MYB29 exert a specific and coordinated control on the regulation of aliphatic glucosinolate biosynthesis in Arabidopsis thaliana. New Phytol. 2008, 177, 627–642. [Google Scholar] [CrossRef] [PubMed]
  45. Frerigmann, H. Glucosinolate regulation in a complex relationship—MYC and MYB—No one can act without each other. Adv. Bot. Res. 2016, 80, 57–97. [Google Scholar]
  46. Sonderby, I.E.; Hansen, B.G.; Bjarnholt, N.; Ticconi, C.; Halkier, B.A.; Kliebenstein, D.J. A systems biology approach identifies a R2R3 MYB gene subfamily with distinct and overlapping functions in regulation of aliphatic glucosinolates. PLoS ONE 2007, 2, e1322. [Google Scholar] [CrossRef] [PubMed]
  47. Burow, M.; Atwell, S.; Francisco, M.; Kerwin, R.E.; Halkier, B.A.; Kliebenstein, D.J. The glucosinolate biosynthetic gene AOP2 mediates feed-back regulation of jasmonic acid signaling in Arabidopsis. Mol. Plant 2015, 8, 1201–1212. [Google Scholar] [CrossRef] [PubMed]
  48. Gigolashvili, T.; Berger, B.; Mock, H.P.; Muller, C.; Weisshaar, B.; Flugge, U.I. The transcription factor HIG1/MYB51 regulates indolic glucosinolate biosynthesis in Arabidopsis thaliana. Plant J. 2007, 50, 886–901. [Google Scholar] [CrossRef] [PubMed]
  49. Yatusevich, R.; Mugford, S.G.; Matthewman, C.; Gigolashvili, T.; Frerigmann, H.; Delaney, S.; Koprivova, A.; Flugge, U.I.; Kopriva, S. Genes of primary sulfate assimilation are part of the glucosinolate biosynthetic network in Arabidopsis thaliana. Plant J. 2010, 62, 1–11. [Google Scholar] [CrossRef] [PubMed]
  50. Stracke, R.; Werber, M.; Weisshaar, B. The R2R3-MYB gene family in Arabidopsis thaliana. Curr. Opin. Plant Biol. 2001, 4, 447–456. [Google Scholar] [CrossRef]
  51. Prouse, M.B.; Campbell, M.M. The interaction between MYB proteins and their target DNA binding sites. Biochim. Biophys. Acta 2012, 1819, 67–77. [Google Scholar] [CrossRef] [PubMed]
  52. Malitsky, S.; Blum, E.; Less, H.; Venger, I.; Elbaz, M.; Morin, S.; Eshed, Y.; Aharoni, A. The transcript and metabolite networks affected by the two clades of Arabidopsis glucosinolate biosynthesis regulators. Plant Physiol. 2008, 148, 2021–2049. [Google Scholar] [CrossRef] [PubMed]
  53. Yi, G.E.; Robin, A.H.; Yang, K.; Park, J.I.; Kang, J.G.; Yang, T.J.; Nou, I.S. Identification and expression analysis of glucosinolate biosynthetic genes and estimation of glucosinolate contents in edible organs of Brassica oleracea subspecies. Molecules 2015, 20, 13089–13111. [Google Scholar] [CrossRef] [PubMed]
  54. Cheng, F.; Liang, J.; Cai, C.; Cai, X.; Wu, J.; Wang, X. Genome sequencing supports a multi-vertex model for Brassiceae species. Curr. Opin. Plant Biol. 2017, 36, 79–87. [Google Scholar] [CrossRef] [PubMed]
  55. Udall, J.A.; Swanson, J.M.; Nettletion, D.; Percifield, R.J.; Wendel, J.F. A novel approach for characterizing expression levels of genes dupicated by polyploidy. Genetics 2006, 173, 1823–1827. [Google Scholar] [CrossRef] [PubMed]
  56. Lukens, L.N.; Pires, J.C.; Leon, E.; Vogelzang, R.; Oslach, L.; Osborn, T. Patterns of sequence loss and cytosine methylation within a population of newly resynthesized Brassica napus allopolyploids. Plant Physiol. 2006, 140, 336–348. [Google Scholar] [CrossRef] [PubMed]
  57. Small, R.I.; Ryburn, J.A.; Wendel, J.F. Low levels of nucleotide diversity at homeologous Adh loci in allotetraploid cotton (Gossypium L.). Mol. Biol. Evolut. 1999, 16, 491–501. [Google Scholar] [CrossRef]
  58. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evolut. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [PubMed]
  59. Augustine, R.; Mukhopadhyay, A.; Bisht, N.C. Targeted silencing of BjMYB28 transcription factor gene directs development of low glucosinolate lines in oilseed Brassica juncea. Plant Biotechnol. J. 2013, 11, 855–866. [Google Scholar] [CrossRef] [PubMed]
  60. Robin, A.H.; Yi, G.E.; Laila, R.; Yang, K.; Park, J.I.; Kim, H.R.; Nou, I.S. Expression profiling of glucosinolate biosynthetic genes in Brassica oleracea l. Var. Capitata inbred lines reveals their association with glucosinolate content. Molecules 2016, 21, 787. [Google Scholar] [CrossRef] [PubMed]
  61. Bisht, N.C.; Gupta, V.; Ramchiary, N.; Sodhi, Y.S.; Mukhopadhyay, A.; Arumugam, N.; Pental, D.; Pradhan, A.K. Fine mapping of loci involved with glucosinolate biosynthesis in oilseed mustard (Brassica juncea) using genomic information from allied species. Theor. Appl. Genet. 2009, 118, 413–421. [Google Scholar] [CrossRef] [PubMed]
  62. Feng, J.; Long, Y.; Shi, L.; Shi, J.; Barker, G.; Meng, J. Characterization of metabolite quantitative trait loci and metabolic networks that control glucosinolate concentration in the seeds and leaves of Brassica napus. New Phytol. 2012, 193, 96–108. [Google Scholar] [CrossRef] [PubMed]
  63. Harper, A.L.; Trick, M.; Higgins, J.; Fraser, F.; Clissold, L.; Wells, R.; Hattori, C.; Werner, P.; Bancroft, I. Associative transcriptomics of traits in the polyploid crop species Brassica napus. Nat. Biotechnol. 2012, 30, 798–802. [Google Scholar] [CrossRef] [PubMed]
  64. Qu, C.M.; Li, S.M.; Duan, X.J.; Fan, J.H.; Jia, L.D.; Zhao, H.Y.; Lu, K.; Li, J.N.; Xu, X.F.; Wang, R. Identification of candidate genes for seed glucosinolate content using association mapping in Brassica napus L. Genes 2015, 6, 1215–1229. [Google Scholar] [CrossRef] [PubMed]
  65. Rout, K.; Sharma, M.; Gupta, V.; Mukhopadhyay, A.; Sodhi, Y.S.; Pental, D.; Pradhan, A.K. Deciphering allelic variations: For seed glucosinolate traits in oilseed mustard (Brassica juncea) using two bi-parental mapping populations. Theor. Appl. Genet. 2015, 128, 657–666. [Google Scholar] [CrossRef] [PubMed]
Sample Availability: Samples of the compounds are available from the authors.
Molecules 22 01549 g001 550
Figure 1. Three major biosynthesis pathways of the glucosinolates in Brassicaceae. Outline of the pathways which can be divided into three parts; (1) amino acid chain elongation; (2) formation of the core structure of GSLs; (3) side chain modification. Blue bold indicates biosynthetic genes involved in each step. The names in parentheses denote common names.
Figure 1. Three major biosynthesis pathways of the glucosinolates in Brassicaceae. Outline of the pathways which can be divided into three parts; (1) amino acid chain elongation; (2) formation of the core structure of GSLs; (3) side chain modification. Blue bold indicates biosynthetic genes involved in each step. The names in parentheses denote common names.
Molecules 22 01549 g001
Molecules 22 01549 g002 550
Figure 2. Amino acid sequence alignment of the R2R3 MYB DNA domains in 55 glucosinolate biosynthetic transcription factors in Brassica species. The R2 and R3 binding domains are boxed in white and black, respectively. At, A. thalinala; Br, B. rapa; Bni, B. nigra; Bol, B. oleracea; Bn, B. napus; Bj, B. juncea.
Figure 2. Amino acid sequence alignment of the R2R3 MYB DNA domains in 55 glucosinolate biosynthetic transcription factors in Brassica species. The R2 and R3 binding domains are boxed in white and black, respectively. At, A. thalinala; Br, B. rapa; Bni, B. nigra; Bol, B. oleracea; Bn, B. napus; Bj, B. juncea.
Molecules 22 01549 g002
Molecules 22 01549 g003 550
Figure 3. Phylogenetic analysis of R2R3 MYB DNA-binding domains in 55 MYB transcription factors involved in the glucosinolate biosynthesis pathway in Brassica species. This tree was constructed using MEGA, version 6, software [58]. Bootstrap values with 1000 replicates are denoted as percentages. At, A. thalinala; Br, B. rapa; Bni, B. nigra; Bol, B. oleracea; Bn, B. napus; Bj, B. juncea.
Figure 3. Phylogenetic analysis of R2R3 MYB DNA-binding domains in 55 MYB transcription factors involved in the glucosinolate biosynthesis pathway in Brassica species. This tree was constructed using MEGA, version 6, software [58]. Bootstrap values with 1000 replicates are denoted as percentages. At, A. thalinala; Br, B. rapa; Bni, B. nigra; Bol, B. oleracea; Bn, B. napus; Bj, B. juncea.
Molecules 22 01549 g003
Table
Table 1. DNA sequence summary of the MYB transcription factors related to glucosinolate biosynthesis in Brassicaceae.
Table 1. DNA sequence summary of the MYB transcription factors related to glucosinolate biosynthesis in Brassicaceae.
SpeciesGene NameGene IDCDS Length (bp)
A. thaliana (diploid)AtMYB28AT5G614201101
AtMYB29AT5G076901011
AtMYB34AT5G60890888
AtMYB51AT1G185701059
AtMYB76AT5G077001017
AtMYB122AT1G740801002
B. rapa (diploid)BrMYB28.1Bra0129611065
BrMYB28.2Bra0359291074
BrMYB28.3Bra0293111119
BrMYB29.1Bra005949993
BrMYB29.2Bra009245267 *
BrMYB34.1Bra035954909
BrMYB34.2Bra013000951
BrMYB34.3Bra029349843
BrMYB34.4Bra029350930
BrMYB51.1Bra031035963
BrMYB51.2Bra0165531002
BrMYB51.3Bra0256661026
BrMYB122.1Bra015939981
BrMYB122.2Bra0081311005
B. oleracea (diploid)BolMYB28.1XM_0137661401059
BolMYB28.2XM_0137508381083
BolMYB28.3XM_0137380781068
BolMYB28.4Bol036743426 *
BolMYB29.1Bol043899324 *
BolMYB29.2XM_013771134993
BolMYB34.1Bol007760843
BolMYB34.2Bol017062951
BolMYB34.3XM_013754685882
BolMYB51.1Bol030761990
BolMYB51.2Bol0132071002
BolMYB122Bol026204981
B. nigra (diploid)BniMYB28.1JX9478411053
BniMYB28.2JX9478421095
B. juncea (allopolyploid)BjMYB28.1JQ6661661053
BjMYB28.2JQ6661671095
BjMYB28.3JQ6661681065
BjMYB28.4JQ6661691065
BjMYB29.1JX316031993
BjMYB29.2JX3160321020
B. napus (allopolyploid)BnMYB28.1GSBRNA2T000409130011011
BnMYB28.2GSBRNA2T00146148001420 *,1
BnMYB28.3GSBRNA2T00146147001555 *
BnMYB28.4GSBRNA2T00113547001987
BnMYB29.1GSBRNA2T00092496001867 *
BnMYB29.2GSBRNA2T00134508001966
BnMYB29.3GSBRNA2T001361490011005 *
BnMYB29.4GSBRNA2T00129660001948
BnMYB34.1GSBRNA2T00052239001843
BnMYB34.2GSBRNA2T00003879001948
BnMYB34.3GSBRNA2T00075530001843
BnMYB34.4GSBRNA2T00075529001918
BnMYB34.5GSBRNA2T00146117001882
BnMYB34.6GSBRNA2T00071741001867
BnMYB34.7GSBRNA2T00094719001909
BnMYB34.8GSBRNA2T00113599001951
BnMYB51.1GSBRNA2T00136837001963
BnMYB51.2GSBRNA2T001395960011002
BnMYB51.3GSBRNA2T00139597001480 *,1
BnMYB51.4GSBRNA2T00006526001990
BnMYB51.5GSBRNA2T00016192001963
BnMYB51.6GSBRNA2T001029300011002
BnMYB51.7GSBRNA2T00035671001309 *
BnMYB51.8GSBRNA2T00035672001477 *,1
BnMYB122.1GSBRNA2T00088839001870 *
BnMYB122.2GSBRNA2T001034850011113
BnMYB122.3GSBRNA2T00082812001960
BnMYB122.4GSBRNA2T00102693001981
BnMYB122.5GSBRNA2T00092372001957
* Partial sequence; 1 Gene contain R2R3 MYB DNA-binding domain of partial sequence.
Table
Table 2. Amino acid sequence identity (%) for the allopolyploidy R2R3 MYB DNA-binding domains of diploidy genes of Brassica. Values represent the percentage of sequence similarity and highest percentage demarcated by the gray blocks.
Table 2. Amino acid sequence identity (%) for the allopolyploidy R2R3 MYB DNA-binding domains of diploidy genes of Brassica. Values represent the percentage of sequence similarity and highest percentage demarcated by the gray blocks.
AllopolyploidDiploid
BrMYB28.1BrMYB28.2BrMYB28.3BniMYB28.1BniMYB28.2
BjMYB28.192.1697.0697.0610097.06
BjMYB28.295.110098.0497.06100
BjMYB28.393.1498.0410097.0698.04
BjMYB28.499.0296.0894.1293.1496.08
BrMYB28.1BrMYB28.2BrMYB28.3BolMYB28.2BolMYB28.3
BnMYB28.199.0296.0894.1296.08100
BnMYB28.295.110098.0410096.08
BnMYB28.499.0296.0894.1296.08100
BrMYB29.1BolMYB29.2
BnMYB29.297.0699.02
BnMYB29.410096.08
BrMYB34.1BrMYB34.2BrMYB34.3BrMYB34.4BolMYB34.1BolMYB34.2BolMYB34.3
BnMYB34.190.294.1296.0891.1810092.1691.18
BnMYB34.293.1498.0493.1490.292.1698.0494.12
BnMYB34.390.294.1298.0490.296.0892.1691.18
BnMYB34.489.2292.1691.1898.0491.1890.290.2
BnMYB34.598.0496.0892.1690.291.1894.12100
BnMYB34.690.293.1492.1699.0292.1691.1891.18
BnMYB34.710095.191.1889.2290.293.1498.04
BnMYB34.895.110095.192.1694.1298.0496.08
BrMYB51.1BrMYB51.2BrMYB51.3BolMYB51.1BolMYB51.2
BnMYB51.199.0297.0694.1294.1297.06
BnMYB51.297.0698.0493.1493.14100
BnMYB51.397.0699.0293.1493.1499.02
BnMYB51.495.193.1410010093.14
BnMYB51.510096.0895.195.197.06
BnMYB51.696.0810093.1493.1498.04
BnMYB51.895.193.1410010093.14
BrMYB122.1BrMYB122.2BolMYB122
BnMYB122.210094.12100
BnMYB122.395.197.0695.1
BnMYB122.410094.12100
BnMYB122.594.1210094.12

Share and Cite

MDPI and ACS Style

Seo, M.-S.; Kim, J.S. Understanding of MYB Transcription Factors Involved in Glucosinolate Biosynthesis in Brassicaceae. Molecules 2017, 22, 1549. https://doi.org/10.3390/molecules22091549

AMA Style

Seo M-S, Kim JS. Understanding of MYB Transcription Factors Involved in Glucosinolate Biosynthesis in Brassicaceae. Molecules. 2017; 22(9):1549. https://doi.org/10.3390/molecules22091549

Chicago/Turabian Style

Seo, Mi-Suk, and Jung Sun Kim. 2017. "Understanding of MYB Transcription Factors Involved in Glucosinolate Biosynthesis in Brassicaceae" Molecules 22, no. 9: 1549. https://doi.org/10.3390/molecules22091549

APA Style

Seo, M. -S., & Kim, J. S. (2017). Understanding of MYB Transcription Factors Involved in Glucosinolate Biosynthesis in Brassicaceae. Molecules, 22(9), 1549. https://doi.org/10.3390/molecules22091549

Article Metrics

Back to TopTop