Cistrome and Epicistrome Features Shape the Regulatory DNA Landscape.

The cistrome is the complete set of transcription factor (TF) binding sites (cis-elements) in an organism, while an epicistrome incorporates tissue-specific DNA chemical modifications and TF-specific chemical sensitivities into these binding profiles. Robust methods to construct comprehensive cistrome and epicistrome maps are critical for elucidating complex transcriptional networks that underlie growth, behavior, and disease. Here, we describe DNA affinity purification sequencing (DAP-seq), a high-throughput TF binding site discovery method that interrogates genomic DNA with in-vitro-expressed TFs. Using DAP-seq, we defined the Arabidopsis cistrome by resolving motifs and peaks for 529 TFs. Because genomic DNA used in DAP-seq retains 5-methylcytosines, we determined that >75% (248/327) of Arabidopsis TFs surveyed were methylation sensitive, a property that strongly impacts the epicistrome landscape. DAP-seq datasets also yielded insight into the biology and binding site architecture of numerous TFs, demonstrating the value of DAP-seq for cost-effective cistromic and epicistromic annotation in any organism.

Determination and inference of eukaryotic transcription factor sequence specificity.

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

The retrograde signaling regulator ANAC017 recruits the MKK9-MPK3/6, ethylene, and auxin signaling pathways to balance mitochondrial dysfunction with growth.

In plant cells, mitochondria are ideally positioned to sense and balance changes in energy metabolism in response to changing environmental conditions. Retrograde signaling from mitochondria to the nucleus is crucial for adjusting the required transcriptional responses. We show that ANAC017, the master regulator of mitochondrial stress, directly recruits a signaling cascade involving the plant hormones ethylene and auxin as well as the MAP KINASE KINASE (MKK) 9-MAP KINASE (MPK) 3/6 pathway in Arabidopsis thaliana. Chromatin immunoprecipitation followed by sequencing and overexpression demonstrated that ANAC017 directly regulates several genes of the ethylene and auxin pathways, including MKK9, 1-AMINO-CYCLOPROPANE-1-CARBOXYLATE SYNTHASE 2, and YUCCA 5, in addition to genes encoding transcription factors regulating plant growth and stress responses such as BASIC REGION/LEUCINE ZIPPER MOTIF (bZIP) 60, bZIP53, ANAC081/ATAF2, and RADICAL-INDUCED CELL DEATH1. A time-resolved RNA-seq experiment established that ethylene signaling precedes the stimulation of auxin signaling in the mitochondrial stress response, with a large part of the transcriptional regulation dependent on ETHYLENE-INSENSITIVE 3. These results were confirmed by mutant analyses. Our findings identify the molecular components controlled by ANAC017, which integrates the primary stress responses to mitochondrial dysfunction with whole plant growth via the activation of regulatory and partly antagonistic feedback loops.

Spatially resolved transcriptomic analysis of the germinating barley grain.

Seeds are a vital source of calories for humans and a unique stage in the life cycle of flowering plants. During seed germination, the embryo undergoes major developmental transitions to become a seedling. Studying gene expression in individual seed cell types has been challenging due to the lack of spatial information or low throughput of existing methods. To overcome these limitations, a spatial transcriptomics workflow was developed for germinating barley grain. This approach enabled high-throughput analysis of spatial gene expression, revealing specific spatial expression patterns of various functional gene categories at a sub-tissue level. This study revealed over 14 000 genes differentially regulated during the first 24 h after imbibition. Individual genes, such as the aquaporin gene family, starch degradation, cell wall modification, transport processes, ribosomal proteins and transcription factors, were found to have specific spatial expression patterns over time. Using spatial autocorrelation algorithms, we identified auxin transport genes that had increasingly focused expression within subdomains of the embryo over time, suggesting their role in establishing the embryo axis. Overall, our study provides an unprecedented spatially resolved cellular map for barley germination and identifies specific functional genomics targets to better understand cellular restricted processes during germination. The data can be viewed at https://spatial.latrobe.edu.au/.

Spatiotemporal deposition of cell wall polysaccharides in oat endosperm during grain development.

