Enough is enough: feedback control of specialized metabolism

Recent advances in our understanding of plant metabolism have highlighted the significance of specialized metabolites in the regulation of gene expression associated with biosynthetic networks. This opinion article focuses on the molecular mechanisms of small-molecule-mediated feedback regulation at the transcriptional level and its potential modes of action, including metabolite signal perception, the nature of the sensor, and the signaling transduction mechanisms leading to transcriptional and post-transcriptional regulation, based on evidence available from plants and other kingdoms of life. We also discuss the challenges associated with identifying the occurrences, effects, and localization of small molecule-protein interactions. Further understanding of small-molecule-controlled metabolic fluxes will enable rational design of transcriptional regulation systems in metabolic engineering to produce high-value specialized metabolites.

Should I stay or should I go? Trafficking of plant extra-nuclear transcription factors

At the heart of all biological processes lies the control of nuclear gene expression, which is primarily achieved through the action of transcription factors (TFs) that generally contain a nuclear localization signal (NLS) to facilitate their transport into the nucleus. However, some TFs reside in the cytoplasm in a transcriptionally inactive state and only enter the nucleus in response to specific signals, which in plants include biotic or abiotic stresses. These extra-nuclear TFs can be found in the cytosol or associated with various membrane systems, including the endoplasmic reticulum and plasma membrane. They may be integral proteins with transmembrane domains or associate peripherally with the lipid bilayer via acylation or membrane-binding domains. Although over 30 plant TFs, most of them involved in stress responses, have been experimentally shown to reside outside the nucleus, computational predictions suggest that this number is much larger. Understanding how extra-nuclear TFs are trafficked into the nucleus is essential for reconstructing transcriptional regulatory networks that govern major cellular pathways in response to biotic and abiotic signals. Here, we provide a perspective on what is known on plant extranuclear-nuclear TF retention, nuclear trafficking, and the post-translational modifications that ultimately enable them to regulate gene expression upon entering the nucleus.

Multiscale physiological responses to nitrogen supplementation of maize hybrids

Maize (Zea mays) production systems are heavily reliant on the provision of managed inputs such as fertilizers to maximize growth and yield. Hence, the effective use of nitrogen (N) fertilizer is crucial to minimize the associated financial and environmental costs, as well as maximize yield. However, how to effectively utilize N inputs for increased grain yields remains a substantial challenge for maize growers that requires a deeper understanding of the underlying physiological responses to N fertilizer application. We report a multiscale investigation of five field-grown maize hybrids under low or high N supplementation regimes that includes the quantification of phenolic and prenyl-lipid compounds, cellular ultrastructural features, and gene expression traits at three developmental stages of growth. Our results reveal that maize perceives the lack of supplemented N as a stress and, when provided with additional N, will prolong vegetative growth. However, the manifestation of the stress and responses to N supplementation are highly hybrid-specific. Eight genes were differentially expressed in leaves in response to N supplementation in all tested hybrids and at all developmental stages. These genes represent potential biomarkers of N status and include two isoforms of Thiamine Thiazole Synthase involved in vitamin B1 biosynthesis. Our results uncover a detailed view of the physiological responses of maize hybrids to N supplementation in field conditions that provides insight into the interactions between management practices and the genetic diversity within maize.

Prioritizing Metabolic Gene Regulators through Multi-Omic Network Integration in Maize

Elucidating gene regulatory networks (GRNs) is a major area of study within plant systems biology. Phenotypic traits are intricately linked to specific gene expression profiles. These expression patterns arise primarily from regulatory connections between sets of transcription factors (TFs) and their target genes. In this study, we integrated publicly available co-expression networks derived from more than 6,000 RNA-seq samples, 283 protein-DNA interaction assays, and 16 million of SNPs used to identify expression quantitative loci (eQTL), to construct TF-target networks. In total, we analyzed ~4.6M interactions to generate four distinct types of TF-target networks: co-expression, protein-DNA interaction (PDI), trans-expression quantitative loci (trans-eQTL), and cis-eQTL combined with PDIs. To improve the functional annotation of TFs based on its target genes, we implemented three different strategies to integrate these four types of networks. We subsequently evaluated the effectiveness of our method through loss-of function mutant and random networks. The multi-network integration allowed us to identify transcriptional regulators of hormone-, metabolic- and development-related processes. Finally, using the topological properties of the fully integrated network, we identified potentially functional redundant TF paralogs. Our findings retrieved functions previously documented for numerous TFs and revealed novel functions that are crucial for informing the design of future experiments. The approach here-described lays the foundation for the integration of multi-omic datasets in maize and other plant systems.

Catalytically impaired chalcone isomerase retains flavonoid biosynthetic capacity

Catalytically impaired chalcone isomerases having no appreciable activity in vitro retain the ability to robustly promote flavonoid biosynthesis in planta.

Transcriptome reprogramming through alternative splicing triggered by apigenin drives cell death in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is characterized by its aggressiveness and resistance to cancer-specific transcriptome alterations. Alternative splicing (AS) is a major contributor to the diversification of cancer-specific transcriptomes. The TNBC transcriptome landscape is characterized by aberrantly spliced isoforms that promote tumor growth and resistance, underscoring the need to identify approaches that reprogram AS circuitry towards transcriptomes, favoring a delay in tumorigenesis or responsiveness to therapy. We have previously shown that flavonoid apigenin is associated with splicing factors, including heterogeneous nuclear ribonucleoprotein A2 (hnRNPA2). Here, we showed that apigenin reprograms TNBC-associated AS transcriptome-wide. The AS events affected by apigenin were statistically enriched in hnRNPA2 substrates. Comparative transcriptomic analyses of human TNBC tumors and non-tumor tissues showed that apigenin can switch cancer-associated alternative spliced isoforms (ASI) to those found in non-tumor tissues. Apigenin preferentially affects the splicing of anti-apoptotic and proliferation factors, which are uniquely observed in cancer cells, but not in non-tumor cells. Apigenin switches cancer-associated aberrant ASI in vivo in TNBC xenograft mice by diminishing proliferation and increasing pro-apoptotic ASI. In accordance with these findings, apigenin increased apoptosis and reduced tumor proliferation, thereby halting TNBC growth in vivo. Our results revealed that apigenin reprograms transcriptome-wide TNBC-specific AS, thereby inducing apoptosis and hindering tumor growth. These findings underscore the impactful effects of nutraceuticals in altering cancer transcriptomes, offering new options to influence outcomes in TNBC treatments.

Mutator transposon insertions within maize genes often provide a novel outward reading promoter

