Novel insights into the function of Arabidopsis R2R3-MYB transcription factors regulating aliphatic glucosinolate biosynthesis

LED Lights Influenced Phytochemical Contents and Biological Activities in Kale (Brassica oleracea L. var. acephala) Microgreens

Light-emitting diodes (LEDs) are regarded as an effective artificial light source for producing sprouts, microgreens, and baby leaves. Thus, this study aimed to investigate the influence of different LED lights (white, red, and blue) on the biosynthesis of secondary metabolites (glucosinolates, carotenoids, and phenolics) and the biological effects on kale microgreens. Microgreens irradiated with white LEDs showed higher levels of carotenoids, including lutein, 13-cis-β-carotene, α-carotene, β-carotene, and 9-cis-β-carotene, than those irradiated with red or blue LEDs. These findings were consistent with higher expression levels of carotenoid biosynthetic genes (BoPDS and BoZDS) in white-irradiated kale microgreens. Similarly, microgreens irradiated with white and blue LEDs showed slightly higher levels of glucosinolates, including glucoiberin, progoitrin, sinigrin, and glucobrassicanapin, than those irradiated with red LEDs. These results agree with the high expression levels of BoMYB28-2, BoMYB28-3, and BoMYB29 in white- and blue-irradiated kale microgreens. In contrast, kale microgreens irradiated with blue LEDs contained higher levels of phenolic compounds (gallic acid, catechin, ferulic acid, sinapic acid, and quercetin). According to the total phenolic content (TPC) and 2,2-diphenyl-1-picrylhydrazyl (DPPH) inhibition assays, the extracts of kale microgreens irradiated with blue LEDs had slightly higher antioxidant activities, and the DPPH inhibition percentage had a positive correlation with TPC in the microgreens. Furthermore, the extracts of kale microgreens irradiated with blue LEDs exhibited stronger antibacterial properties against normal pathogens and multidrug-resistant pathogens than those irradiated with white and red LEDs. These results indicate that white-LED lights are suitable for carotenoid production, whereas blue-LED lights are efficient in increasing the accumulation of phenolics and their biological activities in kale microgreens.

Transcriptome and QTL mapping analyses of major QTL genes controlling glucosinolate contents in vegetable- and oilseed-type Brassica rapa plants

Glucosinolates (GSLs) are secondary metabolites providing defense against pathogens and herbivores in plants, and anti-carcinogenic activity against human cancer cells. Profiles of GSLs vary greatly among members of genus Brassica. In this study, we found that a reference line of Chinese cabbage (B. rapa ssp. pekinensis), 'Chiifu' contains significantly lower amounts of total GSLs than the oilseed-type B. rapa (B. rapa ssp. trilocularis) line 'LP08'. This study aimed to identify the key regulators of the high accumulation of GSLs in Brassica rapa plants using transcriptomic and linkage mapping approaches. Comparative transcriptome analysis showed that, in total, 8,276 and 9,878 genes were differentially expressed between 'Chiifu' and 'LP08' under light and dark conditions, respectively. Among 162 B. rapa GSL pathway genes, 79 were related to GSL metabolism under light conditions. We also performed QTL analysis using a single nucleotide polymorphism-based linkage map constructed using 151 F₅ individuals derived from a cross between the 'Chiifu' and 'LP08' inbred lines. Two major QTL peaks were successfully identified on chromosome 3 using high-performance liquid chromatography to obtain GSL profiles from 97 F₅ recombinant inbred lines. The MYB-domain transcription factor gene BrMYB28.1 (Bra012961) was found in the highest QTL peak region. The second highest peak was located near the 2-oxoacid-dependent dioxygenase gene BrGSL-OH.1 (Bra022920). This study identified major genes responsible for differing profiles of GSLs between 'Chiifu' and 'LP08'. Thus, our study provides molecular insights into differences in GSL profiles between vegetative- and oilseed-type B. rapa plants.

BocODD1 and BocODD2 Regulate the Biosynthesis of Progoitrin Glucosinolate in Chinese Kale

Progoitrin (2-hydroxy-3-butenyl glucosinolate, PRO) is the main source of bitterness of Brassica plants. Research on the biosynthesis of PRO glucosinolate can aid the understanding of the nutritional value in Brassica plants. In this study, four ODD genes likely involved in PRO biosynthesis were cloned from Chinese kale. These four genes, designated as BocODD1-4, shared 75-82% similarities with the ODD sequence of Arabidopsis. The sequences of these four BocODDs were analyzed, and BocODD1 and BocODD2 were chosen for further study. The gene BocODD1,2 showed the highest expression levels in the roots, followed by the leaves, flowers, and stems, which is in accordance with the trend of the PRO content in the same tissues. Both the expression levels of BocODD1,2 and the content of PRO were significantly induced by high- and low-temperature treatments. The function of BocODDs involved in PRO biosynthesis was identified. Compared with the wild type, the content of PRO was increased twofold in the over-expressing BocODD1 or BocODD2 plants. Meanwhile, the content of PRO was decreased in the BocODD1 or BocODD2 RNAi lines more than twofold compared to the wildtype plants. These results suggested that BocODD1 and BocODD2 may play important roles in the biosynthesis of PRO glucosinolate in Chinese kale.

Glucosinolates and Biotic Stress Tolerance in Brassicaceae with Emphasis on Cabbage: A Review

Glucosinolates (GSLs) and GSL-associated genes are receiving increasing attention from molecular biologists due to their multifunctional properties. GSLs are secondary metabolites considered to be highly active in most Brassica species. Their importance has motivated the discovery and functional analysis of the GSLs and GSL hydrolysis products involved in disease development in brassicas and other plants. Comprehensive knowledge of the GSL content of Brassica species and the molecular details of GSL-related genes will help elucidate the molecular control of this plant defense system. This report provides an overview of the current status of knowledge on GSLs, GSL biosynthesis, as well as hydrolysis related genes, and GSL hydrolysis products that regulate fungal, bacterial, and insect resistance in cabbage and other brassicas.

Glutathione degradation activity of γ-glutamyl peptidase 1 manifests its dual roles in primary and secondary sulfur metabolism in Arabidopsis