Oat (Avena sativa) is a cereal crop whose grains are rich in (1,3;1,4)-ß-D-glucan (mixed-linkage glucan or MLG), a soluble dietary fiber. In our study, we analyzed oat endosperm development in 2 Canadian varieties with differing MLG content and nutritional value. We confirmed that oat undergoes a nuclear type of endosperm development but with a shorter cellularization phase than barley (Hordeum vulgare). Callose and cellulose were the first polysaccharides to be detected in the early anticlinal cell walls at 11 days postemergence (DPE) of the panicle. Other polysaccharides such as heteromannan and homogalacturonan were deposited early in cellularization around 12 DPE after the first periclinal walls are laid down. In contrast to barley, heteroxylan deposition coincided with completion of cellularization and was detected from 14 DPE but was only detectable after demasking. Notably, MLG was the last polysaccharide to be laid down at 18 DPE within the differentiation phase, rather than during cellularization. In addition, differences in the spatiotemporal patterning of MLG were also observed between the 2 varieties. The lower MLG-containing cultivar AC Morgan (3.5% w/w groats) was marked by the presence of a discontinuous pattern of MLG labeling, while labeling in the same walls in CDC Morrison (5.6% w/w groats) was mostly even and continuous. RNA-sequencing analysis revealed higher transcript levels of multiple MLG biosynthetic cellulose synthase-like F (CSLF) and CSLH genes during grain development in CDC Morrison compared with AC Morgan that likely contributes to the increased abundance of MLG at maturity in CDC Morrison. CDC Morrison was also observed to have smaller endosperm cells with thicker walls than AC Morgan from cellularization onwards, suggesting the processes controlling cell size and shape are established early in development. This study has highlighted that the molecular processes influencing MLG content and deposition are more complex than previously imagined.

RNA-seq analysis of laser microdissected Arabidopsis thaliana leaf epidermis, mesophyll and vasculature defines tissue-specific transcriptional responses to multiple stress treatments.

Acclimation of plants to adverse conditions requires the coordination of gene expression and signalling pathways between tissues and cell types. As the energy and carbon capturing organs, leaves are significantly affected by abiotic and biotic stresses. However, tissue- or cell type-specific analyses of stress responses have focussed on the Arabidopsis root. Here, we comparatively explore the transcriptomes of three leaf tissues (epidermis, mesophyll, vasculature) after induction of diverse stress pathways by chemical stimuli (antimycin A, 3-amino-1,2,4-triazole, methyl viologen, salicylic acid) and ultraviolet light in Arabidopsis using laser capture microdissection followed by RNA sequencing. Stimulation of stress pathways caused an overall reduction in the number of genes expressed in a tissue-specific manner, though a small subset gained or changed their tissue specificity. We find no evidence of a common stress response, with only a few genes consistently responsive to two or more treatments in the analysed tissues. However, differentially expressed genes overlap between tissues for individual treatments. A focussed analysis provided evidence for an interaction of auxin and ethylene that mediates retrograde signalling during mitochondrial dysfunction specifically in the epidermis, and a gene regulatory network defined the hierarchy of interactions. Taken together, we have generated an extensive reference dataset that will be valuable for future experiments analysing transcriptional responses on a tissue or single-cell level. Our results will enable the tailoring of the tissue-specific engineering of stress-tolerant plants.

A MYC2/MYC3/MYC4-dependent transcription factor network regulates water spray-responsive gene expression and jasmonate levels.

Mechanical stimuli, such as wind, rain, and touch affect plant development, growth, pest resistance, and ultimately reproductive success. Using water spray to simulate rain, we demonstrate that jasmonic acid (JA) signaling plays a key role in early gene-expression changes, well before it leads to developmental changes in flowering and plant architecture. The JA-activated transcription factors MYC2/MYC3/MYC4 modulate transiently induced expression of 266 genes, most of which peak within 30 min, and control 52% of genes induced >100-fold. Chromatin immunoprecipitation-sequencing analysis indicates that MYC2 dynamically binds >1,300 promoters and trans-activation assays show that MYC2 activates these promoters. By mining our multiomic datasets, we identified a core MYC2/MYC3/MYC4-dependent "regulon" of 82 genes containing many previously unknown MYC2 targets, including transcription factors bHLH19 and ERF109 bHLH19 can in turn directly activate the ORA47 promoter, indicating that MYC2/MYC3/MYC4 initiate a hierarchical network of downstream transcription factors. Finally, we also reveal that rapid water spray-induced accumulation of JA and JA-isoleucine is directly controlled by MYC2/MYC3/MYC4 through a positive amplification loop that regulates JA-biosynthesis genes.