The highly active family of Mutator (Mu) DNA transposons has been widely used for forward and reverse genetics in maize. There are examples of Mu-suppressible alleles that result in conditional phenotypic effects based on the activity of Mu. Phenotypes from these Mu-suppressible mutations are observed in Mu-active genetic backgrounds, but absent when Mu activity is lost. For some Mu-suppressible alleles, phenotypic suppression likely results from an outward-reading promoter within Mu that is only active when the autonomous Mu element is silenced or lost. We isolated 35 Mu alleles from the UniformMu population that represent insertions in 24 different genes. Most of these mutant alleles are due to insertions within gene coding sequences, but several 5' UTR and intron insertions were included. RNA-seq and de novo transcript assembly were utilized to document the transcripts produced from 33 of these Mu insertion alleles. For 20 of the 33 alleles, there was evidence of transcripts initiating within the Mu sequence reading through the gene. This outward-reading promoter activity was detected in multiple types of Mu elements and does not depend on the orientation of Mu. Expression analyses of Mu-initiated transcripts revealed the Mu promoter often provides gene expression levels and patterns that are similar to the wild-type gene. These results suggest the Mu promoter may represent a minimal promoter that can respond to gene cis-regulatory elements. Findings from this study have implications for maize researchers using the UniformMu population, and more broadly highlight a strategy for transposons to co-exist with their host.

Metabolic and transcriptomic study of pennycress natural variation identifies targets for oil improvement

Pennycress (Thlaspi arvense L.), a member of the Brassicaceae family, produces seed oil high in erucic acid, suitable for biodiesel and aviation fuel. Although pennycress, a winter annual, could be grown as a dedicated bioenergy crop, an increase in its seed oil content is required to improve its economic competitiveness. The success of crop improvement relies upon finding the right combination of biomarkers and targets, and the best genetic engineering and/or breeding strategies. In this work, we combined biomass composition with metabolomic and transcriptomic studies of developing embryos from 22 pennycress natural variants to identify targets for oil improvement. The selected accession collection presented diverse levels of fatty acids at maturity ranging from 29% to 41%. Pearson correlation analyses, weighted gene co-expression network analysis and biomarker identifications were used as complementary approaches to detect associations between metabolite level or gene expression and oil content at maturity. The results indicated that improving seed oil content can lead to a concomitant increase in the proportion of erucic acid without affecting the weight of embryos. Processes, such as carbon partitioning towards the chloroplast, lipid metabolism, photosynthesis, and a tight control of nitrogen availability, were found to be key for oil improvement in pennycress. Besides identifying specific targets, our results also provide guidance regarding the best timing for their modification, early or middle maturation. Thus, this work lays out promising strategies, specific for pennycress, to accelerate the successful development of lines with increased seed oil content for biofuel applications.

Evolution and diversification of the ACT-like domain associated with plant basic helix-loop-helix transcription factors

Basic helix-loop-helix (bHLH) proteins are one of the largest families of transcription factor (TF) in eukaryotes, and ~30% of all flowering plants' bHLH TFs contain the aspartate kinase, chorismate mutase, and TyrA (ACT)-like domain at variable distances C-terminal from the bHLH. However, the evolutionary history and functional consequences of the bHLH/ACT-like domain association remain unknown. Here, we show that this domain association is unique to the plantae kingdom with green algae (chlorophytes) harboring a small number of bHLH genes with variable frequency of ACT-like domain's presence. bHLH-associated ACT-like domains form a monophyletic group, indicating a common origin. Indeed, phylogenetic analysis results suggest that the association of ACT-like and bHLH domains occurred early in Plantae by recruitment of an ACT-like domain in a common ancestor with widely distributed ACT DOMAIN REPEAT (ACR) genes by an ancestral bHLH gene. We determined the functional significance of this association by showing that Chlamydomonas reinhardtii ACT-like domains mediate homodimer formation and negatively affect DNA binding of the associated bHLH domains. We show that, while ACT-like domains have experienced faster selection than the associated bHLH domain, their rates of evolution are strongly and positively correlated, suggesting that the evolution of the ACT-like domains was constrained by the bHLH domains. This study proposes an evolutionary trajectory for the association of ACT-like and bHLH domains with the experimental characterization of the functional consequence in the regulation of plant-specific processes, highlighting the impacts of functional domain coevolution.

Methodology for Constructing a Knowledgebase for Plant Gene Regulation Information

The amount of biological data is growing at a rapid pace as many high-throughput omics technologies and data pipelines are developed. This is resulting in the growth of databases for DNA and protein sequences, gene expression, protein accumulation, structural, and localization information. The diversity and multi-omics nature of such bioinformatic data requires well-designed databases for flexible organization and presentation. Besides general-purpose online bioinformatic databases, users need narrowly focused online databases to quickly access a meaningful collection of related data for their research. Here, we describe the methodology used to implement a plant gene regulatory knowledgebase, with data, query, and tool features, as well as the ability to expand to accommodate future datasets. We exemplify this methodology for the GRASSIUS knowledgebase, but it is applicable to developing and updating similar plant gene regulatory knowledgebases. GRASSIUS organizes and presents gene regulatory data from grass species with a central focus on maize (Zea mays). The main class of data presented include not only the families of transcription factors (TFs) and co-regulators (CRs) but also protein-DNA interaction data, where available.

Effective Mechanisms for Improving Seed Oil Production in Pennycress (Thlaspi arvense L.) Highlighted by Integration of Comparative Metabolomics and Transcriptomics

Pennycress is a potentially lucrative biofuel crop due to its high content of long-chain unsaturated fatty acids, and because it uses non-conventional pathways to achieve efficient oil production. However, metabolic engineering is required to improve pennycress oilseed content and make it an economically viable source of aviation fuel. Research is warranted to determine if further upregulation of these non-conventional pathways could improve oil production within the species even more, which would indicate these processes serve as promising metabolic engineering targets and could provide the improvement necessary for economic feasibility of this crop. To test this hypothesis, we performed a comparative biomass, metabolomic, and transcriptomic analyses between a high oil accession (HO) and low oil accession (LO) of pennycress to assess potential factors required to optimize oil content. An evident reduction in glycolysis intermediates, improved oxidative pentose phosphate pathway activity, malate accumulation in the tricarboxylic acid cycle, and an anaplerotic pathway upregulation were noted in the HO genotype. Additionally, higher levels of threonine aldolase transcripts imply a pyruvate bypass mechanism for acetyl-CoA production. Nucleotide sugar and ascorbate accumulation also were evident in HO, suggesting differential fate of associated carbon between the two genotypes. An altered transcriptome related to lipid droplet (LD) biosynthesis and stability suggests a contribution to a more tightly-packed LD arrangement in HO cotyledons. In addition to the importance of central carbon metabolism augmentation, alternative routes of carbon entry into fatty acid synthesis and modification, as well as transcriptionally modified changes in LD regulation, are key aspects of metabolism and storage associated with economically favorable phenotypes of the species.