Glutathione (GSH) functions as a major sulfur repository and hence occupies an important position in primary sulfur metabolism. GSH degradation results in sulfur reallocation and is believed to be carried out mainly by γ-glutamyl cyclotransferases (GGCT2;1, GGCT2;2, and GGCT2;3), which, however, do not fully explain the rapid GSH turnover. Here, we discovered that γ-glutamyl peptidase 1 (GGP1) contributes to GSH degradation through a yeast complementation assay. Recombinant proteins of GGP1, as well as GGP3, showed high degradation activity of GSH, but not of oxidized glutathione (GSSG), in vitro. Notably, the GGP1 transcripts were highly abundant in rosette leaves, in agreement with the ggp1 mutants constantly accumulating more GSH regardless of nutritional conditions. Given the lower energy requirements of the GGP- than the GGCT-mediated pathway, the GGP-mediated pathway could be a more efficient route for GSH degradation than the GGCT-mediated pathway. Therefore, we propose a model wherein cytosolic GSH is degraded chiefly by GGP1 and likely also by GGP3. Another noteworthy fact is that GGPs are known to process GSH conjugates in glucosinolate and camalexin synthesis; indeed, we confirmed that the ggp1 mutant contained higher levels of O-acetyl-l-Ser, a signaling molecule for sulfur starvation, and lower levels of glucosinolates and their degradation products. The predicted structure of GGP1 further provided a rationale for this hypothesis. In conclusion, we suggest that GGP1 and possibly GGP3 play vital roles in both primary and secondary sulfur metabolism.

Structural insights into partner selection for MYB and bHLH transcription factor complexes

MYB and basic helix-loop-helix (bHLH) transcription factors form complexes to regulate diverse metabolic and developmental processes in plants. However, the molecular mechanisms responsible for MYB-bHLH interaction and partner selection remain unclear. Here, we report the crystal structures of three MYB-bHLH complexes (WER-EGL3, CPC-EGL3 and MYB29-MYC3), uncovering two MYB-bHLH interaction modes. WER and CPC are R2R3- and R3-type MYBs, respectively, but interact with EGL3 through their N-terminal R3 domain in a similar mode. A single amino acid of CPC, Met49, is crucial for competition with WER to interact with EGL3. MYB29, a R2R3-type MYB transcription factor, interacts with MYC3 by its C-terminal MYC-interaction motif. The WER-EGL3 and MYB29-MYC3 binding modes are widely applied among MYB-bHLH complexes in Arabidopsis and evolve independently in plants.

An aquaporin gene MdPIP1;2 from Malus domestica confers salt tolerance in transgenic Arabidopsis

Aquaporins are known as water channel proteins. In this study, an aquaporin gene MdPIP1;2 was cloned from Malus domestica cv. Qinguan encoding a protein of 289 amino acids that formed the typical structure of aquaporin by six transmembrane domains, two asparagine-proline-alanine motifs, aromatic/arginine filter, and Forger's position. MdPIP1;2 was highly expressed in the water-sensitive or water-requiring tissues, and upregulated by salt and PEG stresses. MdPIP1;2 transgenic Arabidopsis exhibited enhanced salt stress tolerance with less Na ⁺ accumulation, lower malondialdehyde (MDA) content, lower electrolyte leakage (EL) level, and higher superoxide dismutase (SOD) and peroxidase (POD) activities compared with WT plants. Additionally, transcriptome analysis indicated MdPIP1;2 transgenic Arabidopsis could present healthier growth and development condition probably through regulating morphological structures and accumulating specific secondary metabolites under salt stress. Our results are a useful reference for better understanding the biological function of aquaporin in apple tree, especially in plant response to abiotic stress.

IQD1 Involvement in Hormonal Signaling and General Defense Responses Against Botrytis cinerea

IQ Domain 1 (IQD1) is a novel Arabidopsis thaliana calmodulin-binding protein, which was found to be a positive regulator of glucosinolate (GS) accumulation and plant defense responses against insects. We demonstrate here that the IQD1 overexpressing line (IQD1 ^OXP ) was also more resistant also to the necrotrophic fungus Botrytis cinerea, whereas an IQD1 knockout line (iqd1-1) was much more sensitive. Furthermore, we showed that IQD1 is up-regulated by jasmonic acid (JA) and downregulated by salicylic acid (SA). A comparison of whole transcriptome expression between iqd1-1 and wild type plants revealed a substantial downregulation of genes involved in plant defense and hormone regulation. Further examination revealed a marked reduction of SA and increases in the levels of ethylene, JA and abscisic acid response genes in the iqd1-1 line. Moreover, quantification of SA, JA, and abscisic acids in IQD1 ^OXP and iqd1-1 lines relative to the wild type, showed a significant reduction in endogenous JA levels in the knockout line, simultaneously with increased SA levels. Relations between IQD1 ^OXP and mutants defective in plant-hormone response indicated that IQD1 cannot rescue the absence of NPR1 or impaired SA accumulation in the NahG line. IQD1 cannot rescue ein2 or eto1 mutations connected to the ethylene pathway involved in both defense responses against B. cinerea and in regulating GS accumulation. Furthermore, IQD1cannot rescue the aos, coi1 or jar1mutations, all involved in the defense response against B. cinerea and it depends on JAR1 to control indole glucosinolate accumulation. We also found that in the B. cinerea, which infected the iqd1-1 mutant, the most abundant upregulated group of proteins is involved in the degradation of complex carbohydrates, as correlated with the sensitivity of this mutant. In summary, our results suggest that IQD1 is an important A. thaliana defensive protein against B. cinerea that is integrated into several important pathways, such as those involved in plant defense and hormone responses.

De novo transcriptome characterization of Iris atropurpurea (the Royal Iris) and phylogenetic analysis of MADS-box and R2R3-MYB gene families

