Search

The transcriptional effector domains of transcription factors play a key role in controlling gene expression; however, their functional nature is poorly understood, hampering our ability to explore this fundamental dimension of gene regulatory networks. To map the trans-regulatory landscape in a complex eukaryote, we systematically characterized the putative transcriptional effector domains of over 400 Arabidopsis thaliana transcription factors for their capacity to modulate transcription. We demonstrate that transcriptional effector activity can be integrated into gene regulatory networks capable of elucidating the functional dynamics underlying gene expression patterns. We further show how our characterized domains can enhance genome engineering efforts and reveal how plant transcriptional activators share regulatory features conserved across distantly related eukaryotes. Our results provide a framework to systematically characterize the regulatory role of transcription factors at a genome-scale in order to understand the transcriptional wiring of biological systems.

Plant metabolism is more complex relative to individual microbes. In single-celled microbes, transcriptional regulation by single transcription factors (TFs) is sufficient to shift primary metabolism. Corresponding genome-level transcriptional regulatory maps of metabolism reveal the underlying design principles responsible for these shifts as a model in which master regulators largely coordinate specific metabolic pathways. Plant primary and specialized metabolism occur within innumerable cell types, and their reactions shift depending on internal and external cues. Given the importance of plants and their metabolites in providing humanity with food, fiber, and medicine, we set out to develop a genome-scale transcriptional regulatory map of Arabidopsis metabolic genes. A comprehensive set of protein-DNA interactions between Arabidopsis thaliana TFs and gene promoters in primary and specialized metabolic pathways were mapped. To demonstrate the utility of this resource, we identified and functionally validated regulators of the tricarboxylic acid (TCA) cycle. The resulting network suggests that plant metabolic design principles are distinct from those of microbes. Instead, metabolism appears to be transcriptionally coordinated via developmental- and stress-conditional processes that can coordinate across primary and specialized metabolism. These data represent the most comprehensive resource of interactions between TFs and metabolic genes in plants.

Developmental transitions are typically tightly controlled at the transcriptional level. Two of these transitions involve the induction of the embryo maturation program midway through seed development and its repression during the vegetative phase of plant growth. Very little is known about the factors responsible for this regulation during early embryogenesis, and only a couple of transcription factors have been characterized as repressors during the postgerminative phase. Arabidopsis 6b-INTERACTING PROTEIN-LIKE1 (ASIL1), a trihelix transcription factor, has been proposed to repress maturation both embryonically and postembryonically. Preliminary data also suggested that its closest paralog, ASIL2, might play a role as well. We used a transcriptomic approach, coupled with phenotypical observations, to test the hypothesis that ASIL1 and ASIL2 redundantly turn off maturation during both phases of growth. Our results indicate that, contrary to what was previously published, neither of the ASIL genes plays a role in the regulation of maturation, at any point during plant development. Analyses of gene ontology (GO)-enriched terms and published transcriptomic datasets suggest that these genes might be involved in responses during the vegetative phase to certain biotic and abiotic stresses.

Plants integrate internal and external signals to finely coordinate growth and defense for maximal fitness within a complex environment. A common model suggests that growth and defense show a trade-offs relationship driven by energy costs. However, recent studies suggest that the coordination of growth and defense likely involves more conditional and intricate connections than implied by the trade-off model. To explore how a transcription factor (TF) network may coordinate growth and defense, we used a high-throughput phenotyping approach to measure growth and flowering in a set of single and pairwise mutants previously linked to the aliphatic glucosinolate (GLS) defense pathway. Supporting a link between growth and defense, 17 of the 20 tested defense-associated TFs significantly influenced plant growth and/or flowering time. The TFs' effects were conditional upon the environment and age of the plant, and more critically varied across the growth and defense phenotypes for a given genotype. In support of the coordination model of growth and defense, the TF mutant's effects on short-chain aliphatic GLS and growth did not display a simple correlation. We propose that large TF networks integrate internal and external signals and separately modulate growth and the accumulation of the defensive aliphatic GLS.

