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.

Molecular mechanisms underlying gene regulatory variation of maize metabolic traits.

Variation in gene expression levels is pervasive among individuals and races or varieties, and has substantial agronomic consequences, for example, by contributing to hybrid vigor. Gene expression level variation results from mutations in regulatory sequences (cis) and/or transcription factor (TF) activity (trans), but the mechanisms underlying cis- and/or trans-regulatory variation of complex phenotypes remain largely unknown. Here, we investigated gene expression variation mechanisms underlying the differential accumulation of the insecticidal compounds maysin and chlorogenic acid in silks of widely used maize (Zea mays) inbreds, B73 and A632. By combining transcriptomics and cistromics, we identified 1,338 silk direct targets of the maize R2R3-MYB TF Pericarp color1 (P1), consistent with it being a regulator of maysin and chlorogenic acid biosynthesis. Among these P1 targets, 464 showed allele-specific expression (ASE) between B73 and A632 silks. Allelic DNA-affinity purification sequencing identified 34 examples in which P1 allelic specific binding (ASB) correlated with cis-expression variation. From previous yeast one-hybrid studies, we identified 9 TFs potentially implicated in the control of P1 targets, with ASB to 83 out of 464 ASE genes (cis) and differential expression of 4 out of 9 TFs between B73 and A632 silks (trans). These results provide a molecular framework for understanding universal mechanisms underlying natural variation of gene expression levels, and how the regulation of metabolic diversity is established.

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.

Targeted modulation of pennycress lipid droplet proteins impacts droplet morphology and seed oil content.

Lipid droplets (LDs) are unusual organelles that have a phospholipid monolayer surface and a hydrophobic matrix. In oilseeds, this matrix is nearly always composed of triacylglycerols (TGs) for efficient storage of carbon and energy. Various proteins play a role in their assembly, stability and turnover, and even though the major structural oleosin proteins in seed LDs have been known for decades, the factors influencing LD formation and dynamics are still being uncovered mostly in the "model oilseed" Arabidopsis. Here we identified several key LD biogenesis proteins in the seeds of pennycress, a potential biofuel crop, that were correlated previously with seed oil content and characterized here for their participation in LD formation in transient expression assays and stable transgenics. One pennycress protein, the lipid droplet associated protein-interacting protein (LDIP), was able to functionally complement the Arabidopsis ldip mutant, emphasizing the close conservation of lipid storage among these two Brassicas. Moreover, loss-of-function ldip mutants in pennycress exhibited increased seed oil content without compromising plant growth, raising the possibility that LDIP or other LD biogenesis factors may be suitable targets for improving yields in oilseed crops more broadly.

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.

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.

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.

Multi-scale 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 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 multi-scale 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.

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.

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.

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.

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 (RACT) and GL3 (GL3ACT). We show that RACT dimerizes with a dissociation constant around 100 nM, over an order of magnitude stronger than GL3ACT. Structural predictions combined with mutational analyses demonstrated that V568, located in a hydrophobic pocket in RACT, is important: when mutated to the Ser residue present in GL3ACT, dimerization affinity dropped by almost an order of magnitude. The converse S595V mutation in GL3ACT 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.

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.

Unraveling the KNOTTED1 regulatory network in maize meristems.

KNOTTED1 (KN1)-like homeobox (KNOX) transcription factors function in plant meristems, self-renewing structures consisting of stem cells and their immediate daughters. We defined the KN1 cistrome in maize inflorescences and found that KN1 binds to several thousand loci, including 643 genes that are modulated in one or multiple tissues. These KN1 direct targets are strongly enriched for transcription factors (including other homeobox genes) and genes participating in hormonal pathways, most significantly auxin, demonstrating that KN1 plays a key role in orchestrating the upper levels of a hierarchical gene regulatory network that impacts plant meristem identity and function.

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.

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.

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.

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.

Core Promoter Plasticity Between Maize Tissues and Genotypes Contrasts with Predominance of Sharp Transcription Initiation Sites.

Core promoters are crucial for gene regulation, providing blueprints for the assembly of transcriptional machinery at transcription start sites (TSSs). Empirically, TSSs define the coordinates of core promoters and other regulatory sequences. Thus, experimental TSS identification provides an essential step in the characterization of promoters and their features. Here, we describe the application of CAGE (cap analysis of gene expression) to identify genome-wide TSSs used in root and shoot tissues of two maize (Zea mays) inbred lines (B73 and Mo17). Our studies indicate that most TSS clusters are sharp in maize, similar to mice, but distinct from Arabidopsis thaliana, Drosophila melanogaster, or zebra fish, in which a majority of genes have broad-shaped TSS clusters. We established that ∼38% of maize promoters are characterized by a broader TATA-motif consensus, and this motif is significantly enriched in genes with sharp TSSs. A noteworthy plasticity in TSS usage between tissues and inbreds was uncovered, with ∼1500 genes showing significantly different dominant TSSs, sometimes affecting protein sequence by providing alternate translation initiation codons. We experimentally characterized instances in which this differential TSS utilization results in protein isoforms with additional domains or targeted to distinct subcellular compartments. These results provide important insights into TSS selection and gene expression in an agronomically important crop.