The Royal Irises (section Oncocyclus) are a Middle-Eastern group of irises, characterized by extremely large flowers with a huge range of flower colors and a unique pollination system. The Royal Irises are considered to be in the course of speciation and serve as a model for evolutionary processes of speciation and pollination ecology. However, no transcriptomic and genomic data are available for these plants. Transcriptome sequencing is a valuable resource for determining the genetic basis of ecological-meaningful traits, especially in non-model organisms. Here we describe the de novo transcriptome assembly of Iris atropurpurea, an endangered species endemic to Israel's coastal plain. We sequenced and analyzed the transcriptomes of roots, leaves, and three stages of developing flower buds. To identify genes involved in developmental processes we generated phylogenetic gene trees for two major gene families, the MADS-box and MYB transcription factors, which play an important role in plant development. In addition, we identified 1503 short sequence repeats that can be developed for molecular markers for population genetics in irises. This first reported transcriptome for the Royal Irises, and the data generated, provide a valuable resource for this non-model plant that will facilitate gene discovery, functional genomic studies, and development of molecular markers in irises, to complete the intensive eco-evolutionary studies of this group.

Arabidopsis MYB28 and MYB29 transcription factors are involved in ammonium-mediated alterations of root-system architecture

Ammonium (NH₄⁺) is known to produce alterations in root-system architecture, notably, by inhibiting primary root elongation and stimulating lateral root branching. This stimulation is associated with higher auxin transport promoted by apoplast acidification. Recently, we showed that MYB28 and MYB29 transcription factors play a role in ammonium tolerance, since its double mutant (myb28myb29) is highly hypersensitive toward ammonium nutrition in relation to altered Fe homeostasis. In the present work, we observed that primary root elongation was lower in the mutant with respect to wild-type plants under ammonium nutrition. Moreover, ammonium-induced lateral root branching was impaired in myb28myb29 in a Fe-supply dependent manner. Further research is required to decipher the link between MYB28 and MYB29 functions and the signaling pathway leading to root-system architecture modification by NH₄⁺ supply.

Abscisic acid-mediated induction of FLAVIN-CONTAINING MONOOXYGENASE 2 leads to reduced accumulation of methylthioalkyl glucosinolates in Arabidopsis thaliana

Side-chain modification contributes to the structural diversity of aliphatic glucosinolates (GSLs), a class of sulfur-containing secondary metabolites found in Brassicales. The first step in side-chain modification of aliphatic GSLs is the S-oxygenation of the methylthioalkyl (MT) moiety to the methylsulfinylalkyl (MS) moiety. This reaction is catalyzed by flavin-containing monooxygenase (FMO_GS-OX), which is encoded by seven genes in Arabidopsis thaliana. Therefore, the regulation of FMO_GS-OX gene expression is key to controlling side-chain structural diversity. In this study, we demonstrated that the expression of FMO_GS-OX2 and FMO_GS-OX4 was induced by glucose treatment, independent of MYB28/29 and MYC2/3/4, the transcription factors that positively regulate aliphatic GSL biosynthesis. Glucose treatment of the abscisic acid (ABA)-related mutants indicated that glucose-triggered upregulation of FMO_GS-OX2 and FMO_GS-OX4 was partially regulated by ABA through the key negative regulators ABI1 and ABI2, and the positive regulator SnRK2, but not via the transcription factor ABI5. In wild-type plants, glucose treatment drastically reduced the accumulation of 4-methylthiobutyl (4MT) GSL, whereas a decrease in 4MT GSL was not observed in the fmo_gs-ox2, abi1-1, abi2-1, aba2-1, or aba3-1 mutants. This result indicated that the decreased accumulation of 4MT GSL by glucose treatment was attributed to upregulation of FMO_GS-OX2 via the ABA signaling pathway.

Regulation of glucosinolate biosynthesis

Glucosinolates are secondary defense metabolites produced by plants of the order Brassicales, which includes the model species Arabidopsis and many crop species. In the past 13 years, the regulation of glucosinolate synthesis in plants has been intensively studied, with recent research revealing complex molecular mechanisms that connect glucosinolate production with responses to other central pathways. In this review, we discuss how the regulation of glucosinolate biosynthesis is ecologically relevant for plants, how it is controlled by transcription factors, and how this transcriptional machinery interacts with hormonal, environmental, and epigenetic mechanisms. We present the central players in glucosinolate regulation, MYB and basic helix-loop-helix transcription factors, as well as the plant hormone jasmonate, which together with other hormones and environmental signals allow the coordinated and rapid regulation of glucosinolate genes. Furthermore, we highlight the regulatory connections between glucosinolates, auxin, and sulfur metabolism and discuss emerging insights and open questions on the regulation of glucosinolate biosynthesis.

Arabidopsis thaliana transcription factors MYB28 and MYB29 shape ammonium stress responses by regulating Fe homeostasis

Although ammonium (NH₄⁺ ) is a key intermediate of plant nitrogen metabolism, high concentrations of NH₄⁺ in the soil provoke physiological disorders that lead to the development of stress symptoms. Ammonium nutrition was shown to induce the accumulation of glucosinolates (GSLs) in leaves of different Brassicaceae species. To further understand the link between ammonium nutrition and GSLs, we analysed the ammonium stress response of Arabidopsis mutants impaired in GSL metabolic pathway. We showed that the MYB28 and MYB29 double mutant (myb28myb29), which is almost deprived of aliphatic GSLs, is highly hypersensitive to ammonium nutrition. Moreover, we evidenced that the stress symptoms developed were not a consequence of the lack of aliphatic GSLs. Transcriptomic analysis highlighted the induction of an iron (Fe) deficiency response in myb28myb29 under ammonium nutrition. Consistently, ammonium-grown myb28myb29 plants showed altered Fe accumulation and homeostasis. Interestingly, we showed overall that growing Arabidopsis with increased Fe availability relieved ammonium stress symptoms and that this was associated with MYB28 and MYB29 expression. Taken together, our data indicated that the control of Fe homeostasis was crucial for the Arabidopsis response to ammonium nutrition and evidenced that MYB28 and MYB29 play a role in this control.

Effect of selenium-sulfur interaction on the anabolism of sulforaphane in broccoli