Temporal tissue-specific regulation of transcriptomes during barley (Hordeum vulgare) seed germination.

The distinct functions of individual cell types require cells to express specific sets of genes. The germinating seed is an excellent model to study genome regulation between cell types since the majority of the transcriptome is differentially expressed in a short period, beginning from a uniform, metabolically inactive state. In this study, we applied laser-capture microdissection RNA-sequencing to small numbers of cells from the plumule, radicle tip and scutellum of germinating barley seeds every 8 h, over a 48 h time course. Tissue-specific gene expression was notably common; 25% (910) of differentially expressed transcripts in plumule, 34% (1876) in radicle tip and 41% (2562) in scutellum were exclusive to that organ. We also determined that tissue-specific storage of transcripts occurs during seed development and maturation. Co-expression of genes had strong spatiotemporal structure, with most co-expression occurring within one organ and at a subset of specific time points during germination. Overlapping and distinct enrichment of functional categories were observed in the tissue-specific profiles. We identified candidate transcription factors amongst these that may be regulators of spatiotemporal gene expression programs. Our findings contribute to the broader goal of generating an integrative model that describes the structure and function of individual cells within seeds during germination.

The JA-pathway MYC transcription factors regulate photomorphogenic responses by targeting HY5 gene expression.

Jasmonates are key regulators of the balance between defence and growth in plants. However, the molecular mechanisms by which activation of defence reduces growth are not yet fully understood. Here, we analyze the role of MYC transcription factors (TFs) and jasmonic acid (JA) in photomorphogenic growth. We found that multiple myc mutants share light-associated phenotypes with mutants of the phytochrome B photoreceptor, such as delayed seed germination in the dark and long hypocotyl growth. Overexpression of MYC2 in a phyB background partially suppressed its long hypocotyl phenotype. Transcriptomic analysis of multiple myc mutants confirmed that MYCs are required for full expression of red (R) light-regulated genes, including the master regulator HY5. ChIP-seq analyses revealed that MYC2 and MYC3 bind directly to the promoter of HY5 and that HY5 gene expression and protein levels are compromised in multiple myc mutants. Altogether, our results pinpoint MYCs as photomorphogenic TFs that control phytochrome responses by activating HY5 expression. This has important implications in understanding the trade-off between growth and defence as the same TFs that activate defence responses are photomorphogenic growth regulators.

Recent advances in Cannabis sativa genomics research.

Cannabis (Cannabis sativa L.) is one of the oldest cultivated plants purported to have unique medicinal properties. However, scientific research of cannabis has been restricted by the Single Convention on Narcotic Drugs of 1961, an international treaty that prohibits the production and supply of narcotic drugs except under license. Legislation governing cannabis cultivation for research, medicinal and even recreational purposes has been relaxed recently in certain jurisdictions. As a result, there is now potential to accelerate cultivar development of this multi-use and potentially medically useful plant species by application of modern genomics technologies. Whilst genomics has been pivotal to our understanding of the basic biology and molecular mechanisms controlling key traits in several crop species, much work is needed for cannabis. In this review we provide a comprehensive summary of key cannabis genomics resources and their applications. We also discuss prospective applications of existing and emerging genomics technologies for accelerating the genetic improvement of cannabis.

ANAC017 Coordinates Organellar Functions and Stress Responses by Reprogramming Retrograde Signaling.