Normalizing and Correcting Variable and Complex LC-MS Metabolomic Data with the R Package pseudoDrift

In biological research domains, liquid chromatography-mass spectroscopy (LC-MS) has prevailed as the preferred technique for generating high quality metabolomic data. However, even with advanced instrumentation and established data acquisition protocols, technical errors are still routinely encountered and can pose a significant challenge to unveiling biologically relevant information. In large-scale studies, signal drift and batch effects are how technical errors are most commonly manifested. We developed pseudoDrift, an R package with capabilities for data simulation and outlier detection, and a new training and testing approach that is implemented to capture and to optionally correct for technical errors in LC-MS metabolomic data. Using data simulation, we demonstrate here that our approach performs equally as well as existing methods and offers increased flexibility to the researcher. As part of our study, we generated a targeted LC-MS dataset that profiled 33 phenolic compounds from seedling stem tissue in 602 genetically diverse non-transgenic maize inbred lines. This dataset provides a unique opportunity to investigate the dynamics of specialized metabolism in plants.

Natural variation and improved genome annotation of the emerging biofuel crop field pennycress (Thlaspi arvense)

The Brassicaceae family comprises more than 3,700 species with a diversity of phenotypic characteristics, including seed oil content and composition. Recently, the global interest in Thlaspi arvense L. (pennycress) has grown as the seed oil composition makes it a suitable source for biodiesel and aviation fuel production. However, many wild traits of this species need to be domesticated to make pennycress ideal for cultivation. Molecular breeding and engineering efforts require the availability of an accurate genome sequence of the species. Here, we describe pennycress genome annotation improvements, using a combination of long- and short-read transcriptome data obtained from RNA derived from embryos of 22 accessions, in addition to public genome and gene expression information. Our analysis identified 27,213 protein-coding genes, as well as on average 6,188 biallelic SNPs. In addition, we used the identified SNPs to evaluate the population structure of our accessions. The data from this analysis support that the accession Ames 32872, originally from Armenia, is highly divergent from the other accessions, while the accessions originating from Canada and the United States cluster together. When we evaluated the likely signatures of natural selection from alternative SNPs, we found 7 candidate genes under likely recent positive selection. These genes are enriched with functions related to amino acid metabolism and lipid biosynthesis and highlight possible future targets for crop improvement efforts in pennycress.

Exploring Camelina sativa lipid metabolism regulation by combining gene co-expression and DNA affinity purification analyses

Camelina (Camelina sativa) is an annual oilseed plant that is gaining momentum as a biofuel cover crop. Understanding gene regulatory networks is essential to deciphering plant metabolic pathways, including lipid metabolism. Here, we take advantage of a growing collection of gene expression datasets to predict transcription factors (TFs) associated with the control of Camelina lipid metabolism. We identified approximately 350 TFs highly co-expressed with lipid-related genes (LRGs). These TFs are highly represented in the MYB, AP2/ERF, bZIP, and bHLH families, including a significant number of homologs of well-known Arabidopsis lipid and seed developmental regulators. After prioritizing the top 22 TFs for further validation, we identified DNA-binding sites and predicted target genes for 16 out of the 22 TFs tested using DNA affinity purification followed by sequencing (DAP-seq). Enrichment analyses of targets supported the co-expression prediction for most TF candidates, and the comparison to Arabidopsis revealed some common themes, but also aspects unique to Camelina. Within the top potential lipid regulators, we identified CsaMYB1, CsaABI3AVP1-2, CsaHB1, CsaNAC2, CsaMYB3, and CsaNAC1 as likely involved in the control of seed fatty acid elongation and CsaABI3AVP1-2 and CsabZIP1 as potential regulators of the synthesis and degradation of triacylglycerols (TAGs), respectively. Altogether, the integration of co-expression data and DNA-binding assays permitted us to generate a high-confidence and short list of Camelina TFs involved in the control of lipid metabolism during seed development.

Microautophagy Mediates Vacuolar Delivery of Storage Proteins in Maize Aleurone Cells

The molecular machinery orchestrating microautophagy, whereby eukaryotic cells sequester autophagic cargo by direct invagination of the vacuolar/lysosomal membrane, is still largely unknown, especially in plants. Here, we demonstrate microautophagy of storage proteins in the maize aleurone cells of the endosperm and analyzed proteins with potential regulatory roles in this process. Within the cereal endosperm, starchy endosperm cells accumulate storage proteins (mostly prolamins) and starch whereas the peripheral aleurone cells store oils, storage proteins, and specialized metabolites. Although both cell types synthesize prolamins, they employ different pathways for their subcellular trafficking. Starchy endosperm cells accumulate prolamins in protein bodies within the endoplasmic reticulum (ER), whereas aleurone cells deliver prolamins to vacuoles via an autophagic mechanism, which we show is by direct association of ER prolamin bodies with the tonoplast followed by engulfment via microautophagy. To identify candidate proteins regulating this process, we performed RNA-seq transcriptomic comparisons of aleurone and starchy endosperm tissues during seed development and proteomic analysis on tonoplast-enriched fractions of aleurone cells. From these datasets, we identified 10 candidate proteins with potential roles in membrane modification and/or microautophagy, including phospholipase-Dα5 and a possible EUL-like lectin. We found that both proteins increased the frequency of tonoplast invaginations when overexpressed in Arabidopsis leaf protoplasts and are highly enriched at the tonoplast surface surrounding ER protein bodies in maize aleurone cells, thus supporting their potential connections to microautophagy. Collectively, this candidate list now provides useful tools to study microautophagy in plants.

Modeling temporal and hormonal regulation of plant transcriptional response to wounding