The effects of S (as sulphate) and Se (as selenite) treatment (S mM/Se μM: 1/0, 1/50, 1/100, 1/150, 4/0, 4/50, 4/100, and 4/150) on the production of sulforaphane (an anticancer compound), the accumulation of its precursor substance, and the expression of genes related to glucoraphanin biosynthesis in broccoli were examined. Sulforaphane yield and myrosinase activity increased significantly with the combined application of 4 mM S and 100 μM Se on broccoli. Furthermore, the concentrations of glucoraphanin (a sulforaphane precursor) and methionine (a glucoraphanin substrate) slightly changed after Se application. And the strong anticancer activity of compound Se-SMC was further improved. Analysis of related gene expression showed that MY, which encodes myrosinase, was strongly induced by Se treatment. Thus, the myrosinase activity induced by Se treatment is the dominant factor affecting sulforaphane yield from glucoraphanin hydrolyzation. Selenium-sulfur biofortification provides a technical support for the cultivation of broccoli with high sulforaphane and high anti-cancer selenium compounds.

The Arabidopsis lncRNA ASCO modulates the transcriptome through interaction with splicing factors

Alternative splicing (AS) is a major source of transcriptome diversity. Long noncoding RNAs (lncRNAs) have emerged as regulators of AS through different molecular mechanisms. In Arabidopsis thaliana, the AS regulators NSRs interact with the ALTERNATIVE SPLICING COMPETITOR (ASCO) lncRNA. Here, we analyze the effect of the knock-down and overexpression of ASCO at the genome-wide level and find a large number of deregulated and differentially spliced genes related to flagellin responses and biotic stress. In agreement, ASCO-silenced plants are more sensitive to flagellin. However, only a minor subset of deregulated genes overlaps with the AS defects of the nsra/b double mutant, suggesting an alternative way of action for ASCO. Using biotin-labeled oligonucleotides for RNA-mediated ribonucleoprotein purification, we show that ASCO binds to the highly conserved spliceosome component PRP8a. ASCO overaccumulation impairs the recognition of specific flagellin-related transcripts by PRP8a. We further show that ASCO also binds to another spliceosome component, SmD1b, indicating that it interacts with multiple splicing factors. Hence, lncRNAs may integrate a dynamic network including spliceosome core proteins, to modulate transcriptome reprogramming in eukaryotes.

Function, dynamics and evolution of network motif modules in integrated gene regulatory networks of worm and plant

Gene regulatory networks (GRNs) consist of different molecular interactions that closely work together to establish proper gene expression in time and space. Especially in higher eukaryotes, many questions remain on how these interactions collectively coordinate gene regulation. We study high quality GRNs consisting of undirected protein–protein, genetic and homologous interactions, and directed protein–DNA, regulatory and miRNA–mRNA interactions in the worm Caenorhabditis elegans and the plant Arabidopsis thaliana. Our data-integration framework integrates interactions in composite network motifs, clusters these in biologically relevant, higher-order topological network motif modules, overlays these with gene expression profiles and discovers novel connections between modules and regulators. Similar modules exist in the integrated GRNs of worm and plant. We show how experimental or computational methodologies underlying a certain data type impact network topology. Through phylogenetic decomposition, we found that proteins of worm and plant tend to functionally interact with proteins of a similar age, while at the regulatory level TFs favor same age, but also older target genes. Despite some influence of the duplication mode difference, we also observe at the motif and module level for both species a preference for age homogeneity for undirected and age heterogeneity for directed interactions. This leads to a model where novel genes are added together to the GRNs in a specific biological functional context, regulated by one or more TFs that also target older genes in the GRNs. Overall, we detected topological, functional and evolutionary properties of GRNs that are potentially universal in all species.

Enhancement of Glucosinolate Production in Watercress ( Nasturtium officinale) Hairy Roots by Overexpressing Cabbage Transcription Factors

Glucosinolates are secondary metabolites that play important roles in plant defense and human health, as their production in plants is enhanced by overexpressing transcription factors. Here, four cabbage transcription factors (IQD1-1, IQD1-2, MYB29-1, and MYB29-2) affecting genes in both aliphatic and indolic glucosinolates biosynthetic pathways and increasing glucosinolates accumulation were overexpressed in watercress. Five IQD1-1, six IQD1-2, five MYB29-1, six MYB29-2, and one GUS hairy root lines were created. The expression of all genes involved in glucosinolates biosynthesis was higher in transgenic lines than in the GUS hairy root line, in agreement with total glucosinolates contents, determined by high-performance liquid chromatography. In transgenic IQD1-1 (1), IQD1-2 (4), MYB29-1 (2), and MYB29-2 (1) hairy root lines, total glucosinolates were 3.39-, 3.04-, 2.58-, and 4.69-fold higher than those in the GUS hairy root lines, respectively. These results suggest a central regulatory function for IQD1-1, IQD1-2, MYB29-1, and MYB29-2 transcription factors in glucosinolates biosynthesis in watercress hairy roots.

The Lamin-Like LITTLE NUCLEI 1 (LINC1) Regulates Pattern-Triggered Immunity and Jasmonic Acid Signaling

Pathogen-associated molecular pattern (PAMP) recognition occurs by plasma membrane located receptors that induce among other processes nuclear gene expression. However, signaling to the nuclear compartment is restricted by the nuclear envelope and nuclear pore complexes. We show here that among the four Arabidopsis lamin homologs LITTLE NUCLEI/CROWDED NUCLEI (LINC/CRWN), LINC1 plays an important role in PTI and jasmonic acid (JA) signaling. We show that linc1 knock out mutants affect PAMP-triggered MAPK activation and growth inhibition, but not reactive oxygen species or callose accumulation. We also demonstrate that linc1 mutants are compromised in regulating PAMP-triggered pathogen-related genes, in particular encoding factors involved in JA signaling and responses. Expression of a number of JAZ domain proteins, the key JA-related transcription factor MYC2 as well as key MYB transcription factors and biosynthesis genes of both the indole and aliphatic glucosinolate pathways are changed in linc1 mutants. Moreover, PAMP triggers JA and JA-Ile accumulation in linc1 mutants, whereas salicylic acid levels are unchanged. Despite impairment in PAMP-triggered immunity, linc1 mutants still show basal immunity towards Pseudomonas syringae DC3000 strains. High JA levels usually render plants resistant to necrotrophic pathogen. Thus, linc1 mutants show enhanced resistance to Botrytis cinerea infection. In accordance with a general role of LINC1 in JA signaling, linc1 mutants are hypersensitive to growth inhibition to external JA. In summary, our findings show that the lamin-like LINC1 protein plays a key role in JA signaling and regulation of PTI responses in Arabidopsis.