Mitochondria adjust their activities in response to external and internal stimuli to optimize growth via the mitochondrial retrograde response signaling pathway. The Arabidopsis (Arabidopsis thaliana) NAC domain transcription factor ANAC017 has previously been identified as a regulator of the mitochondrial retrograde response. We show here that overexpression of ANAC017 in Arabidopsis leads to growth retardation, altered leaf development with decreased cell size and viability, and early leaf senescence. RNA sequencing analyses revealed that increased ANAC017 expression leads to higher expression of genes related to mitochondrial stress, cell death/autophagy, and leaf senescence under nonlimiting growth conditions as well as extensive repression of chloroplast function. Gene regulatory network analysis indicated that a complex hierarchy of transcription factors exists downstream of ANAC017. These involve a set of up-regulated ANAC and WRKY transcription factors associated with organellar signaling and senescence. The network also includes a number of ethylene- and gibberellic acid-related transcription factors with established functions in stress responses and growth regulation, which down-regulate their target genes. A number of BASIC LEUCINE-ZIPPER MOTIF transcription factors involved in the endoplasmic reticulum unfolded protein response or balancing of energy homeostasis via the SNF1-RELATED PROTEIN KINASE1 were also down-regulated by ANAC017 overexpression. Our results show that the endoplasmic reticulum membrane tethering of the constitutively expressed ANAC017, and its controlled release, are crucial to fine-tune a fast reactive but potentially harmful signaling cascade. Thus, ANAC017 is a master regulator of cellular responses with mitochondria acting as central sensors.

Analysis of Spatio-Temporal Transcriptome Profiles of Soybean (Glycine max) Tissues during Early Seed Development.

Soybean (Glycine max) is an important crop providing oil and protein for both human and animal consumption. Knowing which biological processes take place in specific tissues in a temporal manner will enable directed breeding or synthetic approaches to improve seed quantity and quality. We analyzed a genome-wide transcriptome dataset from embryo, endosperm, endothelium, epidermis, hilum, outer and inner integument and suspensor at the global, heart and cotyledon stages of soybean seed development. The tissue specificity of gene expression was greater than stage specificity, and only three genes were differentially expressed in all seed tissues. Tissues had both unique and shared enriched functional categories of tissue-specifically expressed genes associated with them. Strong spatio-temporal correlation in gene expression was identified using weighted gene co-expression network analysis, with the most co-expression occurring in one seed tissue. Transcription factors with distinct spatiotemporal gene expression programs in each seed tissue were identified as candidate regulators of expression within those tissues. Gene ontology (GO) enrichment of orthogroup clusters revealed the conserved functions and unique roles of orthogroups with similar and contrasting expression patterns in transcript abundance between soybean and Arabidopsis during embryo proper and endosperm development. Key regulators in each seed tissue and hub genes connecting those networks were characterized by constructing gene regulatory networks. Our findings provide an important resource for describing the structure and function of individual soybean seed compartments during early seed development.

Mobile small RNAs regulate genome-wide DNA methylation.

RNA silencing at the transcriptional and posttranscriptional levels regulates endogenous gene expression, controls invading transposable elements (TEs), and protects the cell against viruses. Key components of the mechanism are small RNAs (sRNAs) of 21-24 nt that guide the silencing machinery to their nucleic acid targets in a nucleotide sequence-specific manner. Transcriptional gene silencing is associated with 24-nt sRNAs and RNA-directed DNA methylation (RdDM) at cytosine residues in three DNA sequence contexts (CG, CHG, and CHH). We previously demonstrated that 24-nt sRNAs are mobile from shoot to root in Arabidopsis thaliana and confirmed that they mediate DNA methylation at three sites in recipient cells. In this study, we extend this finding by demonstrating that RdDM of thousands of loci in root tissues is dependent upon mobile sRNAs from the shoot and that mobile sRNA-dependent DNA methylation occurs predominantly in non-CG contexts. Mobile sRNA-dependent non-CG methylation is largely dependent on the DOMAINS REARRANGED METHYLTRANSFERASES 1/2 (DRM1/DRM2) RdDM pathway but is independent of the CHROMOMETHYLASE (CMT)2/3 DNA methyltransferases. Specific superfamilies of TEs, including those typically found in gene-rich euchromatic regions, lose DNA methylation in a mutant lacking 22- to 24-nt sRNAs (dicer-like 2, 3, 4 triple mutant). Transcriptome analyses identified a small number of genes whose expression in roots is associated with mobile sRNAs and connected to DNA methylation directly or indirectly. Finally, we demonstrate that sRNAs from shoots of one accession move across a graft union and target DNA methylation de novo at normally unmethylated sites in the genomes of root cells from a different accession.

AgriSeqDB: an online RNA-Seq database for functional studies of agriculturally relevant plant species.