Plants respond to wounding stress by changing gene expression patterns and inducing the production of hormones including jasmonic acid. This wounding transcriptional response activates specialized metabolism pathways such as the glucosinolate pathways in Arabidopsis thaliana. While the regulatory factors and sequences controlling a subset of wound-response genes are known, it remains unclear how wound response is regulated globally. Here, we how these responses are regulated by incorporating putative cis-regulatory elements, known transcription factor binding sites, in vitro DNA affinity purification sequencing, and DNase I hypersensitive sites to predict genes with different wound-response patterns using machine learning. We observed that regulatory sites and regions of open chromatin differed between genes upregulated at early and late wounding time-points as well as between genes induced by jasmonic acid and those not induced. Expanding on what we currently know, we identified cis-elements that improved model predictions of expression clusters over known binding sites. Using a combination of genome editing, in vitro DNA-binding assays, and transient expression assays using native and mutated cis-regulatory elements, we experimentally validated four of the predicted elements, three of which were not previously known to function in wound-response regulation. Our study provides a global model predictive of wound response and identifies new regulatory sequences important for wounding without requiring prior knowledge of the transcriptional regulators.

Cis-regulatory sequences in plants: Their importance, discovery, and future challenges

The identification and characterization of cis-regulatory DNA sequences and how they function to coordinate responses to developmental and environmental cues is of paramount importance to plant biology. Key to these regulatory processes are cis-regulatory modules (CRMs), which include enhancers and silencers. Despite the extraordinary advances in high-quality sequence assemblies and genome annotations, the identification and understanding of CRMs, and how they regulate gene expression, lag significantly behind. This is especially true for their distinguishing characteristics and activity states. Here, we review the current knowledge on CRMs and breakthrough technologies enabling identification, characterization, and validation of CRMs; we compare the genomic distributions of CRMs with respect to their target genes between different plant species, and discuss the role of transposable elements harboring CRMs in the evolution of gene expression. This is an exciting time to study cis-regulomes in plants; however, significant existing challenges need to be overcome to fully understand and appreciate the role of CRMs in plant biology and in crop improvement.

Prediction of conserved and variable heat and cold stress response in maize using cis-regulatory information

Changes in gene expression are important for responses to abiotic stress. Transcriptome profiling of heat- or cold-stressed maize genotypes identifies many changes in transcript abundance. We used comparisons of expression responses in multiple genotypes to identify alleles with variable responses to heat or cold stress and to distinguish examples of cis- or trans-regulatory variation for stress-responsive expression changes. We used motifs enriched near the transcription start sites (TSSs) for thermal stress-responsive genes to develop predictive models of gene expression responses. Prediction accuracies can be improved by focusing only on motifs within unmethylated regions near the TSS and vary for genes with different dynamic responses to stress. Models trained on expression responses in a single genotype and promoter sequences provided lower performance when applied to other genotypes but this could be improved by using models trained on data from all three genotypes tested. The analysis of genes with cis-regulatory variation provides evidence for structural variants that result in presence/absence of transcription factor binding sites in creating variable responses. This study provides insights into cis-regulatory motifs for heat- and cold-responsive gene expression and defines a framework for developing models to predict expression responses across multiple genotypes.

A hydrophobic residue stabilizes dimers of regulatory ACT-like domains in plant basic helix-loop-helix transcription factors

About a third of the plant basic helix–loop–helix (bHLH) transcription factors harbor a C-terminal aspartate kinase, chorismate mutase, and TyrA (ACT)-like domain, which was originally identified in the maize R regulator of anthocyanin biosynthesis, where it modulates the ability of the bHLH to dimerize and bind DNA. Characterization of other bHLH ACT-like domains, such as the one in the Arabidopsis R ortholog, GL3, has not definitively confirmed dimerization, raising the question of the overall role of this potential regulatory domain. To learn more, we compared the dimerization of the ACT-like domains of R (R^ACT) and GL3 (GL3^ACT). We show that R^ACT dimerizes with a dissociation constant around 100 nM, over an order of magnitude stronger than GL3^ACT. Structural predictions combined with mutational analyses demonstrated that V568, located in a hydrophobic pocket in R^ACT, is important: when mutated to the Ser residue present in GL3^ACT, dimerization affinity dropped by almost an order of magnitude. The converse S595V mutation in GL3^ACT significantly increased the dimerization strength. We cloned and assayed dimerization for all identified maize ACT-like domains and determined that 12 of 42 formed heterodimers in yeast two-hybrid assays, irrespective of whether they harbored V568, which was often replaced by other aliphatic amino acids. Moreover, we determined that the presence of polar residues at that position occurs only in a small subset of anthocyanin regulators. The combined results provide new insights into possibly regulatory mechanisms and suggest that many of the other plant ACT-like domains associate to modulate fundamental cellular processes.

Rhamnose in plants - from biosynthesis to diverse functions

In plants, the deoxy sugar l-rhamnose is widely present as rhamnose-containing polymers in cell walls and as part of the decoration of various specialized metabolites. Here, we review the current knowledge on the distribution of rhamnose, highlighting the differences between what is known in dicotyledoneuos compared to commelinid monocotyledoneous (grasses) plants. We discuss the biosynthesis and transport of UDP-rhamnose, as well as the transfer of rhamnose from UDP-rhamnose to various primary and specialized metabolites. This is carried out by rhamnosyltransferases, enzymes that can use a large variety of substrates. Some unique characteristics of rhamnose synthases, the multifunctional enzymes responsible for the conversion of UDP-glucose into UDP-rhamnose, are considered, particularly from the perspective of their ability to convert glucose present in flavonoids. Finally, we discuss how little is still known with regards to how plants rescue rhamnose from the many compounds to which it is linked, or how rhamnose is catabolized.

CamRegBase: a gene regulation database for the biofuel crop, Camelina sativa