Changes in gene expression in Camelina sativa roots and vegetative tissues in response to salinity stress

The response of Camelina sativa to salt stress was examined. Salt reduced shoot, but not root length. Root and shoot weight were affected by salt, as was photosynthetic capacity. Salt did not alter micro-element concentration in shoots, but increased macro-element (Ca and Mg) levels. Gene expression patterns in shoots indicated that salt stress may have led to shuttling of Na+ from the cytoplasm to the tonoplast and to an increase in K+ and Ca+2 import into the cytoplasm. In roots, gene expression patterns indicated that Na+ was exported from the cytoplasm by the SOS pathway and that K+ was imported in response to salt. Genes involved in chelation and storage were up-regulated in shoots, while metal detoxification appeared to involve various export mechanisms in roots. In shoots, genes involved in secondary metabolism leading to lignin, anthocyanin and wax production were up-regulated. Partial genome partitioning was observed in roots and shoots based on the expression of homeologous genes from the three C. sativa sub-genomes. Sub-genome I and II were involved in the response to salinity stress to about the same degree, while about 10% more differentially-expressed genes were associated with sub-genome III.

Jasmonic Acid-Mediated Aliphatic Glucosinolate Metabolism Is Involved in Clubroot Disease Development in Brassica napus L

Glucosinolate (GSL) is associated with clubroot disease, which is caused by the obligate biotrophic protist Plasmodiophora brassicae. Due to the complicated composition of GSLs, their exact role in clubroot disease development remains unclear. By investigating clubroot disease resistance in cruciferous plants and characterizing the GSL content in seeds, we can determine if clubroot disease development is related to the components of GSLs. The difference in the infection process between Matthiola incana L. (resistant) and Brassica napus L. (susceptible) was determined. Root hair infection was definitely observed in both resistant and susceptible hosts, but no infection was observed during the cortical infection stage in resistant roots; this finding was verified by molecular detection of P. brassicae via PCR amplification at various times after inoculation. Based on the time course detection of the contents and compositions of GSLs after P. brassicae inoculation, susceptible roots exhibited increased accumulation of aliphatic, indolic, and aromatic GSLs in B. napus, but only aromatic GSLs were significantly increased in M. incana. Gluconapin, which was the main aliphatic GSL in B. napus and present only in B. napus, was significantly increased during the secondary infection stage. Quantification of the internal jasmonic acid (JA) concentration showed that both resistant and susceptible plants exhibited an enhanced level of JA, particularly in susceptible roots. The exogenous JA treatment induced aliphatic GSLs in B. napus and aromatic GSLs in M. incana. JA-induced aromatic GSLs may be involved in the defense against P. brassicae, whereas aliphatic GSLs induced by JA in B. napus likely play a role during the secondary infection stage. Three candidate MYB28 genes regulate the content of aliphatic GSLs identified in B. napus; one such gene was BnMYB28.1, which was significantly increased following both the treatment with exogenous JA and P. brassicae inoculation. In summary, the increased content of JA during the secondary infection stage may induce the expression of BnMYB28.1, which caused the accumulation of aliphatic GSLs in clubroot disease development.

MYB76 Inhibits Seed Fatty Acid Accumulation in Arabidopsis

The MYB family of transcription factors is important in regulatory networks controlling development, metabolism and responses to biotic and abiotic stresses in Arabidopsis. However, their role in regulating fatty acid accumulation in seeds is still largely unclear. Here, we found that MYB76, localized in the nucleus, was predominantly expressed in developing seeds during maturation. The myb76 mutation caused a significant increase in the amounts of total fatty acids and several major fatty acid compositions in mature seeds, suggesting that MYB76 functioned as an important repressor during seed oil biosynthesis. RNA sequencing and quantitative real-time PCR analysis revealed remarkable alteration of numerous genes involved in photosynthesis, fatty acid biosynthesis, modification, and degradation, and oil body formation in myb76 seeds at 12 days after pollination. These results help us to understand the novel function of MYB76 and provide new insights into the regulatory network of MYB transcriptional factors controlling seed oil accumulation in Arabidopsis.

Membrane Proteomics of Arabidopsis Glucosinolate Mutants cyp79B2/B3 and myb28/29

Glucosinolates (Gls) constitute a major group of natural metabolites represented by three major classes (aliphatic, indolic and aromatic) of more than 120 chemical structures. In our previous work, soluble proteins and metabolites in Arabidopsis mutants deficient of aliphatic (myb28/29) and indolic Gls (cyp79B2B3) were analyzed. Here we focus on investigating the changes at the level of membrane proteins in these mutants. Our LC/MS-MS analyses of tandem mass tag (TMT) labeled peptides derived from the cyp79B2/B3 and myb28/29 relative to wild type resulted in the identification of 4,673 proteins, from which 2,171 are membrane proteins. Fold changes and statistical analysis showed 64 increased and 74 decreased in cyp79B2/B3, while 28 increased and 17 decreased in myb28/29. As to the shared protein changes between the mutants, one protein was increased and eight were decreased. Bioinformatics analysis of the changed proteins led to the discovery of three cytochromes in glucosinolate molecular network (GMN): cytochrome P450 86A7 (At1g63710), cytochrome P450 71B26 (At3g26290), and probable cytochrome c (At1g22840). CYP86A7 and CYP71B26 may play a role in hydroxyl-indolic Gls production. In addition, flavone 3'-O-methyltransferase 1 represents an interesting finding as it is likely to participate in the methylation process of the hydroxyl-indolic Gls to form methoxy-indolic Gls. The analysis also revealed additional new nodes in the GMN related to stress and defense activity, transport, photosynthesis, and translation processes. Gene expression and protein levels were found to be correlated in the cyp79B2/B3, but not in the myb28/29.