BACKGROUND: The genome-wide expression profile of genes in different tissues/cell types and developmental stages is a vital component of many functional genomic studies. Transcriptome data obtained by RNA-sequencing (RNA-Seq) is often deposited in public databases that are made available via data portals. Data visualization is one of the first steps in assessment and hypothesis generation. However, these databases do not typically include visualization tools and establishing one is not trivial for users who are not computational experts. This, as well as the various formats in which data is commonly deposited, makes the processes of data access, sharing and utility more difficult. Our goal was to provide a simple and user-friendly repository that meets these needs for data-sets from major agricultural crops. DESCRIPTION: AgriSeqDB ( https://expression.latrobe.edu.au/agriseqdb ) is a database for viewing, analysing and interpreting developmental and tissue/cell-specific transcriptome data from several species, including major agricultural crops such as wheat, rice, maize, barley and tomato. The disparate manner in which public transcriptome data is often warehoused and the challenge of visualizing raw data are both major hurdles to data reuse. The popular eFP browser does an excellent job of presenting transcriptome data in an easily interpretable view, but previous implementation has been mostly on a case-by-case basis. Here we present an integrated visualisation database of transcriptome data-sets from six species that did not previously have public-facing visualisations. We combine the eFP browser, for gene-by-gene investigation, with the Degust browser, which enables visualisation of all transcripts across multiple samples. The two visualisation interfaces launch from the same point, enabling users to easily switch between analysis modes. The tools allow users, even those without bioinformatics expertise, to mine into data-sets and understand the behaviour of transcripts of interest across samples and time. We have also incorporated an additional graphic download option to simplify incorporation into presentations or publications. CONCLUSION: Powered by eFP and Degust browsers, AgriSeqDB is a quick and easy-to-use platform for data analysis and visualization in five crops and Arabidopsis. Furthermore, it provides a tool that makes it easy for researchers to share their data-sets, promoting research collaborations and data-set reuse.

Regulation of genome-wide DNA methylation by mobile small RNAs.

Contents Summary 540 I. Introduction 540 II. There are different types of sRNA mobility 541 III. Mechanisms of sRNA movement 541 IV. Long-distance, shoot-root, mobile siRNAs influence DNA methylation in recipient tissues 541 V. Classes of interactions between shoot-root mobile siRNAs and DNA methylation 542 VI. Loci targeted directly and indirectly by shoot-root mobile siRNAs are associated with different histone modifications 543 VII. Is mobile siRNA-regulated DNA methylation important in specific tissues or under specific conditions? 543 VIII. Mobile sRNAs can be used to modify plant traits 544 IX. Conclusions 544 Acknowledgements 544 References 544 SUMMARY: RNA-directed DNA methylation (RdDM) at cytosine residues regulates gene expression, silences transposable elements and influences genome stability. The mechanisms responsible for RdDM are guided to target loci by small RNAs (sRNAs) that can move within plants cell to cell and long distance. Here we discuss recent advances in the understanding of interactions between mobile sRNAs and DNA methylation. We describe the mechanisms of sRNA movement, the differences between known classes of mobile sRNA-DNA methylation interactions and the limits of current knowledge. Finally, we discuss potential applications of mobile sRNAs in modifying plant traits.

JAZ2 controls stomata dynamics during bacterial invasion.

Coronatine (COR) facilitates entry of bacteria into the plant apoplast by stimulating stomata opening. COR-induced signaling events at stomata remain unclear. We found that the COR and jasmonate isoleucine (JA-Ile) co-receptor JAZ2 is constitutively expressed in guard cells and modulates stomatal dynamics during bacterial invasion We analyzed tissue expression patterns of AtJAZ genes and measured stomata opening and pathogen resistance in loss- and gain-of-function mutants. Arabidopsis jaz2 mutants are partially impaired in pathogen-induced stomatal closing and more susceptible to Pseudomonas. Gain-of-function mutations in JAZ2 prevent stomatal reopening by COR and are highly resistant to bacterial penetration. The JAZ2 targets MYC2, MYC3 and MYC4 directly regulate the expression of ANAC19, ANAC55 and ANAC72 to modulate stomata aperture. Due to the antagonistic interactions between the salicylic acid (SA) and JA defense pathways, efforts to increase resistance to biotrophs result in enhanced susceptibility to necrotrophs, and vice versa. Remarkably, dominant jaz2Δjas mutants are resistant to Pseudomonas syringae but retain unaltered resistance against necrotrophs. Our results demonstrate the existence of a COI1-JAZ2-MYC2,3,4-ANAC19,55,72 module responsible for the regulation of stomatal aperture that is hijacked by bacterial COR to promote infection. They also provide novel strategies for crop protection against biotrophs without compromising resistance to necrotrophs.

Salicylic acid treatment and expression of an RNA-dependent RNA polymerase 1 transgene inhibit lethal symptoms and meristem invasion during tobacco mosaic virus infection in Nicotiana benthamiana.

BACKGROUND: Host RNA-dependent RNA polymerases (RDRs) 1 and 6 contribute to antiviral RNA silencing in plants. RDR6 is constitutively expressed and was previously shown to limit invasion of Nicotiana benthamiana meristem tissue by potato virus X and thereby inhibit disease development. RDR1 is inducible by salicylic acid (SA) and several other phytohormones. But although it contributes to basal resistance to tobacco mosaic virus (TMV) it is dispensable for SA-induced resistance in inoculated leaves. The laboratory accession of N. benthamiana is a natural rdr1 mutant and highly susceptible to TMV. However, TMV-induced symptoms are ameliorated in transgenic plants expressing Medicago truncatula RDR1. RESULTS: In MtRDR1-transgenic N. benthamiana plants the spread of TMV expressing the green fluorescent protein (TMV.GFP) into upper, non-inoculated, leaves was not inhibited. However, in these plants exclusion of TMV.GFP from the apical meristem and adjacent stem tissue was greater than in control plants and this exclusion effect was enhanced by SA. TMV normally kills N. benthamiana plants but although MtRDR1-transgenic plants initially displayed virus-induced necrosis they subsequently recovered. Recovery from disease was markedly enhanced by SA treatment in MtRDR1-transgenic plants whereas in control plants SA delayed but did not prevent systemic necrosis and death. Following SA treatment of MtRDR1-transgenic plants, extractable RDR enzyme activity was increased and Western blot analysis of RDR extracts revealed a band cross-reacting with an antibody raised against MtRDR1. Expression of MtRDR1 in the transgenic N. benthamiana plants was driven by a constitutive 35S promoter derived from cauliflower mosaic virus, confirmed to be non-responsive to SA. This suggests that the effects of SA on MtRDR1 are exerted at a post-transcriptional level. CONCLUSIONS: MtRDR1 inhibits severe symptom development by limiting spread of virus into the growing tips of infected plants. Thus, RDR1 may act in a similar fashion to RDR6. MtRDR1 and SA acted additively to further promote recovery from disease symptoms in MtRDR1-transgenic plants. Thus it is possible that SA promotes MtRDR1 activity and/or stability through post-transcriptional effects.

Arabidopsis basic helix-loop-helix transcription factors MYC2, MYC3, and MYC4 regulate glucosinolate biosynthesis, insect performance, and feeding behavior.

Arabidopsis thaliana plants fend off insect attack by constitutive and inducible production of toxic metabolites, such as glucosinolates (GSs). A triple mutant lacking MYC2, MYC3, and MYC4, three basic helix-loop-helix transcription factors that are known to additively control jasmonate-related defense responses, was shown to have a highly reduced expression of GS biosynthesis genes. The myc2 myc3 myc4 (myc234) triple mutant was almost completely devoid of GS and was extremely susceptible to the generalist herbivore Spodoptera littoralis. On the contrary, the specialist Pieris brassicae was unaffected by the presence of GS and preferred to feed on wild-type plants. In addition, lack of GS in myc234 drastically modified S. littoralis feeding behavior. Surprisingly, the expression of MYB factors known to regulate GS biosynthesis genes was not altered in myc234, suggesting that MYC2/MYC3/MYC4 are necessary for direct transcriptional activation of GS biosynthesis genes. To support this, chromatin immunoprecipitation analysis showed that MYC2 binds directly to the promoter of several GS biosynthesis genes in vivo. Furthermore, yeast two-hybrid and pull-down experiments indicated that MYC2/MYC3/MYC4 interact directly with GS-related MYBs. This specific MYC-MYB interaction plays a crucial role in the regulation of defense secondary metabolite production and underlines the importance of GS in shaping plant interactions with adapted and nonadapted herbivores.