Camelina is an annual oilseed plant from the Brassicaceae family that is gaining momentum as a biofuel winter cover crop. However, a significant limitation in further enhancing its utility as a producer of oils that can be used as biofuels, jet fuels or bio-based products is the absence of a repository for all the gene expression and regulatory information that is being rapidly generated by the community. Here, we provide CamRegBase (https://camregbase.org/) as a one-stop resource to access Camelina information on gene expression and co-expression, transcription factors, lipid associated genes and genome-wide orthologs in the close-relative reference plant Arabidopsis. We envision this as a resource of curated information for users, as well as a repository of new gene regulation information.

Meta Gene Regulatory Networks in Maize Highlight Functionally Relevant Regulatory Interactions

The regulation of gene expression is central to many biological processes. Gene regulatory networks (GRNs) link transcription factors (TFs) to their target genes and represent maps of potential transcriptional regulation. Here, we analyzed a large number of publically available maize (Zea mays) transcriptome data sets including >6000 RNA sequencing samples to generate 45 coexpression-based GRNs that represent potential regulatory relationships between TFs and other genes in different populations of samples (cross-tissue, cross-genotype, and tissue-and-genotype samples). While these networks are all enriched for biologically relevant interactions, different networks capture distinct TF-target associations and biological processes. By examining the power of our coexpression-based GRNs to accurately predict covarying TF-target relationships in natural variation data sets, we found that presence/absence changes rather than quantitative changes in TF gene expression are more likely associated with changes in target gene expression. Integrating information from our TF-target predictions and previous expression quantitative trait loci (eQTL) mapping results provided support for 68 TFs underlying 74 previously identified trans-eQTL hotspots spanning a variety of metabolic pathways. This study highlights the utility of developing multiple GRNs within a species to detect putative regulators of important plant pathways and provides potential targets for breeding or biotechnological applications.

Discovery of modules involved in the biosynthesis and regulation of maize phenolic compounds

Phenolic compounds are among the most diverse and widespread of specialized plant compounds and underly many important agronomic traits. Our comprehensive analysis of the maize genome unraveled new aspects of the genes involved in phenylpropanoid, monolignol, and flavonoid production in this important crop. Remarkably, just 19 genes accounted for 70 % of the overall mRNA accumulation of these genes across 95 tissues, indicating that these are the main contributors to the flux of phenolic metabolites. Eighty genes with intermediate to low expression play minor and more specialized roles. Remaining genes are likely undergoing loss of function or are expressed in limited cell types. Phylogenetic and expression analyses revealed which members of gene families governing metabolic entry and branch points exhibit duplication, subfunctionalization, or loss of function. Co-expression analysis applied to genes in sequential biosynthetic steps revealed that certain isoforms are highly co-expressed and are candidates for metabolic complexes that ensure metabolite delivery to correct cellular compartments. Co-expression of biosynthesis genes with transcription factors discovered connections that provided candidate components for regulatory modules governing this pathway. Our study provides a comprehensive analysis of maize phenylpropanoid related genes, identifies major pathway contributors, and novel candidate enzymatic and regulatory modules of the metabolic network.

Diversity of genetic lesions characterizes new Arabidopsis flavonoid pigment mutant alleles from T-DNA collections

The visual phenotypes afforded by flavonoid pigments have provided invaluable tools for modern genetics. Many Arabidopsis transparent testa (tt) mutants lacking the characteristic proanthocyanidin (PA) seed coat pigmentation and often failing to accumulate anthocyanins in vegetative tissues have been characterized. These mutants have significantly contributed to our understanding of flavonoid biosynthesis, regulation, and transport. A comprehensive screening for tt mutants in available large T-DNA collection lines resulted in the identification of 16 independent lines lacking PAs and anthocyanins, or with seed coat pigmentation clearly distinct from wild type. Segregation analyses and the characterization of second alleles in the genes disrupted by the indexed T-DNA insertions demonstrated that all the lines contained at least one additional mutation responsible for the tt phenotypes. Using a combination of RNA-Seq and whole genome re-sequencing and confirmed through complementation, we show here that these mutations correspond to novel alleles of ttg1 (two alleles), tt3 (two alleles), tt5 (two alleles), ban (two alleles), tt1 (two alleles), and tt8 (six alleles), which harbored additional T-DNA insertions, indels, missense mutations, and large genomic deletion. Several of the identified alleles offer interesting perspectives on flavonoid biosynthesis and regulation.

Challenges of Translating Gene Regulatory Information into Agronomic Improvements

Improvement of agricultural species has exploited the genetic variation responsible for complex quantitative traits. Much of the functional variation is regulatory, in cis-regulatory elements and trans-acting factors that ultimately contribute to gene expression differences. However, the identification of gene regulatory network components that, when modulated, will increase plant productivity or resilience, is challenging, yet essential to provide increased predictive power for genome engineering approaches that are likely to benefit useful traits. Here, we discuss the opportunities and limitations of using data obtained from gene coexpression, transcription factor binding, and genome-wide association mapping analyses to predict regulatory interactions that impact crop improvement. It is apparent that a combination of information from these data types is necessary for the reliable identification and utilization of important regulatory interactions that underlie complex agronomic traits.

Arabidopsis EMSY-like (EML) histone readers are necessary for post-fertilization seed development, but prevent fertilization-independent seed formation

Seed development is a complex regulatory process that includes a transition from gametophytic to sporophytic program. The synchronized development of different seed compartments (seed coat, endosperm and embryo) depends on a balance in parental genome contributions. In the most severe cases, the disruption of the balance leads to seed abortion. This represents one of the main obstacles for the engineering of asexual reproduction through seeds (apomixis), and for generating new interspecies hybrids. The repression of auxin synthesis by the Polycomb Repressive Complex 2 (PRC2) is a major mechanism contributing to sensing genome balance. However, current efforts focusing on decreasing PRC2 or elevating auxin levels have not yet resulted in the production of apomictic seed. Here, we show that EMSY-Like Tudor/Agenet H3K36me3 histone readers EML1 and EML3 are necessary for early stages of seed development to proceed at a normal rate in Arabidopsis. We further demonstrate that both EML1 and EML3 are required to prevent the initiation of seed development in the absence of fertilization. Based on the whole genome expression analysis, we identify auxin transport and signaling genes as the most enriched downstream targets of EML1 and EML3. We hypothesize that EML1 and EML3 function to repress and soften paternal gene expression by fine-tuning auxin responses. Discovery of this pathway may contribute to the engineering of apomixis and interspecies hybrids.

Apigenin produced by maize flavone synthase I and II protects plants against UV-B-induced damage

Flavones, one of the largest groups of flavonoids, have beneficial effects on human health and are considered of high nutritional value. Previously, we demonstrated that maize type I flavone synthase (ZmFNSI) is one of the enzymes responsible for the synthesis of O-glycosyl flavones in floral tissues. However, in related species such as rice and sorghum, type II FNS enzymes also contribute to flavone biosynthesis. In this work, we provide evidence that maize has both one FNSI and one FNSII flavone synthases. Arabidopsis transgenic plants expressing each FNS enzyme were generated to validate the role of flavones in protecting plants against UV-B radiation. Here, we demostrate that ZmCYP93G7 (FNSII) has flavone synthase activity and is able to complement the Arabidopsis dmr6 mutant, restoring the susceptibility to Pseudomonas syringae. ZmFNSII expression is controlled by the C1/PL1 + R/B anthocyanin transcriptional complexes, and both ZmFNSI and ZmFNSII are regulated by UV-B. Arabidopsis transgenic plants expressing ZmFNSI or ZmFNSII that accumulate apigenin exhibit less UV-B-induced damage than wild-type plants. Together, we show that maize has two FNS-type enzymes that participate in the synthesis of apigenin, conferring protection against UV-B radiation.

Tomato floral induction and flower development are orchestrated by the interplay between gibberellin and two unrelated microRNA-controlled modules

Age-regulated microRNA156 (miR156) and targets similarly control the competence to flower in diverse species. By contrast, the diterpene hormone gibberellin (GA) and the microRNA319-regulated TEOSINTE BRANCHED/CYCLOIDEA/PCF (TCP) transcription factors promote flowering in the facultative long-day Arabidopsis thaliana, but suppress it in the day-neutral tomato (Solanum lycopersicum). We combined genetic and molecular studies and described a new interplay between GA and two unrelated miRNA-associated pathways that modulates tomato transition to flowering. Tomato PROCERA/DELLA activity is required to promote flowering along with the miR156-targeted SQUAMOSA PROMOTER BINDING-LIKE (SPL/SBP) transcription factors by activating SINGLE FLOWER TRUSS (SFT) in the leaves and the MADS-Box gene APETALA1(AP1)/MC at the shoot apex. Conversely, miR319-targeted LANCEOLATE represses floral transition by increasing GA concentrations and inactivating SFT in the leaves and AP1/MC at the shoot apex. Importantly, the combination of high GA concentrations/responses with the loss of SPL/SPB function impaired canonical meristem maturation and flower initiation in tomato. Our results reveal a cooperative regulation of tomato floral induction and flower development, integrating age cues (miR156 module) with GA responses and miR319-controlled pathways. Importantly, this study contributes to elucidate the mechanisms underlying the effects of GA in controlling flowering time in a day-neutral species.

Arabidopsis JMJD5/JMJ30 Acts Independently of LUX ARRHYTHMO Within the Plant Circadian Clock to Enable Temperature Compensation

The circadian system ensures that plants respond appropriately to environmental change by predicting regular transitions that occur during diel cycles. In order to be most useful, the circadian system needs to be compensated against daily and seasonal changes in temperature that would otherwise alter the pace of this biological oscillator. We demonstrate that an evening-phased protein, the putative histone demethylase JMJD5, contributes to temperature compensation. JMJD5 is co-expressed with components of the Evening Complex, an agglomeration of proteins including EARLY FLOWERING3 (ELF3), ELF4, and LUX ARRHYTHYMO (LUX), which also integrates temperature changes into the molecular clockwork. One role of the Evening Complex is to regulate expression of PSEUDORESPONSE REGULATOR9 (PRR9) and PRR7, important components of the temperature compensation mechanism. Surprisingly we find that LUX, but not other Evening Complex components, is dispensable for clock function at low temperatures. Further genetic analysis suggests JMJD5 acts in a parallel pathway to LUX within the circadian system. Although an intact JMJD5 catalytic domain is required for its function within the clock, our findings suggest JMJD5 does not directly regulate H3K36 methylation at circadian loci. Such data refine our understanding of how JMDJ5 acts within the Arabidopsis circadian system.

Beyond the wall: High-throughput quantification of plant soluble and cell-wall bound phenolics by liquid chromatography tandem mass spectrometry

Plants accumulate several thousand of phenolic compounds, including lignins and flavonoids, which are mainly synthesized through the phenylpropanoid pathway, and play important roles in plant growth and adaptation. A novel high-throughput ultra-high performance liquid chromatography tandem mass spectrometry (UHPLC-MS/MS) method was established to quantify the levels of 19 flavonoids and 15 other phenolic compounds, including acids, aldehydes, and alcohols. The chromatographic separation was performed in 10 min, allowing for the resolution of isomers such as 3-, 4-, and 5-chlorogenic acids, 4-hydroxybenzoic and salicylic acids, isoorientin and orientin, and luteolin and kaempferol. The linearity range for each compound was found to be in the low fmol to the high pmol. Furthermore, this UHPLC-MS/MS approach was shown to be very sensitive with limits of detection between 1.5 amol to 300 fmol, and limits of quantification between 5 amol to 1000 fmol. Extracts from maize seedlings were used to assess the robustness of the method in terms of recovery efficiency, matrix effect, and accuracy. The biological matrix did not suppress the signal for 32 out of the 34 metabolites under investigation. Additionally, the majority of the analytes were recovered from the biological samples with an efficiency above 75%. All flavonoids and other phenolic compounds had an intra- and inter-day accuracy within a ±20% range, except for coniferyl alcohol and vanillic acid. Finally, the quantification of flavonoids, free and cell wall-bound phenolics in seedlings from two maize lines with contrasting phenolic content was successfully achieved using this methodology.

Arabidopsis MATE45 antagonizes local abscisic acid signaling to mediate development and abiotic stress responses

Anthocyanins provide ideal visual markers for the identification of mutations that disrupt molecular responses to abiotic stress. We screened Arabidopsis mutants of ABC (ATP-Binding Cassette) and MATE (Multidrug And Toxic compound Extrusion) transporter genes under nutritional stress and identified four genes (ABCG25,ABCG9,ABCG5, and MATE45) required for normal anthocyanin pigmentation. ABCG25 was previously demonstrated to encode a vascular-localized cellular exporter of abscisic acid (ABA). Our results show that MATE45 encodes an aerial meristem- and a vascular-localized transporter associated with the trans-Golgi, and that it plays an important role in controlling the levels and distribution of ABA in growing aerial meristems and non-meristematic tissues. MATE45 promoter-GUS reporter fusions revealed the activity localized to the leaf and influorescence meristems and the vasculature. Loss-of-function mate45 mutants exhibited accelerated rates of aerial organ initiation suggesting at least partial functional conservation with the maize ortholog bige1. The aba2-1 mutant, which is deficient in ABA biosynthesis, exhibited a number of phenotypes that were rescued in the mate45-1 aba2-1 double mutant. mate45 exhibited enhanced the seed dormancy, and germination was hypersensitive to ABA. Enhanced frequency of leaf primordia growth in mate45 seedlings grown in nutrient imbalance stress was ABA-dependent. The ABA signaling reporter construct pRD29B::GUS revealed elevated levels of ABA signaling in the true leaf primordia of mate45 seedlings grown under nutritional stress, and gradually reduced signaling in surrounding cotyledon and hypocotyl tissues concomitant with reduced expressions of ABCG25. Our results suggest a role of MATE45 in reducing meristematic ABA and in maintaining ABA distribution in adjacent non-meristematic tissues.

Topological and statistical analyses of gene regulatory networks reveal unifying yet quantitatively different emergent properties

Understanding complexity in physical, biological, social and information systems is predicated on describing interactions amongst different components. Advances in genomics are facilitating the high-throughput identification of molecular interactions, and graphs are emerging as indispensable tools in explaining how the connections in the network drive organismal phenotypic plasticity. Here, we describe the architectural organization and associated emergent topological properties of gene regulatory networks (GRNs) that describe protein-DNA interactions (PDIs) in several model eukaryotes. By analyzing GRN connectivity, our results show that the anticipated scale-free network architectures are characterized by organism-specific power law scaling exponents. These exponents are independent of the fraction of the GRN experimentally sampled, enabling prediction of properties of the complete GRN for an organism. We further demonstrate that the exponents describe inequalities in transcription factor (TF)-target gene recognition across GRNs. These observations have the important biological implication that they predict the existence of an intrinsic organism-specific trans and/or cis regulatory landscape that constrains GRN topologies. Consequently, architectural GRN organization drives not only phenotypic plasticity within a species, but is also likely implicated in species-specific phenotype.

Imaging Vacuolar Anthocyanins with Fluorescence Lifetime Microscopy (FLIM)

Anthocyanins are intrinsically fluorescent pigments that accumulate in plant vacuoles. We have developed a platform to analyze the fluorescence decay of anthocyanins by fluorescence lifetime imaging microscopy (FLIM), under in vitro and in vivo conditions. Fluorescence lifetime of a fluorophore can be influenced by temperature, pH, oxygen concentration, and other environmental conditions. Within plant cells, the anthocyanin fluorescence lifetime correlates with distinct subcellular compartments. Vacuolar anthocyanins exhibit shorter fluorescence lifetime than the cytoplasmic pool. Consistent with these observations, lower pH of anthocyanins solutions correlated with shorter fluorescence lifetimes. We discuss here the use of FLIM as a tool for analyzing the subcellular distribution of anthocyanins and estimating variation in vacuolar pH in intact cells.

The Next Generation of Training for Arabidopsis Researchers: Bioinformatics and Quantitative Biology

Training for experimental plant biologists needs to combine bioinformatics, quantitative approaches, computational biology, and training in the art of collaboration, best achieved through fully integrated curriculum development.

Standardized Method for High-throughput Sterilization of Arabidopsis Seeds

Arabidopsis thaliana (Arabidopsis) seedlings often need to be grown on sterile media. This requires prior seed sterilization to prevent the growth of microbial contaminants present on the seed surface. Currently, Arabidopsis seeds are sterilized using two distinct sterilization techniques in conditions that differ slightly between labs and have not been standardized, often resulting in only partially effective sterilization or in excessive seed mortality. Most of these methods are also not easily scalable to a large number of seed lines of diverse genotypes. As technologies for high-throughput analysis of Arabidopsis continue to proliferate, standardized techniques for sterilizing large numbers of seeds of different genotypes are becoming essential for conducting these types of experiments. The response of a number of Arabidopsis lines to two different sterilization techniques was evaluated based on seed germination rate and the level of seed contamination with microbes and other pathogens. The treatments included different concentrations of sterilizing agents and times of exposure, combined to determine optimal conditions for Arabidopsis seed sterilization. Optimized protocols have been developed for two different sterilization methods: bleach (liquid-phase) and chlorine (Cl₂) gas (vapor-phase), both resulting in high seed germination rates and minimal microbial contamination. The utility of these protocols was illustrated through the testing of both wild type and mutant seeds with a range of germination potentials. Our results show that seeds can be effectively sterilized using either method without excessive seed mortality, although detrimental effects of sterilization were observed for seeds with lower than optimal germination potential. In addition, an equation was developed to enable researchers to apply the standardized chlorine gas sterilization conditions to airtight containers of different sizes. The protocols described here allow easy, efficient, and inexpensive seed sterilization for a large number of Arabidopsis lines.

The challenges faced by living stock collections in the USA

Many discoveries in the life sciences have been made using material from living stock collections. These collections provide a uniform and stable supply of living organisms and related materials that enhance the reproducibility of research and minimize the need for repetitive calibration. While collections differ in many ways, they all require expertise in maintaining living organisms and good logistical systems for keeping track of stocks and fulfilling requests for specimens. Here, we review some of the contributions made by living stock collections to research across all branches of the tree of life, and outline the challenges they face.

DOI:http://dx.doi.org/10.7554/eLife.24611.001

A Maize Gene Regulatory Network for Phenolic Metabolism

The translation of the genotype into phenotype, represented for example by the expression of genes encoding enzymes required for the biosynthesis of phytochemicals that are important for interaction of plants with the environment, is largely carried out by transcription factors (TFs) that recognize specific cis-regulatory elements in the genes that they control. TFs and their target genes are organized in gene regulatory networks (GRNs), and thus uncovering GRN architecture presents an important biological challenge necessary to explain gene regulation. Linking TFs to the genes they control, central to understanding GRNs, can be carried out using gene- or TF-centered approaches. In this study, we employed a gene-centered approach utilizing the yeast one-hybrid assay to generate a network of protein-DNA interactions that participate in the transcriptional control of genes involved in the biosynthesis of maize phenolic compounds including general phenylpropanoids, lignins, and flavonoids. We identified 1100 protein-DNA interactions involving 54 phenolic gene promoters and 568 TFs. A set of 11 TFs recognized 10 or more promoters, suggesting a role in coordinating pathway gene expression. The integration of the gene-centered network with information derived from TF-centered approaches provides a foundation for a phenolics GRN characterized by interlaced feed-forward loops that link developmental regulators with biosynthetic genes.

WIND1 Promotes Shoot Regeneration through Transcriptional Activation of ENHANCER OF SHOOT REGENERATION1 in Arabidopsis

Many plant species display remarkable developmental plasticity and regenerate new organs after injury. Local signals produced by wounding are thought to trigger organ regeneration but molecular mechanisms underlying this control remain largely unknown. We previously identified an AP2/ERF transcription factor WOUND INDUCED DEDIFFERENTIATION1 (WIND1) as a central regulator of wound-induced cellular reprogramming in plants. In this study, we demonstrate that WIND1 promotes callus formation and shoot regeneration by upregulating the expression of the ENHANCER OF SHOOT REGENERATION1 (ESR1) gene, which encodes another AP2/ERF transcription factor in Arabidopsis thaliana The esr1 mutants are defective in callus formation and shoot regeneration; conversely, its overexpression promotes both of these processes, indicating that ESR1 functions as a critical driver of cellular reprogramming. Our data show that WIND1 directly binds the vascular system-specific and wound-responsive cis-element-like motifs within the ESR1 promoter and activates its expression. The expression of ESR1 is strongly reduced in WIND1-SRDX dominant repressors, and ectopic overexpression of ESR1 bypasses defects in callus formation and shoot regeneration in WIND1-SRDX plants, supporting the notion that ESR1 acts downstream of WIND1. Together, our findings uncover a key molecular pathway that links wound signaling to shoot regeneration in plants.

Using fluorescence lifetime microscopy to study the subcellular localization of anthocyanins

Anthocyanins are flavonoid pigments that accumulate in most seed plants. They are synthesized in the cytoplasm but accumulate inside the vacuoles. Anthocyanins are pigmented at the lower vacuolar pH, but in the cytoplasm they can be visualized based on their fluorescence properties. Thus, anthocyanins provide an ideal system for the development of new methods to investigate cytoplasmic pools and association with other molecular components. We have analyzed the fluorescence decay of anthocyanins by fluorescence lifetime imaging microscopy (FLIM), in both in vitro and in vivo conditions, using wild-type and mutant Arabidopsis thaliana seedlings. Within plant cells, the amplitude-weighted mean fluorescence lifetime (τ_m ) correlated with distinct subcellular localizations of anthocyanins. The vacuolar pool of anthocyanins exhibited shorter τ_m than the cytoplasmic pool. Consistently, lowering the pH of anthocyanins in solution shortened their fluorescence decay. We propose that FLIM is a useful tool for understanding the trafficking of anthocyanins and, potentially, for estimating vacuolar pH inside intact plant cells.

Combinatorial control of plant gene expression

Combinatorial gene regulation provides a mechanism by which relatively small numbers of transcription factors can control the expression of a much larger number of genes with finely tuned temporal and spatial patterns. This is achieved by transcription factors assembling into complexes in a combinatorial fashion, exponentially increasing the number of genes that they can target. Such an arrangement also increases the specificity and affinity for the cis-regulatory sequences required for accurate target gene expression. Superimposed on this transcription factor combinatorial arrangement is the increasing realization that histone modification marks expand the regulatory information, which is interpreted by histone readers and writers that are part of the regulatory apparatus. Here, we review the progress in these areas from the perspective of plant combinatorial gene regulation, providing examples of different regulatory solutions and comparing them to other metazoans. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Identification and Characterization of Maize salmon silks Genes Involved in Insecticidal Maysin Biosynthesis

The century-old maize (Zea mays) salmon silks mutation has been linked to the absence of maysin. Maysin is a C-glycosyl flavone that, when present in silks, confers natural resistance to the maize earworm (Helicoverpa zea), which is one of the most damaging pests of maize in America. Previous genetic analyses predicted Pericarp Color1 (P1; R2R3-MYB transcription factor) to be epistatic to the sm mutation. Subsequent studies identified two loci as being capable of conferring salmon silks phenotypes, salmon silks1 (sm1) and sm2 Benefitting from available sm1 and sm2 mapping information and from knowledge of the genes regulated by P1, we describe here the molecular identification of the Sm1 and Sm2 gene products. Sm2 encodes a rhamnosyl transferase (UGT91L1) that uses isoorientin and UDP-rhamnose as substrates and converts them to rhamnosylisoorientin. Sm1 encodes a multidomain UDP-rhamnose synthase (RHS1) that converts UDP-glucose into UDP-l-rhamnose. Here, we demonstrate that RHS1 shows unexpected substrate plasticity in converting the glucose moiety in rhamnosylisoorientin to 4-keto-6-deoxy glucose, resulting in maysin. Both Sm1 and Sm2 are direct targets of P1, as demonstrated by chromatin immunoprecipitation experiments. The molecular characterization of Sm1 and Sm2 described here completes the maysin biosynthetic pathway, providing powerful tools for engineering tolerance to maize earworm in maize and other plants.

Dietary Apigenin Exerts Immune-Regulatory Activity in Vivo by Reducing NF-κB Activity, Halting Leukocyte Infiltration and Restoring Normal Metabolic Function

The increasing prevalence of inflammatory diseases and the adverse effects associated with the long-term use of current anti-inflammatory therapies prompt the identification of alternative approaches to reestablish immune balance. Apigenin, an abundant dietary flavonoid, is emerging as a potential regulator of inflammation. Here, we show that apigenin has immune-regulatory activity in vivo. Apigenin conferred survival to mice treated with a lethal dose of Lipopolysaccharide (LPS) restoring normal cardiac function and heart mitochondrial Complex I activity. Despite the adverse effects associated with high levels of splenocyte apoptosis in septic models, apigenin had no effect on reducing cell death. However, we found that apigenin decreased LPS-induced apoptosis in lungs, infiltration of inflammatory cells and chemotactic factors’ accumulation, re-establishing normal lung architecture. Using NF-κB luciferase transgenic mice, we found that apigenin effectively modulated NF-κB activity in the lungs, suggesting the ability of dietary compounds to exert immune-regulatory activity in an organ-specific manner. Collectively, these findings provide novel insights into the underlying immune-regulatory mechanisms of dietary nutraceuticals in vivo.

50 years of Arabidopsis research: highlights and future directions

922 I. 922 II. 922 III. 925 IV. 925 V. 926 VI. 927 VII. 928 VIII. 929 IX. 930 X. 931 XI. 932 XII. 933 XIII. Natural variation and genome-wide association studies 934 XIV. 934 XV. 935 XVI. 936 XVII. 937 937 References 937 SUMMARY: The year 2014 marked the 25(th) International Conference on Arabidopsis Research. In the 50 yr since the first International Conference on Arabidopsis Research, held in 1965 in Göttingen, Germany, > 54 000 papers that mention Arabidopsis thaliana in the title, abstract or keywords have been published. We present herein a citational network analysis of these papers, and touch on some of the important discoveries in plant biology that have been made in this powerful model system, and highlight how these discoveries have then had an impact in crop species. We also look to the future, highlighting some outstanding questions that can be readily addressed in Arabidopsis. Topics that are discussed include Arabidopsis reverse genetic resources, stock centers, databases and online tools, cell biology, development, hormones, plant immunity, signaling in response to abiotic stress, transporters, biosynthesis of cells walls and macromolecules such as starch and lipids, epigenetics and epigenomics, genome-wide association studies and natural variation, gene regulatory networks, modeling and systems biology, and synthetic biology.