Genome-wide identification and expression analysis of sulfate transporter (SULTR) genes in potato (Solanum tuberosum L.)

Solanum tuberosum genome analysis revealed 12 StSULTR genes encoding 18 transcripts. Among genes annotated at group level ( StSULTR I-IV), group III members formed the largest SULTRs-cluster and were potentially involved in biotic/abiotic stress responses via various regulatory factors, and stress and signaling proteins. Employing bioinformatics tools, this study performed genome-wide identification and expression analysis of SULTR (StSULTR) genes in potato (Solanum tuberosum L.). Very strict homology search and subsequent domain verification with Hidden Markov Model revealed 12 StSULTR genes encoding 18 transcripts. StSULTR genes were mapped on seven S. tuberosum chromosomes. Annotation of StSULTR genes was also done as StSULTR I-IV at group level based mainly on the phylogenetic distribution with Arabidopsis SULTRs. Several tandem and segmental duplications were identified between StSULTR genes. Among these duplications, Ka/Ks ratios indicated neutral nature of mutations that might not be causing any selection. Two segmental and one-tandem duplications were calculated to occur around 147.69, 180.80 and 191.00 million years ago (MYA), approximately corresponding to the time of monocot/dicot divergence. Two other segmental duplications were found to occur around 61.23 and 67.83 MYA, which is very close to the origination of monocotyledons. Most cis-regulatory elements in StSULTRs were found associated with major hormones (such as abscisic acid and methyl jasmonate), and defense and stress responsiveness. The cis-element distribution in duplicated gene pairs indicated the contribution of duplication events in conferring the neofunctionalization/s in StSULTR genes. Notably, RNAseq data analyses unveiled expression profiles of StSULTR genes under different stress conditions. In particular, expression profiles of StSULTR III members suggested their involvement in plant stress responses. Additionally, gene co-expression networks of these group members included various regulatory factors, stress and signaling proteins, and housekeeping and some other proteins with unknown functions.

Fire and Brimstone: Molecular Interactions between Sulfur and Glucosinolate Biosynthesis in Model and Crop Brassicaceae

Glucosinolates (GSLs) represent one of the most widely studied classes of plant secondary metabolite, and have a wide range of biological activities. Their unique properties also affect livestock and human health, and have been harnessed for food and other end-uses. Since GSLs are sulfur (S)-rich there are many lines of evidence suggesting that plant S status plays a key role in determining plant GSL content. However, there is still a need to establish a detailed knowledge of the distribution and remobilization of S and GSLs throughout the development of Brassica crops, and to represent this in terms of primary and secondary sources and sinks. The increased genome complexity, gene duplication and divergence within brassicas, together with their ontogenetic plasticity during crop development, appear to have a marked effect on the regulation of S and GSLs. Here, we review the current understanding of inorganic S (sulfate) assimilation into organic S forms, including GSLs and their precursors, the intracellular and inter-organ transport of inorganic and organic S forms, and the accumulation of GSLs in specific tissues. We present this in the context of overlapping sources and sinks, transport processes, signaling molecules and their associated molecular interactions. Our analysis builds on recent insights into the molecular regulation of sulfate uptake and transport by different transporters, transcription factors and miRNAs, and the role that these may play in GSL biosynthesis. We develop a provisional model describing the key processes that could be targeted in crop breeding programs focused on modifying GSL content.

Sulfur deficiency-induced repressor proteins optimize glucosinolate biosynthesis in plants

Glucosinolates (GSLs) in the plant order of the Brassicales are sulfur-rich secondary metabolites that harbor antipathogenic and antiherbivory plant-protective functions and have medicinal properties, such as carcinopreventive and antibiotic activities. Plants repress GSL biosynthesis upon sulfur deficiency (-S); hence, field performance and medicinal quality are impaired by inadequate sulfate supply. The molecular mechanism that links -S to GSL biosynthesis has remained understudied. We report here the identification of the -S marker genes sulfur deficiency induced 1 (SDI1) and SDI2 acting as major repressors controlling GSL biosynthesis in Arabidopsis under -S condition. SDI1 and SDI2 expression negatively correlated with GSL biosynthesis in both transcript and metabolite levels. Principal components analysis of transcriptome data indicated that SDI1 regulates aliphatic GSL biosynthesis as part of -S response. SDI1 was localized to the nucleus and interacted with MYB28, a major transcription factor that promotes aliphatic GSL biosynthesis, in both yeast and plant cells. SDI1 inhibited the transcription of aliphatic GSL biosynthetic genes by maintaining the DNA binding composition in the form of an SDI1-MYB28 complex, leading to down-regulation of GSL biosynthesis and prioritization of sulfate usage for primary metabolites under sulfur-deprived conditions.

New nodes and edges in the glucosinolate molecular network revealed by proteomics and metabolomics of Arabidopsis myb28/29 and cyp79B2/B3 glucosinolate mutants

Glucosinolates present in Brassicales are important for human health and plant defense against insects and pathogens. Here we investigate the proteomes and metabolomes of Arabidopsis myb28/29 and cyp79B2/B3 mutants deficient in aliphatic glucosinolates and indolic glucosinolates, respectively. Quantitative proteomics of the myb28/29 and cyp79B2/B3 mutants led to the identification of 2785 proteins, of which 142 proteins showed significant changes in the two mutants compared to wild type (WT). By mapping the differential proteins using STRING, we detected 59 new edges in the glucosinolate metabolic network. These connections can be classified as primary with direct roles in glucosinolate metabolism, secondary related to plant stress responses, and tertiary involved in other biological processes. Gene Ontology analysis of the differential proteins showed high level of enrichment in the nodes belonging to metabolic process including glucosinolate biosynthesis and response to stimulus. Using metabolomics, we quantified 292 metabolites covering a broad spectrum of metabolic pathways, and 89 exhibited differential accumulation patterns between the mutants and WT. The changing metabolites (e.g., γ-glutamyl amino acids, auxins and glucosinolate hydrolysis products) complement our proteomics findings. This study contributes toward engineering and breeding of glucosinolate profiles in plants in efforts to improve human health, crop quality and productivity.