Nitrogen is an essential macronutrient for plant growth and basic metabolic processes. The application of nitrogen-containing fertilizer increases yield, which has been a substantial factor in the green revolution1. Ecologically, however, excessive application of fertilizer has disastrous effects such as eutrophication2. A better understanding of how plants regulate nitrogen metabolism is critical to increase plant yield and reduce fertilizer overuse. Here we present a transcriptional regulatory network and twenty-one transcription factors that regulate the architecture of root and shoot systems in response to changes in nitrogen availability. Genetic perturbation of a subset of these transcription factors revealed coordinate transcriptional regulation of enzymes involved in nitrogen metabolism. Transcriptional regulators in the network are transcriptionally modified by feedback via genetic perturbation of nitrogen metabolism. The network, genes and gene-regulatory modules identified here will prove critical to increasing agricultural productivity.

Plants have evolved tightly regulated signaling networks to respond and adapt to environmental perturbations, but the nature of the signaling hub(s) involved have remained an enigma. We have previously established that methylerythritol cyclodiphosphate (MEcPP), a precursor of plastidial isoprenoids and a stress-specific retrograde signaling metabolite, enables cellular readjustments for high-order adaptive functions. Here, we specifically show that MEcPP promotes two Brassicaceae-specific traits, namely endoplasmic reticulum (ER) body formation and induction of indole glucosinolate (IGs) metabolism selectively, via transcriptional regulation of key regulators NAI1 for ER body formation and MYB51/122 for IGs biosynthesis). The specificity of MEcPP is further confirmed by the lack of induction of wound-inducible ER body genes as well as IGs by other altered methylerythritol phosphate pathway enzymes. Genetic analyses revealed MEcPP-mediated COI1-dependent induction of these traits. Moreover, MEcPP signaling integrates the biosynthesis and hydrolysis of IGs through induction of nitrile-specifier protein1 and reduction of the suppressor, ESM1, and production of simple nitriles as the bioactive end product. The findings position the plastidial metabolite, MEcPP, as the initiation hub, transducing signals to adjust the activity of hard-wired gene circuitry to expand phytochemical diversity and alter the associated subcellular structure required for functionality of the secondary metabolites, thereby tailoring plant stress responses.

Despite the growing number of studies showing that genotype × environment and epistatic interactions control fitness, the influences of epistasis × environment interactions on adaptive trait evolution remain largely uncharacterized. Across three field trials, we quantified aliphatic glucosinolate (GSL) defense chemistry, leaf damage, and relative fitness using mutant lines of Arabidopsis thaliana varying at pairs of causal aliphatic GSL defense genes to test the impact of epistatic and epistasis × environment interactions on adaptive trait variation. We found that aliphatic GSL accumulation was primarily influenced by additive and epistatic genetic variation, leaf damage was primarily influenced by environmental variation and relative fitness was primarily influenced by epistasis and epistasis × environment interactions. Epistasis × environment interactions accounted for up to 48% of the relative fitness variation in the field. At a single field site, the impact of epistasis on relative fitness varied significantly over 2 yr, showing that epistasis × environment interactions within a location can be temporally dynamic. These results suggest that the environmental dependency of epistasis can profoundly influence the response to selection, shaping the adaptive trajectories of natural populations in complex ways, and deserves further consideration in future evolutionary studies.