BIOLOGICAL SIGNIFICANCE

Glucosinolates in Brassicales constitute an important group of natural metabolites important for plant defense and human health. Its biosynthetic pathways and transcriptional regulation have been well-studied. Using Arabidopsis mutants of important genes in glucosinolate biosynthesis, quantitative proteomics and metabolomics led to identification of many proteins and metabolites that are potentially related to glucosinolate metabolism. This study provides a comprehensive insight into the molecular networks of glucosinolate metabolism, and will facilitate efforts toward engineering and breeding of glucosinolate profiles for enhanced crop defense, and nutritional value.

Jasmonate-responsive transcription factors regulating plant secondary metabolism

Plants produce a large variety of secondary metabolites including alkaloids, glucosinolates, terpenoids and phenylpropanoids. These compounds play key roles in plant-environment interactions and many of them have pharmacological activity in humans. Jasmonates (JAs) are plant hormones which induce biosynthesis of many secondary metabolites. JAs-responsive transcription factors (TFs) that regulate the JAs-induced accumulation of secondary metabolites belong to different families including AP2/ERF, bHLH, MYB and WRKY. Here, we give an overview of the types and functions of TFs that have been identified in JAs-induced secondary metabolite biosynthesis, and highlight their similarities and differences in regulating various biosynthetic pathways. We review major recent developments regarding JAs-responsive TFs mediating secondary metabolite biosynthesis, and provide suggestions for further studies.

The Evolutionary History of R2R3-MYB Proteins Across 50 Eukaryotes: New Insights Into Subfamily Classification and Expansion

R2R3-MYB proteins (2R-MYBs) are one of the main transcription factor families in higher plants. Since the evolutionary history of this gene family across the eukaryotic kingdom remains unknown, we performed a comparative analysis of 2R-MYBs from 50 major eukaryotic lineages, with particular emphasis on land plants. A total of 1548 candidates were identified among diverse taxonomic groups, which allowed for an updated classification of 73 highly conserved subfamilies, including many newly identified subfamilies. Our results revealed that the protein architectures, intron patterns, and sequence characteristics were remarkably conserved in each subfamily. At least four subfamilies were derived from early land plants, 10 evolved from spermatophytes, and 19 from angiosperms, demonstrating the diversity and preferential expansion of this gene family in land plants. Moreover, we determined that their remarkable expansion was mainly attributed to whole genome and segmental duplication, where duplicates were preferentially retained within certain subfamilies that shared three homologous intron patterns (a, b, and c) even though up to 12 types of patterns existed. Through our integrated distributions, sequence characteristics, and phylogenetic tree analyses, we confirm that 2R-MYBs are old and postulate that 3R-MYBs may be evolutionarily derived from 2R-MYBs via intragenic domain duplication.

Distinct patterns of the histone marks associated with recruitment of the methionine chain-elongation pathway from leucine biosynthesis

Aliphatic glucosinolates (GLSs) are derived from chain-elongated methionine produced by an iterative three-step process, known to be evolutionarily recruited from leucine biosynthesis. The divergence of homologous genes between two pathways is mainly linked to the alterations in biochemical features. In this study, it was discovered that a distinct pattern of histone modifications is associated with and/or contributes to the divergence of the two pathways. In general, genes involved in leucine biosynthesis were robustly associated with H3k4me2 and H3K4me3. In contrast, despite the considerable abundances of H3K4me2 observed in some of genes involved in methionine chain elongation, H3K4me3 was completely missing. This H3K4m3-depleted pattern had no effect on gene transcription, whereas it seemingly co-evolved with the entire pathway of aliphatic GLS biosynthesis. The results reveal a novel association of the epigenetic marks with plant secondary metabolism, and may help to understand the recruitment of the methionine chain-elongation pathway from leucine biosynthesis.

Associative transcriptomics study dissects the genetic architecture of seed glucosinolate content in Brassica napus

Breeding new varieties with low seed glucosinolate (GS) concentrations has long been a prime target in Brassica napus. In this study, a novel association mapping methodology termed 'associative transcriptomics' (AT) was applied to a panel of 101 B. napus lines to define genetic regions and also candidate genes controlling total seed GS contents. Over 100,000 informative single-nucleotide polymorphisms (SNPs) and gene expression markers (GEMs) were developed for AT analysis, which led to the identification of 10 SNP and 7 GEM association peaks. Within these peaks, 26 genes were inferred to be involved in GS biosynthesis. A weighted gene co-expression network analysis provided additional 40 candidate genes. The transcript abundance in leaves of two candidate genes, BnaA.GTR2a located on chromosome A2 and BnaC.HAG3b on C9, was correlated with seed GS content, explaining 18.8 and 16.8% of phenotypic variation, respectively. Resequencing of genomic regions revealed six new SNPs in BnaA.GTR2a and four insertions or deletions in BnaC.HAG3b. These deletion polymorphisms were then successfully converted into polymerase chain reaction-based diagnostic markers that can, due to high linkage disequilibrium observed in these regions of the genome, be used for marker-assisted breeding for low seed GS lines.

The multi-protein family of sulfotransferases in plants: composition, occurrence, substrate specificity, and functions

All members of the sulfotransferase (SOT, EC 2.8.2.-) protein family transfer a sulfuryl group from the donor 3'-phosphoadenosine 5'-phosphosulfate (PAPS) to an appropriate hydroxyl group of several classes of substrates. The primary structure of these enzymes is characterized by a histidine residue in the active site, defined PAPS binding sites and a longer SOT domain. Proteins with this SOT domain occur in all organisms from all three domains, usually as a multi-protein family. Arabidopsis thaliana SOTs, the best characterized SOT multi-protein family, contains 21 members. The substrates for several plant enzymes have already been identified, such as glucosinolates, brassinosteroids, jasmonates, flavonoids, and salicylic acid. Much information has been gathered on desulfo-glucosinolate (dsGl) SOTs in A. thaliana. The three cytosolic dsGl SOTs show slightly different expression patterns. The recombinant proteins reveal differences in their affinity to indolic and aliphatic dsGls. Also the respective recombinant dsGl SOTs from different A. thaliana ecotypes differ in their kinetic properties. However, determinants of substrate specificity and the exact reaction mechanism still need to be clarified. Probably, the three-dimensional structures of more plant proteins need to be solved to analyze the mode of action and the responsible amino acids for substrate binding. In addition to A. thaliana, more plant species from several families need to be investigated to fully elucidate the diversity of sulfated molecules and the way of biosynthesis catalyzed by SOT enzymes.

Will chemical defenses become more effective against specialist herbivores under elevated CO2?

Elevated atmospheric CO2 is known to affect plant-insect herbivore interactions. Elevated CO2 causes leaf nitrogen to decrease, the ostensible cause of herbivore compensatory feeding. CO2 may also affect herbivore consumption by altering chemical defenses via changes in plant hormones. We considered the effects of elevated CO2, in conjunction with soil fertility and damage (simulated herbivory), on glucosinolate concentrations of mustard (Brassica nigra) and collard (B. oleracea var. acephala) and the effects of leaf nitrogen and glucosinolate groups on specialist Pieris rapae consumption. Elevated CO2 affected B. oleracea but not B. nigra glucosinolates; responses to soil fertility and damage were also species-specific. Soil fertility and damage also affected B. oleracea glucosinolates differently under elevated CO2. Glucosinolates did not affect P. rapae consumption at either CO2 concentration in B. nigra, but had CO2-specific effects on consumption in B. oleracea. At ambient CO2, leaf nitrogen had strong effects on glucosinolate concentrations and P. rapae consumption but only gluconasturtiin was a feeding stimulant. At elevated CO2, direct effects of leaf nitrogen were weaker, but glucosinolates had stronger effects on consumption. Gluconasturtiin and aliphatic glucosinolates were feeding stimulants and indole glucosinolates were feeding deterrents. These results do not support the compensatory feeding hypothesis as the sole driver of changes in P. rapae consumption under elevated CO2. Support for hormone-mediated CO2 response (HMCR) was mixed; it explained few treatment effects on constitutive or induced glucosinolates, but did explain patterns in SEMs. Further, the novel feeding deterrent effect of indole glucosinolates under elevated CO2 in B. oleracae underscores the importance of defensive chemistry in CO2 response. We speculate that P. rapae indole glucosinolate detoxification mechanisms may have been overwhelmed under elevated CO2 forcing slowed consumption. Specialists may have to contend with hosts with poorer nutritional quality and more effective chemical defenses under elevated CO2.

Update on the role of R2R3-MYBs in the regulation of glucosinolates upon sulfur deficiency

To balance the flux of sulfur (S) into glucosinolates (GSL) and primary metabolites plants exploit various regulatory mechanisms particularly important upon S deficiency (-S). The role of MYB34, MYB51 and MYB122 controlling the production of indolic glucosinolates (IGs) and MYB28, MYB29, and MYB76 regulating the biosynthesis of aliphatic glucosinolates (AGs) in Arabidopsis thaliana has not been fully addressed at -S conditions yet. We show that the decline in the concentrations of GSL during S depletion does not coincide with the globally decreased transcription of R2R3-MYBs. Whereas the levels of GSL are diminished, the expression of MYB34, MYB51, MYB122, and MYB28 is hardly changed in early phase of S limitation. Furthermore, the mRNA levels of these MYBs start to raise under prolonged S starvation. In parallel, we found that SLIM1 can downregulate the MYBs in vitro as demonstrated in trans-activation assays in cultured Arabidopsis cells with SLIM1 as effector and ProMYB51:uidA as a reporter construct. However, in vivo, only the mRNA of MYB29 and MYB76 correlated with the levels of GSL at -S. We propose that the negative effect of SLIM1 on GSL regulatory genes can be overridden by a "low GSL signal" inducing the transcription of MYBs in a feedback regulatory loop. In accordance with this hypothesis, the expression of MYB34, MYB51, MYB122, and CYP83B1 was further induced in cyp79b2 cyp79b3 mutant exposed to -S conditions vs. cyp79b2 cyp79b3 plants grown on control medium. In addition, the possible role of MYBs in the regulation of essential S assimilation enzymes, in the regulation of GSL biosynthesis upon accelerated termination of life cycles, in the mobilization of auxin and lateral root formation at S deficiency is discussed.

Novel bioresources for studies of Brassica oleracea: identification of a kale MYB transcription factor responsible for glucosinolate production

Plants belonging to the Brassicaceae family exhibit species-specific profiles of glucosinolates (GSLs), a class of defence compounds against pathogens and insects. GSLs also exhibit various human health-promoting properties. Among them, glucoraphanin (aliphatic 4-methylsulphinylbutyl GSL) has attracted the most attention because it hydrolyses to form a potent anticancer compound. Increased interest in developing commercial varieties of Brassicaceae crops with desirable GSL profiles has led to attempts to identify genes that are potentially valuable for controlling GSL biosynthesis. However, little attention has been focused on genes of kale (Brassica oleracea var. acephala). In this study, we established full-length kale cDNA libraries containing 59 904 clones, which were used to generate an expressed sequence tag (EST) data set with 119 204 entries. The EST data set clarified genes related to the GSL biosynthesis pathway in kale. We specifically focused on BoMYB29, a homolog of Arabidopsis MYB29/PMG2/HAG3, not only to characterize its function but also to demonstrate its usability as a biological resource. BoMYB29 overexpression in wild-type Arabidopsis enhanced the expression of aliphatic GSL biosynthetic genes and the accumulation of aliphatic GSLs. When expressed in the myb28myb29 mutant, which exhibited no detectable aliphatic GSLs, BoMYB29 restored the expression of biosynthetic genes and aliphatic GSL accumulation. Interestingly, the ratio of methylsulphinyl GSL content, including glucoraphanin, to that of methylthio GSLs was greatly increased, indicating the suitability of BoMYB29 as a regulator for increasing methylsulphinyl GSL content. Our results indicate that these biological resources can facilitate further identification of genes useful for modifications of GSL profiles and accumulation in kale.