BREVIPEDICELLUS (BP) encodes a class-I KNOTTED1-like homeobox (KNOX) transcription factor that plays a critical role in conditioning a replication competent state in the apical meristem, and it also governs growth and cellular differentiation in internodes and pedicels. To search for factors that modify BP signaling, we conducted a suppressor screen on bp er (erecta) plants and identified a mutant that ameliorates many of the pleiotropic defects of the parent line. Map based cloning and complementation studies revealed that the defect lies in the FILAMENTOUS FLOWER (FIL) gene, a member of the YABBY family of transcriptional regulators that contribute to meristem organization and function, phyllotaxy, leaf and floral organ growth and polarity, and are also known to repress KNOX gene expression. Genetic and cytological analyses of the fil-10 suppressor line indicate that the role of FIL in promoting growth is independent of its previously characterized influences on meristem identity and lateral organ polarity, and likely occurs non-cell-autonomously from superior floral organs. Transcription profiling of inflorescences revealed that FIL downregulates numerous transcription factors which in turn may subordinately regulate inflorescence architecture. In addition, FIL, directly or indirectly, activates over a dozen genes involved in glucosinolate production in part by activating MYB28, a known activator of many aliphatic glucosinolate biosynthesis genes. In the bp er fil-10 suppressor mutant background, enhanced expression of CYP71A13, AMIDASE1 (AMI) and NITRILASE genes suggest that auxin levels can be modulated by shunting glucosinolate metabolites into the IAA biosynthetic pathway, and increased IAA levels in the bp er fil-10 suppressor accompany enhanced internode and pedicel elongation. We propose that FIL acts to oppose KNOX1 gene function through a complex regulatory network that involves changes in secondary metabolites and auxin.

Studying the genetic basis of variation in plant metabolism has been greatly facilitated by genomic and metabolic profiling advances. In this study, we use metabolomics and growth measurements to map QTL in rice, a major staple crop. Previous rice metabolism studies have largely focused on identifying genes controlling major effect loci. To complement these studies, we conducted a replicated metabolomics analysis on a japonica (Lemont) by indica (Teqing) rice recombinant inbred line population and focused on the genetic variation for primary metabolism. Using independent replicated studies, we show that in contrast to other rice studies, the heritability of primary metabolism is similar to Arabidopsis. The vast majority of metabolic QTLs had small to moderate effects with significant polygenic epistasis. Two metabolomics QTL hotspots had opposing effects on carbon and nitrogen rich metabolites suggesting that they may influence carbon and nitrogen partitioning, with one locus co-localizing with SUSIBA2 (WRKY78). Comparing QTLs for metabolomic and a variety of growth related traits identified few overlaps. Interestingly, the rice population displayed fewer loci controlling stochastic variation for metabolism than was found in Arabidopsis. Thus, it is possible that domestication has differentially impacted stochastic metabolite variation more than average metabolite variation.

Plant cell walls are important barriers against microbial pathogens. Cell walls of Arabidopsis thaliana leaves contain three major types of polysaccharides: cellulose, various hemicelluloses, and pectins. UDP-D-galacturonic acid, the key building block of pectins, is produced from the precursor UDP-D-glucuronic acid by the action of glucuronate 4-epimerases (GAEs). Pseudomonas syringae pv maculicola ES4326 (Pma ES4326) repressed expression of GAE1 and GAE6 in Arabidopsis, and immunity to Pma ES4326 was compromised in gae6 and gae1 gae6 mutant plants. These plants had brittle leaves and cell walls of leaves had less galacturonic acid. Resistance to specific Botrytis cinerea isolates was also compromised in gae1 gae6 double mutant plants. Although oligogalacturonide (OG)-induced immune signaling was unaltered in gae1 gae6 mutant plants, immune signaling induced by a commercial pectinase, macerozyme, was reduced. Macerozyme treatment or infection with B. cinerea released less soluble uronic acid, likely reflecting fewer OGs, from gae1 gae6 cell walls than from wild-type Col-0. Although both OGs and macerozyme-induced immunity to B. cinerea in Col-0, only OGs also induced immunity in gae1 gae6. Pectin is thus an important contributor to plant immunity, and this is due at least in part to the induction of immune responses by soluble pectin, likely OGs, that are released during plant-pathogen interactions.

A key limitation in modern biology is the ability to rapidly identify genes underlying newly identified complex phenotypes. Genome wide association studies (GWAS) have become an increasingly important approach for dissecting natural variation by associating phenotypes with genotypes at a genome wide level. Recent work is showing that the Arabidopsis thaliana defense metabolite, allyl glucosinolate (GSL), may provide direct feedback regulation, linking defense metabolism outputs to the growth, and defense responses of the plant. However, there is still a need to identify genes that underlie this process. To start developing a deeper understanding of the mechanism(s) that modulate the ability of exogenous allyl GSL to alter growth and defense, we measured changes in plant biomass and defense metabolites in a collection of natural 96 A. thaliana accessions fed with 50 μM of allyl GSL. Exogenous allyl GSL was introduced exclusively to the roots and the compound transported to the leaf leading to a wide range of heritable effects upon plant biomass and endogenous GSL accumulation. Using natural variation we conducted GWAS to identify a number of new genes which potentially control allyl responses in various plant processes. This is one of the first instances in which this approach has been successfully utilized to begin dissecting a novel phenotype to the underlying molecular/polygenic basis.

Natural populations persist in complex environments, where biotic stressors, such as pathogen and insect communities, fluctuate temporally and spatially. These shifting biotic pressures generate heterogeneous selective forces that can maintain standing natural variation within a species. To directly test if genes containing causal variation for the Arabidopsis thaliana defensive compounds, glucosinolates (GSL) control field fitness and are therefore subject to natural selection, we conducted a multi-year field trial using lines that vary in only specific causal genes. Interestingly, we found that variation in these naturally polymorphic GSL genes affected fitness in each of our environments but the pattern fluctuated such that highly fit genotypes in one trial displayed lower fitness in another and that no GSL genotype or genotypes consistently out-performed the others. This was true both across locations and within the same location across years. These results indicate that environmental heterogeneity may contribute to the maintenance of GSL variation observed within Arabidopsis thaliana.

The evolution of plant metabolic diversity is largely driven by gene duplication and ensuing sub-functionalization and/or neo-functionalization to generate new enzymatic activities. However, it is not clear whether the transcription factors (TFs) regulating these new enzyme encoding genes were required to co-evolve with these genes in a similar fashion or if these new genes can be captured by existing conserved TFs to provide the appropriate expression pattern. In this study, we found two conserved TFs, MYB115, and MYB118, co-expressed with the key enzyme encoding genes in the newly evolved benzoyloxy glucosinolate (GLS) pathway. These TFs interacted with the promoters of the GLS biosynthetic genes and negatively influenced their expression. Similarly, the GLS profiles of these two TFs knockouts showed that they influenced the aliphatic GLS accumulation within seed, leaf and flower, while they mainly expressed in seeds. Further studies indicated that they are functionally redundant and epistatically interact to control the transcription of GLS genes. Complementation study confirmed their roles in regulating the aliphatic GLS biosynthesis. These results suggest that the newly evolved enzyme encoding genes for novel metabolites can be regulated by conserved TFs, which helps to improve our model for newly evolved genes regulation.

A key unanswered question in plant biology is how a plant regulates metabolism to maximize performance across an array of biotic and abiotic environmental stresses. In this study, we addressed the potential breadth of transcriptional regulation that can alter accumulation of the defensive glucosinolate metabolites in Arabidopsis (Arabidopsis thaliana). A systematic yeast one-hybrid study was used to identify hundreds of unique potential regulatory interactions with a nearly complete complement of 21 promoters for the aliphatic glucosinolate pathway. Conducting high-throughput phenotypic validation, we showed that >75% of tested transcription factor (TF) mutants significantly altered the accumulation of the defensive glucosinolates. These glucosinolate phenotypes were conditional upon the environment and tissue type, suggesting that these TFs may allow the plant to tune its defenses to the local environment. Furthermore, the pattern of TF/promoter interactions could partially explain mutant phenotypes. This work shows that defense chemistry within Arabidopsis has a highly intricate transcriptional regulatory system that may allow for the optimization of defense metabolite accumulation across a broad array of environments.

A central goal of systems biology is to develop models that are both predictive and accurately describe the biological system. One complexity to this endeavor is that it is possible to develop models that appear predictive even if they use far fewer components than the biological system itself uses for the same process. This problem also occurs in quantitative genetics where it is often possible to describe the variation in a system using fewer genes than are actually variable due to the complications of linkage between causal polymorphisms and population structure. Thus, there is a crucial need to begin an empirical investigation into the true number of components that are used by biological systems to determine a phenotypic outcome. In this study, we use a meta-analysis of directly comparable metabolomics quantitative studies using quantitative trait locus mapping and genome wide association mapping to show that it is currently not possible to estimate how many genetic loci are truly polymorphic within Arabidopsis thaliana. Our analysis shows that it would require the analysis of at least a 1000 line bi-parental population to begin to estimate how many polymorphic loci control metabolic variation within Arabidopsis. Understanding the base number of loci that are actually involved in determining variation in metabolic systems is fundamental to developing systems models that are truly reflective of how metabolism is modulated within a living organism.

Regulation of primary metabolism is a central mechanism by which plants coordinate their various responses to biotic and abiotic challenge. To identify genes responsible for natural variation in primary metabolism, we focused on cloning a locus from Arabidopsis thaliana that influences the level of TCA cycle metabolites in planta. We found that the Met.V.67 locus was controlled by natural variation in METABOLIC NETWORK MODULATOR 1 (MNM1), which encoded an AT-hook motif-containing protein that was unique to the Brassicales lineage. MNM1 had wide ranging effects on plant metabolism and displayed a tissue expression pattern that was suggestive of a function in sink tissues. Natural variation within MNM1 had differential effects during a diurnal time course, and this temporal dependency was supported by analysis of T-DNA insertion and over-expression lines for MNM1. Thus, the cloning of a natural variation locus specifically associated with primary metabolism allowed us to identify MNM1 as a lineage-specific modulator of primary metabolism, suggesting that the regulation of primary metabolism can change during evolution.

Understanding genome to phenotype linkages has been greatly enabled by genomic sequencing. However, most genome analysis is typically confined to the nuclear genome. We conducted a metabolomic QTL analysis on a reciprocal RIL population structured to examine how variation in the organelle genomes affects phenotypic variation. This showed that the cytoplasmic variation had effects similar to, if not larger than, the largest individual nuclear locus. Inclusion of cytoplasmic variation into the genetic model greatly increased the explained phenotypic variation. Cytoplasmic genetic variation was a central hub in the epistatic network controlling the plant metabolome. This epistatic influence manifested such that the cytoplasmic background could alter or hide pairwise epistasis between nuclear loci. Thus, cytoplasmic genetic variation plays a central role in controlling natural variation in metabolomic networks. This suggests that cytoplasmic genomes must be included in any future analysis of natural variation. DOI: http://dx.doi.org/10.7554/eLife.00776.001.

Aberrant expression of classical cadherins has been observed in tumor invasion and metastasis, but its involvement in cervical carcinogenesis and cancer progression is not clear. We investigated E-, P- and N-cadherin expression and its significance in cervical squamous cell carcinoma (SCC) and cervical intraepithelial neoplasia (CIN). This retrospective study enrolled 508 patients admitted to Women's Hospital, School of Medicine, Zhejiang University with cervical lesions between January 2006 and December 2010. Immunochemical staining was performed in 98 samples of normal cervical epithelium (NC), 283 of CIN, and 127 of early-stage SCC. The association of cadherin staining with clinical characteristics and survival of the patients was evaluated by univariate and multivariate analysis. We found gradients of decreasing E-cadherin expression and increasing P-cadherin expression from NC through CIN to SCC. Aberrant E-cadherin and P-cadherin expression were significantly associated with clinical parameters indicating poor prognosis and shorter patient survival. Interestingly, we found very low levels of positive N-cadherin expression in CIN and SCC tissues that were not related to CIN or cancer. Pearson chi-square tests showed that E-cadherin expression in SCC was inversely correlated with P-cadherin expression (E-P switch), and was not correlated with N-cadherin expression. More important, patients with tissues exhibiting an E-P switch in expression had highly aggressive phenotypes and poorer prognosis than those without E-P switch expression. Our findings suggest that E-cadherin and P-cadherin, but not N-cadherin staining, might be useful in diagnosing CIN and for predicting prognosis in patients with early-stage SCC. [ABSTRACT FROM AUTHOR]