Transcriptional regulation in plants: the importance of combinatorial control

Genome-wide analysis of bromodomain gene family in Arabidopsis and rice

The bromodomain-containing proteins (BRD-proteins) belongs to family of 'epigenetic mark readers', integral to epigenetic regulation. The BRD-members contain a conserved 'bromodomain' (BRD/BRD-fold: interacts with acetylated-lysine in histones), and several additional domains, making them structurally/functionally diverse. Like animals, plants also contain multiple Brd-homologs, however the extent of their diversity and impact of molecular events (genomic duplications, alternative splicing, AS) therein, is relatively less explored. The present genome-wide analysis of Brd-gene families of Arabidopsis thaliana and Oryza sativa showed extensive diversity in structure of genes/proteins, regulatory elements, expression pattern, domains/motifs, and the bromodomain (w.r.t. length, sequence, location) among the Brd-members. Orthology analysis identified thirteen ortholog groups (OGs), three paralog groups (PGs) and four singleton members (STs). While more than 40% Brd-genes were affected by genomic duplication events in both plants, AS-events affected 60% A. thaliana and 41% O. sativa genes. These molecular events affected various regions (promoters, untranslated regions, exons) of different Brd-members with potential impact on expression and/or structure-function characteristics. RNA-Seq data analysis indicated differences in tissue-specificity and stress response of Brd-members. Analysis by RT-qPCR revealed differential abundance and salt stress response of duplicate A. thaliana and O. sativa Brd-genes. Further analysis of AtBrd gene, AtBrdPG1b showed salinity-induced modulation of splicing pattern. Bromodomain (BRD)-region based phylogenetic analysis placed the A. thaliana and O. sativa homologs into clusters/sub-clusters, mostly consistent with ortholog/paralog groups. The bromodomain-region displayed several conserved signatures in key BRD-fold elements (α-helices, loops), along with variations (1-20 sites) and indels among the BRD-duplicates. Homology modeling and superposition identified structural variations in BRD-folds of divergent and duplicate BRD-members, which might affect their interaction with the chromatin histones, and associated functions. The study also showed contribution of various duplication events in Brd-gene family expansion among diverse plants, including several monocot and dicot plant species.

Transcriptome analysis and exploration of genes involved in the biosynthesis of secoiridoids in Gentiana rhodantha

Gentiana rhodantha is a medicinally important perennial herb used as traditional Chinese and ethnic medicines. Secoiridoids are one of the major bioactive compounds in G. rhodantha. To better understand the secoiridoid biosynthesis pathway, we generated transcriptome sequences from four organs (root, leaf, stem and flower), followed by the de novo sequence assembly. We verified 8-HGO (8-hydroxygeraniol oxidoreductase), which may encode key enzymes of the secoiridoid biosynthesis by qRT-PCR. The mangiferin, swertiamarin and loganic acid contents in root, stem, leaf, and flower were determined by HPLC. The results showed that there were 47,871 unigenes with an average length of 1,107.38 bp. Among them, 1,422 unigenes were involved in 25 standard secondary metabolism-related pathways in the KEGG database. Furthermore, we found that 1,005 unigenes can be divided into 66 transcription factor (TF) families, with no family members exhibiting significant organ-specificity. There were 54 unigenes in G. rhodantha that encoded 17 key enzymes of the secoiridoid biosynthetic pathway. The qRT-PCR of the 8-HGO and HPLC results showed that the relative expression and the mangiferin, swertiamarin, and loganic acid contents of the aerial parts were higher than in the root. Six types of SSR were identified by SSR analysis of unigenes: mono-nucleoside repeat SSR, di-nucleoside repeat SSR, tri-nucleoside repeat SSR, tetra-nucleoside repeat SSR, penta-nucleoside repeat SSR, and hexa-nucleoside repeat SSR. This report not only enriches the Gentiana transcriptome database but helps further study the function and regulation of active component biosynthesis of G. rhodantha.

Uroprotective and Hepatoprotective Potential of Anagallis arvensis against the Experimental Animal Model

Anagallis arvensis (A. arvensis) belonging to the family Primulaceae is traditionally used for liver and kidney diseases. The aim of the study was to evaluate the uroprotective and hepatoprotective potentials of A. arvensis in cyclophosphamide-induced interstitial cystitis and paracetamol-induced hepatotoxicity rat model, respectively. Nociception, bladder weight, vesical vascular permeability, Gray's criteria for edema and hemorrhage, and levels of nitric oxide, catalase, and glutathione were estimated and studied in the cystitis model. Liver function test, lipid profile, and histopathological evaluation were carried out in the hepatoprotective activity. Oral administration of methanol extract of A. arvensis significantly reduced bladder weight, vesical vascular permeability, edema, hemorrhage, nitric oxide, IL-6, and TNF-α, while the level of catalase and glutathione peroxide was increased. In hepatoprotective activity, pretreatment with A. arvensis significantly decreased the level of liver markers (Bilirubin, ALT, AST, and ALP) and lipid profile (cholesterol, TG, LDL, and VLDL). Histopathological studies confirmed the biochemical findings of both studies. GC-MS analysis presented the presence of antioxidant phytoconstituents. Thus, it was concluded that A. arvensis might act as uroprotective and hepatoprotective due to the presence of antioxidant phytochemicals in the rodent model. Isolation and identification of phytochemicals present in the methanol extract of A. arvensis and evaluation of their exact mechanism of action become mandatory in future studies.

An overview of bioinformatics, genomics, and transcriptomics resources for bryophytes

Bryophytes are useful models for the study of plant evolution, development, plant-fungal symbiosis, stress responses, and gametogenesis. Additionally, their dominant haploid gametophytic phase makes them great models for functional genomics research, allowing straightforward genome editing and gene knockout via CRISPR or homologous recombination. Until 2016, however, the only bryophyte genome sequence published was that of Physcomitrium patens. Throughout recent years, several other bryophyte genomes and transcriptome datasets became available, enabling better comparative genomics in evolutionary studies. The increase in the number of bryophyte genome and transcriptome resources available has yielded a plethora of annotations, databases, and bioinformatics tools to access the new data, which covers the large diversity of this clade and whose biology comprises features such as association with arbuscular mycorrhiza fungi, sex chromosomes, low gene redundancy, or loss of RNA editing genes for organellar transcripts. Here we provide a guide to resources available for bryophytes with regards to genome and transcriptome databases and bioinformatics tools.

Genome-wide analysis of AAAG and ACGT cis-elements in Arabidopsis thaliana reveals their involvement with genes downregulated under jasmonic acid response in an orientation independent manner

Cis-regulatory elements are regions of noncoding DNA that regulate the transcription of neighboring genes. The study of cis-element architecture that functions in transcription regulation are essential. AAAG and ACGT are a class of cis-regulatory elements, known to interact with Dof and bZIP transcription factors respectively, and are known to regulate the expression of auxin response, gibberellin response, floral development, light response, seed storage proteins genes, biotic and abiotic stress genes in plants. Analysis of the frequency of occurrence of AAAG and ACGT motifs from varying spacer lengths (0–30 base pair) between these 2 motifs in both possible orientations—AAAG _(N) ACGT and ACGT _(N) AAAG, in the promoters and genome of Arabidopsis thaliana which indicated preferred orientation of AAAG _(N) ACGT over ACGT _(N) AAAG across the genome and in promoters. Further, microarray analysis revealed the involvement of these motifs in the genes downregulated under jasmonic acid response in an orientation-independent manner. These results were further confirmed by the transient expression studies with promoter-reporter cassettes carrying AAAG and ACGT motifs in both orientations. Furthermore, cluster analysis on genes with AAAG _(N) ACGT and ACGT _(N) AAAG motifs orientations revealed clusters of genes to be involved in ABA signaling, transcriptional regulation, DNA binding, and metal ion binding. These findings can be utilized in designing synthetic promoters for the development of stress-tolerant transgenic plants and also provides an insight into the roles of these motifs in transcriptional regulation.

Cis Elements: Added Boost to the Directed Evolution of Plant Genes

To increase the expression of a native/foreign plant/bacterial gene, the complete network of cis-elements must be excavated to increase its biosynthetic yield, especially under industrial stress conditions. For selecting the best set of cis-elements for a foreign gene and aiding the workflow of researchers, often untrained in bioinformatics methodologies, we developed a modular PERL script for their identification and localization. The script is functional on any operating system. It localizes the cis element network of a gene. It aids an easy customization, as per the required analysis, and provides robust strategy, unlike the usually used databases where several applied calculations often become a tricky task. The script allows an uncomplicated analysis of multiplicity of cis elements along with their relative distances, making it easier for designing the more beneficial network of genes for directed evolution experiments. Through a batched scrutiny of several functionally similar genes, it would aid an easy extraction of their evolutionarily favored network of cis elements. It would be extremely helpful to develop the crop plants that are better adapted to the stressful conditions.

Transcription Factors Interact with ABA through Gene Expression and Signaling Pathways to Mitigate Drought and Salinity Stress

Among abiotic stressors, drought and salinity seriously affect crop growth worldwide. In plants, research has aimed to increase stress-responsive protein synthesis upstream or downstream of the various transcription factors (TFs) that alleviate drought and salinity stress. TFs play diverse roles in controlling gene expression in plants, which is necessary to regulate biological processes, such as development and environmental stress responses. In general, plant responses to different stress conditions may be either abscisic acid (ABA)-dependent or ABA-independent. A detailed understanding of how TF pathways and ABA interact to cause stress responses is essential to improve tolerance to drought and salinity stress. Despite previous progress, more active approaches based on TFs are the current focus. Therefore, the present review emphasizes the recent advancements in complex cascades of gene expression during drought and salinity responses, especially identifying the specificity and crosstalk in ABA-dependent and -independent signaling pathways. This review also highlights the transcriptional regulation of gene expression governed by various key TF pathways, including AP2/ERF, bHLH, bZIP, DREB, GATA, HD-Zip, Homeo-box, MADS-box, MYB, NAC, Tri-helix, WHIRLY, WOX, WRKY, YABBY, and zinc finger, operating in ABA-dependent and -independent signaling pathways.

Identification and characterization of long non-coding RNA and their response against citrus bark cracking viroid infection in Humulus lupulus

Long non-coding RNAs (lncRNAs) are a highly heterogeneous class of non-protein-encoding transcripts that play an essential regulatory role in diverse biological processes, including stress responses. The severe stunting disease caused by Citrus bark cracking viroid (CBCVd) poses a major threat to the production of Humulus lupulus (hop) plants. In this study, we systematically investigate the characteristics of the lncRNAs in hop and their role in CBCVd-infection using RNA-sequencing data. Following a stringent filtration criterion, a total of 3598 putative lncRNAs were identified with a high degree of certainty, of which 19% (684) of the lncRNAs were significantly differentially expressed (DE) in CBCVd-infected hop, which were predicted to be mainly involved in plant-pathogen interactions, kinase cascades, secondary metabolism and phytohormone signal transduction. Besides, several lncRNAs and CBCVd-responsive lncRNAs were identified as the precursor of microRNAs and predicted as endogenous target mimics (eTMs) for hop microRNAs involved in CBCVd-infection.

Fine mapping identified the gibberellin 2-oxidase gene CpDw leading to a dwarf phenotype in squash (Cucurbita pepo L.)

Dwarfism is an important agronomic trait in pumpkin that can increase yield. In this study, the dwarf Cucurbita pepo L. line X10 exhibited significantly longitudinally shorter cell length in the stem than did the normal-vine line JIN234. The dwarf stature of X10 was recovered with exogenous gibberellin (GA₃) application, suggesting that X10 might be sensitive to GA biosynthesis. Genetic analysis revealed that this dwarf trait is controlled by a single completely dominant locus: CpDw (Cucurbita pepo L. Dwarf). Using 1,300 F₂ individuals derived from a cross between X10 and JIN234, we mapped the CpDw locus to a region of approximately 24.6 kb on chromosome 10 that contain 5 annotated genes. The high expression level of Cp4.1LG10g05910.1 and high GA2ox enzyme activity in X10 revealed that the GA 2-oxidase gene Cp4.1LG10g05910.1 is a candidate gene for CpDw. Alignment of the Cp4.1LG10g05910.1 gene revealed two nonsynonymous single nucleotide polymorphism (SNP) mutations in the two exons, as well as several SNPs and InDels in the important functional elements of promoter between parental lines. Further allelic diversity analysis of the Cucurbita spp. germplasm resources indicated that Cp4.1LG10g05910.1 may be involved in vine growth during the early developmental stage in C. pepo but not in C. maxima or C. moschata. This study provides an important theoretical basis for the genetic regulation of vine length and crop breeding in pumpkin.

Unravelling the Complex Interplay of Transcription Factors Orchestrating Seed Oil Content in Brassica napus L

Transcription factors (TFs) and their complex interplay are essential for directing specific genetic programs, such as responses to environmental stresses, tissue development, or cell differentiation by regulating gene expression. Knowledge regarding TF-TF cooperations could be promising in gaining insight into the developmental switches between the cultivars of Brassica napus L., namely Zhongshuang11 (ZS11), a double-low accession with high-oil- content, and Zhongyou821 (ZY821), a double-high accession with low-oil-content. In this regard, we analysed a time series RNA-seq data set of seed tissue from both of the cultivars by mainly focusing on the monotonically expressed genes (MEGs). The consideration of the MEGs enables the capturing of multi-stage progression processes that are orchestrated by the cooperative TFs and, thus, facilitates the understanding of the molecular mechanisms determining seed oil content. Our findings show that TF families, such as NAC, MYB, DOF, GATA, and HD-ZIP are highly involved in the seed developmental process. Particularly, their preferential partner choices as well as changes in their gene expression profiles seem to be strongly associated with the differentiation of the oil content between the two cultivars. These findings are essential in enhancing our understanding of the genetic programs in both cultivars and developing novel hypotheses for further experimental studies.

Genome-wide discovery of OsHOX24-binding sites and regulation of desiccation stress response in rice

OsHOX24 mediates regulation of desiccation stress response via complex regulatory network as indicated by its binding to several target genes including transcription factors in rice. HD-ZIP I subfamily of homeobox transcription factors (TFs) are involved in abiotic stress responses and plant development. Previously, we demonstrated the role of OsHOX24, a member of HD-ZIP I subfamily, in abiotic stress responses. In this study, we identified downstream targets of OsHOX24 under control and desiccation stress conditions via chromatin immunoprecipitation-sequencing (ChIP-seq) approach in wild-type and OsHOX24 over-expression transgenic in rice. OsHOX24-binding sites in each sample and differential binding sites between the samples were detected at various genomic locations, including genic and intergenic regions. Gene ontology enrichment analysis revealed that OsHOX24 direct target genes were involved in several biological processes, including plant development, ABA-mediated signalling pathway, ubiquitin-dependent protein catabolic process, ion transport, abiotic and biotic stress responses besides transcriptional and translational regulation. The enrichment of several cis-regulatory motifs representing binding sites of other TFs, such as ABFs, ERF1, MYB1, LTREs and SORLIP2, suggested the involvement of OsHOX24 in a complex regulatory network. These findings indicate that OsHOX24-mediated desiccation stress regulation involves modulation of a plethora of target genes, which participate in diverse pathways in rice.

Genome-wide analysis of PmTCP4 transcription factor binding sites by ChIP-Seq during pistil abortion in Japanese apricot

The TCP4 transcription factor plays an important role in plant growth and development, especially in flower development. PmTCP4 is involved in the process of pistil abortion in Japanese apricot, but its molecular mechanism, particularly the DNA binding sites and co-regulatory genes, are quite unknown. Therefore, to identify the genome-wide binding sites of PmTCP4 transcription factors and their co-regulatory genes, chromatin immunoprecipitation sequencing (ChIP-Seq) was carried out. ChIP-Seq data produced the maximum enriched peaks in two Japanese apricot cultivars 'Daqiandi' (DQD) and 'Longyan' (LY), which showed that the majority of DNA-protein interactions are relevant and have a significant function in binding sites. Moreover, 720 and 251 peak-associated genes regulated by PmTCP4 were identified in DQD and LY, respectively, and most of them were involved in the flower and pistil development process. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that photosynthesis and oxidative phosphorylation were the most enriched pathways in both cultivars and all identified genes related to these pathways were down-regulated. This study will provide a reference for a better understanding of the PmTCP4 regulatory mechanism during pistil abortion in Japanese apricot.

Molecular characterization of a RING E3 ligase SbHCI1 in sorghum under heat and abscisic acid stress

Molecular function ofRING E3 ligase SbHCI1is involved in ABA-mediated basal heat stress tolerancein sorghum. Global warming generally reduces plant survival, owing to the negative effects of high temperatures on plant development. However, little is known about the role of Really Interesting New Gene (RING) E3 ligase in the heat stress responses of plants. As such, the aim of the present study was to characterize the molecular functions of the Sorghum bicolor ortholog of the Oryza sativa gene for Heat- and Cold-Induced RING finger protein 1 (SbHCI1). Subcellular localization revealed that SbHCI1 was mainly associated with the cytosol and that it moved to the Golgi apparatus under heat stress conditions. The fluorescent signals of SbHCI1 substrate proteins were observed to migrate to the cytoplasm under heat stress conditions. Bimolecular fluorescence complementation (BiFC) and yeast two-hybrid (Y2H) assays revealed that SbHCI1 physically interacted with OsHCI1 ortholog partner proteins in the cytoplasm. Moreover, an in vitro ubiquitination assay revealed that SbHCI1 polyubiquitinated each of the three interacting proteins. The ectopic overexpression of SbHCI1 in Arabidopsis revealed that the protein was capable of inducing abscisic acid (ABA)-hypersensitivity and basal heat stress tolerance. Therefore, SbHCI1 possesses E3 ligase activity and may function as a positive regulator of heat stress responses through the modulation of interacting proteins.

Prediction of condition-specific regulatory genes using machine learning

Recent advances in genomic technologies have generated data on large-scale protein-DNA interactions and open chromatin regions for many eukaryotic species. How to identify condition-specific functions of transcription factors using these data has become a major challenge in genomic research. To solve this problem, we have developed a method called ConSReg, which provides a novel approach to integrate regulatory genomic data into predictive machine learning models of key regulatory genes. Using Arabidopsis as a model system, we tested our approach to identify regulatory genes in data sets from single cell gene expression and from abiotic stress treatments. Our results showed that ConSReg accurately predicted transcription factors that regulate differentially expressed genes with an average auROC of 0.84, which is 23.5-25% better than enrichment-based approaches. To further validate the performance of ConSReg, we analyzed an independent data set related to plant nitrogen responses. ConSReg provided better rankings of the correct transcription factors in 61.7% of cases, which is three times better than other plant tools. We applied ConSReg to Arabidopsis single cell RNA-seq data, successfully identifying candidate regulatory genes that control cell wall formation. Our methods provide a new approach to define candidate regulatory genes using integrated genomic data in plants.

Parallel Transcriptional Regulation of Artemisinin and Flavonoid Biosynthesis

Plants regulate the synthesis of specialized compounds through the actions of individual transcription factors (TFs) or sets of TFs. One such compound, artemisinin from Artemisia annua, is widely used as a pharmacological product in the first-line treatment of malaria. However, the emergence of resistance to artemisinin in Plasmodium species, as well as its low production rates, have required innovative treatments such as exploiting the synergistic effects of flavonoids with artemisinin. We overview current knowledge about flavonoid and artemisinin transcriptional regulation in A. annua, and review the dual action of TFs and structural genes that can regulate both pathways simultaneously. Understanding the concerted action of these TFs and their associated structural genes can guide the development of strategies to further improve flavonoid and artemisinin production.

Genomics, molecular and evolutionary perspective of NAC transcription factors

NAC (NAM, ATAF1,2, and CUC2) transcription factors are one of the largest transcription factor families found in the plants and are involved in diverse developmental and signalling events. Despite the availability of comprehensive genomic information from diverse plant species, the basic genomic, biochemical, and evolutionary details of NAC TFs have not been established. Therefore, NAC TFs family proteins from 160 plant species were analyzed in the current study. Study revealed, Brassica napus (410) encodes highest number and Klebsormidium flaccidum (3) encodes the lowest number of TFs. The study further revealed the presence of NAC TF in the Charophyte algae K. flaccidum. On average, the monocot plants encode higher number (141.20) of NAC TFs compared to the eudicots (125.04), gymnosperm (75), and bryophytes (22.66). Furthermore, our analysis revealed that several NAC TFs are membrane bound and contain monopartite, bipartite, and multipartite nuclear localization signals. NAC TFs were also found to encode several novel chimeric proteins and regulate a complex interactome network. In addition to the presence of NAC domain, several NAC proteins were found to encode other functional signature motifs as well. Relative expression analysis of NAC TFs in A. thaliana revealed root tissue treated with urea and ammonia showed higher level of expression and leaf tissues treated with urea showed lower level of expression. The synonymous codon usage is absent in the NAC TFs and it appears that they have evolved from orthologous ancestors and undergone vivid duplications to give rise to paralogous NAC TFs. The presence of novel chimeric NAC TFs are of particular interest and the presence of chimeric NAC domain with other functional signature motifs in the NAC TF might encode novel functional properties in the plants.

A role for Arabidopsis growth-regulating factors 1 and 3 in growth-stress antagonism

Growth-regulating factors (GRFs) belong to a small family of transcription factors that are highly conserved in plants. GRFs regulate many developmental processes and plant responses to biotic and abiotic stimuli. Despite the importance of GRFs, a detailed mechanistic understanding of their regulatory functions is still lacking. In this study, we used ChIP sequencing (ChIP-seq) to identify genome-wide binding sites of Arabidopsis GRF1 and GRF3, and correspondingly their direct downstream target genes. RNA-sequencing (RNA-seq) analysis revealed that GRF1 and GRF3 regulate the expression of a significant number of the identified direct targets. The target genes unveiled broad regulatory functions of GRF1 and GRF3 in plant growth and development, phytohormone biosynthesis and signaling, and the cell cycle. Our analyses also revealed that clock core genes and genes with stress- and defense-related functions are most predominant among the GRF1- and GRF3-bound targets, providing insights into a possible role for these transcription factors in mediating growth-defense antagonism and integrating environmental stimuli into developmental programs. Additionally, GRF1 and GRF3 target molecular nodes of growth-defense antagonism and modulate the levels of defense- and development-related hormones in opposite directions. Taken together, our results point to GRF1 and GRF3 as potential key determinants of plant fitness under stress conditions.

OsNAC2 integrates auxin and cytokinin pathways to modulate rice root development

The rice root system is important for growth. The crosstalk between auxin and cytokinin mediates root initiation and elongation. However, it remains unclear how the transcriptional network upstream of the auxin and cytokinin signalling pathways determines root development. Here, we observed that the knockdown of OsNAC2, which encodes a NAC transcription factor, increased the primary root length and the number of crown roots. OsNAC2 predominantly expressed in primary root tips, crown roots and lateral root primordia, implying it influences root development. Molecular analyses revealed that the expressions of auxin- and cytokinin-responsive genes were affected in OsNAC2-overexpressing (OsNAC2-OX; ON7 and ON11), RNA interference (OsNAC2-RNAi; RNAi25 and RNAi31) and CRISPR/Cas9 plants. Additionally, OsNAC2 can directly bind to the promoters of IAA inactivation-related genes (GH3.6 and GH3.8), an IAA signalling-related gene (OsARF25), and a cytokinin oxidase gene (OsCKX4). Furthermore, genetic analysis of ON11/osgh3.6 and RNAi31/osckx4 homozygote confirmed that OsCKX4 and OsGH3.6 functioned downstream of OsNAC2. The mRNA levels of CROWN ROOTLESS (CRL) genes and cyclin-dependent protein kinase (CDK) genes increased in OsNAC2-RNAi and OsNAC2-cas9 lines while reduced in OsNAC2-OX lines. Thus, we describe that OsNAC2 functions as an upstream integrator of auxin and cytokinin signals that affect CRL and CDK production to regulate cell division during root development. This novel auxin-OsNAC2-cytokinin model should provide a new insight into the understanding of NAC TFs and crosstalk of auxin and cytokinin pathway, and can be potentially applied in agriculture to enhance rice yields by genetic approaches.

Phytocystatins and their Potential Application in the Development of Drought Tolerance Plants in Soybeans (Glycine max L.)

Plant cystatins, also called phytocystatins constitute a family of specific cysteine protease inhibitors found in several monocots and dicots. In soybean, phytocystatins regulate several endogenous processes contributing immensely to this crop's tolerance to abiotic stress factors. Soybeans offer numerous nutritional, pharmaceutical and industrial benefits; however, their growth and yields is hampered by drought, which causes more than 10% yield losses recorded every harvest period worldwide. This review analyses the role of papain-like cysteine proteases and their inhibitors in soybean plant growth and development under drought stress. It also describes their localisation, regulation, target organs and tissues, and the overall impact of cystatins on generating drought tolerance soybean plants. These proteins have many functions that remain poorly characterized, particularly under abiotic stress. Although much information is available on the utilisation of proteases for industrial applications, very few reports have focused on the impact of proteases on plant stress responses. The exploitation of cystatins in plant engineering, as competitive proteases inhibitors is one of the means that will guarantee the continued utilisation of soybeans as an important oilseed crop.

cis-trans Engineering: Advances and Perspectives on Customized Transcriptional Regulation in Plants

Coordinated transcriptional control employing synthetic promoters and transcription factors (TFs) can be used to achieve customized regulation of gene expression in planta. Synthetic promoter technology has yielded a series of promoters with modified cis-regulatory elements that provide useful tools for efficient modulation of gene expression. In addition, the use of zinc fingers (ZFs), transcription activator-like effectors (TALEs), and catalytically inactive clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein 9 (dCas9) has made it feasible to engineer TFs that can produce targeted gene expression regulation; these approaches are particularly effective when artificial TFs are coupled with transcriptional activators or repressors. This review focuses on strategies used to engineer both promoters and TFs in the context of targeted transcriptional regulation. We also discuss the creation of synthetic inducible platforms, which can be used to impart stress tolerance to plants. We propose that combinatorial "cis-trans engineering" using a CRISPR-dCas9-based bipartite module could be used to regulate the expression of multiple target genes. This approach provides an attractive tool for introduction of specific qualitative traits into plants, thus enhancing their overall environmental adaptability.

Combinatorial Transcriptional Control of Plant Specialized Metabolism

Plants produce countless specialized compounds of diverse chemical nature and biological activities. Their biosynthesis often exclusively occurs either in response to environmental stresses or is limited to dedicated anatomical structures. In both scenarios, regulation of biosynthesis appears to be mainly controlled at the transcriptional level, which is generally dependent on a combined interplay of DNA-related mechanisms and the activity of transcription factors that may act in a combinatorial manner. How environmental and developmental cues are integrated into a coordinated cell type-specific stress response has only partially been unraveled so far. Building on the available examples from (metabolic) gene expression, here we propose theoretical models of how this integration of signals may occur at the level of transcriptional control.

TRANSPARENT TESTA GLABRA 1-Dependent Regulation of Flavonoid Biosynthesis

The flavonoid composition of various tissues throughout plant development is of biological relevance and particular interest for breeding. Arabidopsis thaliana TRANSPARENT TESTA GLABRA 1 (AtTTG1) is an essential regulator of late structural genes in flavonoid biosynthesis. Here, we provide a review of the regulation of the pathway's core enzymes through AtTTG1-containing R2R3-MYELOBLASTOSIS-basic HELIX-LOOP-HELIX-WD40 repeat (MBW(AtTTG1)) complexes embedded in an evolutionary context. We present a comprehensive collection of A. thalianattg1 mutants and AtTTG1 orthologs. A plethora of MBW(AtTTG1) mechanisms in regulating the five major TTG1-dependent traits is highlighted.

Six NAC transcription factors involved in response to TYLCV infection in resistant and susceptible tomato cultivars

NAC transcription factors (TFs) belong to plant-specific TFs, which have been identified in many plant species. The NAC TFs act as the nodes of a regulatory network in plant's response to abiotic and biotic stresses. Till now, response of tomato NAC TFs involved in Tomato yellow leaf curl virus (TYLCV) infection is unknown. In the present study, six NAC TFs were identified to respond to TYLCV infection in tomato. We observed that transcripts of four NAC genes (SlNAC20, SlNAC24, SlNAC47, and SlNAC61) were induced after TYLCV infection in resistant tomato cultivar. Virus-induced gene silencing analysis (VIGS) indicated that SlNAC61 played positive roles in response to TYLCV infection. Tomato NAC TFs were not only involved in defense regulation but in development and stress progress. These NAC TFs interacted with other proteins, including protein phosphatase and mitogen-activated protein kinase. Some defense response TFs, such as WRKY, TGA, MYB, NAC, could interact with NAC proteins by binding cis-elements in promoter regions of NAC TFs. These identified tomato NAC TFs cooperated with other TFs and proteins, indicating the complex response mechanism of described NAC TFs involved in TYLCV infection. The results will offer new evidence to further understand the NAC TFs involved in response to TYLCV infection in tomato.

Combinatorial requirement of W- and WT-boxes in microbe-associated molecular pattern-responsive synthetic promoters

The WT-box GGACTTTC belongs to a novel class of MAMP-responsive cis-regulatory sequences that are part of combinatorial elements. Microbe-associated molecular pattern (MAMP)-responsive synthetic promoters were generated with two cis-regulatory modules (CRM1 and CRM2) from the Arabidopsis thaliana WRKY30 promoter. Both modules harbour two W-boxes and one WT-box. Mutation analysis of the synthetic promoters and transient gene expression analysis in parsley protoplasts underline the importance of the W- and WT-boxes for MAMP-responsive gene expression and reveal the combinatorial requirement of at least two boxes for full MAMP responsivity. In the context of the native promoter, CRM1 is required for MAMP responsivity, while CRM2 alone is not sufficient. Yeast one-hybrid screenings using CRM1 with a transcription factor (TF) only prey library select only WRKY factors. Selection of WRKY26, 40, 41, and 70 requires the W-boxes. The WT-box is also required for selection of WRKY26 and 41 in yeast. In plant cells, WRKY26, 40, and 41 act as repressors of MAMP-responsive gene expression, whereas WRKY70 is an activator. To investigate whether the WT-box is also required for WRKY26 and 41 mediated gene expression in plant cells, both were converted into transcriptional activators by adding the GAL4 activating domain (AD). In contrast to yeast, transient gene expression in parsley protoplasts shows that only the W-boxes from CRM1 are required for WRKY41AD-activated reporter gene activity but not the WT-box. In addition, WRKY70-activated reporter gene activity in parsley cells does not require the WT-box of CRM1. The results demonstrate the importance of the WT-box as a new cis-regulatory sequence for MAMP-responsive gene expression. Based on these and earlier results, two types of WT-boxes are proposed.

Functional characterization of tomato membrane-bound NAC transcription factors

Genome-wide analysis was carried out to identify and analyze differential expression pattern of tomato membrane bound NAC transcription factors (SlNACMTFs) during stresses. Two biotic-stress-related SlNACMTFs have been characterized to elucidate their regulatory function. NAC transcription factors are known regulators of stress-related gene expression. As Stresses are perceived and transmitted by membrane-bound proteins, functional characterization of membrane-associated NAC transcription factors in tomato can reveal valuable insight about membrane-mediated stress-signalling. Tomato genome encodes 13 NAC genes which have predicted transmembrane domain(s) (SlNACMTFs). mRNA of 12 SlNACMTFs were readily detected in multiple tissues, and also in polysome isolated from leaf, confirming active transcription and translation from these genes occur under normal physiological condition. Additionally, most of the SlNACMTFs were differentially regulated during stresses and stress-related transcription factor binding sites are prevalent in their promoters. SlNACMTF3 and 8 were majorly regulated in biotic and abiotic stresses. Like other MTFs, SlNACMTF3 was translocated to the plasma membrane, whereas the C-terminus truncated (ΔC) form localized in the cytoplasm and the nucleus. Accordingly, the ΔC forms significantly influenced the activity of promoters harbouring NAC binding sites (NACbs). Furthermore, the NAC domain of these transcription factors could directly interact with an NACbs, and the proteins failed to regulate a promoter lacking a crucial NACbs. Interestingly, the type of influence to an NACbs containing promoter was dependent on the context of the NACbs, as the same SlNACMTF showed an alternative mode of regulation on different promoters, as well as the same promoter activity was oppositely regulated by two different SlNACMTF. Finally, both SlNACMTFs demonstrated the differential regulatory effect on the expression of several stress-related genes by interacting with the putative NACbs in their promoter region, suggesting their direct role in plant stress response.

Deep sequencing and transcriptome analyses to identify genes involved in secoiridoid biosynthesis in the Tibetan medicinal plant Swertia mussotii

Swertia mussotii Franch. is an important traditional Tibetan medicinal plant with pharmacological properties effective in the treatment of various ailments including hepatitis. Secoiridoids are the major bioactive compounds in S. mussotii. To better understand the secoiridoid biosynthesis pathway, we generated transcriptome sequences from the root, leaf, stem, and flower tissues, and performed de novo sequence assembly, yielding 98,613 unique transcripts with an N50 of 1,085 bp. Putative functions could be assigned to 35,029 transcripts (35.52%) based on BLAST searches against annotation databases including GO and KEGG. The expression profiles of 39 candidate transcripts encoding the key enzymes for secoiridoid biosynthesis were examined in different S. mussotii tissues, validated by qRT-PCR, and compared with the homologous genes from S. japonica, a species in the same family, unveiling the gene expression, regulation, and conservation of the pathway. The examination of the accumulated levels of three bioactive compounds, sweroside, swertiamarin, and gentiopicroside, revealed their considerable variations in different tissues, with no significant correlation with the expression profiles of key genes in the pathway, suggesting complex biological behaviours in the coordination of metabolite biosynthesis and accumulation. The genomic dataset and analyses presented here lay the foundation for further research on this important medicinal plant.

The evolution of plant microRNAs: insights from a basal eudicot sacred lotus

microRNAs (miRNAs) are important noncoding small RNAs that regulate mRNAs in eukaryotes. However, under which circumstances different miRNAs/miRNA families exhibit different evolutionary trajectories in plants remains unclear. In this study, we sequenced the small RNAs and degradome from a basal eudicot, sacred lotus (Nelumbo nucifera or lotus), to identify miRNAs and their targets. Combining with public miRNAs, we predicted 57 pre-eudicot miRNA families from different evolutionary stages. We found that miRNA families featuring older age, higher copy and target number tend to show lower propensity for miRNA family loss (PGL) and stronger signature of purifying selection during divergence of temperate and tropical lotus. Further analyses of lotus genome revealed that there is an association between loss of miRNA families in descendent plants and in duplicated genomes. Gene dosage balance is crucial in maintaining those preferentially retained MIRNA duplicates by imposing stronger purifying selection. However, these factors and selection influencing miRNA family evolution are not applicable to the putative MIRNA-likes. Additionally, the MIRNAs participating in lotus pollen-pistil interaction, a conserved process in angiosperms, also have a strong signature of purifying selection. Functionally, sequence divergence in MIRNAs escalates expression divergence of their target genes between temperate and tropical lotus during rhizome and leaf growth. Overall, our study unravels several important factors and selection that determine the miRNA family distribution in plants and duplicated genomes, and provides evidence for functional impact of MIRNA sequence evolution.

Identification of Norway Spruce MYB-bHLH-WDR Transcription Factor Complex Members Linked to Regulation of the Flavonoid Pathway

Transcription factors (TFs) forming MYB-bHLH-WDR complexes are known to regulate the biosynthesis of specialized metabolites in angiosperms through an intricate network. These specialized metabolites participate in a wide range of biological processes including plant growth, development, reproduction as well as in plant immunity. Studying the regulation of their biosynthesis is thus essential. While MYB (TFs) have been previously shown to control specialized metabolism (SM) in gymnosperms, the identity of their partners, in particular bHLH or WDR members, has not yet been revealed. To gain knowledge about MYB-bHLH-WDR transcription factor complexes in gymnosperms and their regulation of SW, we identified two bHLH homologs of AtTT8, six homologs of the MYB transcription factor AtTT2 and one WDR ortholog of AtTTG1 in Norway spruce. We investigated the expression levels of these genes in diverse tissues and upon treatments with various stimuli including methyl-salicylate, methyl-jasmonate, wounding or fungal inoculation. In addition, we also identified protein-protein interactions among different homologs of MYB, bHLH and WDR. Finally, we generated transgenic spruce cell lines overexpressing four of the Norway spruce AtTT2 homologs and observed differential regulation of genes in the flavonoid pathway and flavonoid contents.

A R2R3-MYB transcription factor gene in common wheat (namely TaMYBsm1) involved in enhancement of drought tolerance in transgenic Arabidopsis

We isolated the TaMYBsm1 genes, encoding R2R3-type MYB proteins in common wheat, aimed to uncover the possible molecular mechanisms related to drought response. The TaMYBsm1 genes, TaMYBsm1-A, TaMYBsm1-B and TaMYBsm1-D, were isolated and analyzed from the common wheat cultivar Shimai 15. Their expression patterns under PEG 6000 and mannitol were monitored by semi-quantitative RT-PCR and β-glucuronidase (Gus) assay. The function of TaMYBsm1-D under drought stress in transgenic Arabidopsis plants was investigated, and the germination rate, water loss rate, as well as the proline and malondialdehyde (MDA) content were compared with that in wild type (WT) plants. The expression of three downstream genes (DREB2A, P5CS1 and RD29A) in TaMYBsm1-D transgenic plants was analyzed. The R2R3-MYB domains of the MYBsm1 proteins were highly conserved in plants. In addition, the TaMYBsm1 proteins were targeted to the nucleus and contained transcriptional activation domains (TADs). Gus assay and semi-quantitative RT-PCR analysis demonstrated that the TaMYBsm1 genes were up-regulated when the wheat was treated by PEG and mannitol. Compared with WT plants, the germination rates were much higher, but the water loss rates were much lower in TaMYBsm1-D overexpression plants. TaMYBsm1-D transgenic plants showed distinct higher proline contents but a lower MDA content than the WT plants. The three downstream genes were highly expressed in TaMYBsm1-D transgenic plants. We concluded from these results that TaMYBsm1 genes play an important role in plant drought stress tolerance through up-regulation of DREB2A, P5CS1 and RD29A. The increase of proline content and decrease of MDA content may also be involved in the drought response.

The "putative" role of transcription factors from HlWRKY family in the regulation of the final steps of prenylflavonid and bitter acids biosynthesis in hop (Humulus lupulus L.)

Lupulin glands localized in female hop (Humulus lupulus L.) cones are valuable source of bitter acids, essential oils and polyphenols. These compounds are used in brewing industry and are important for biomedical applications. In this study we describe the potential effect of transcription factors from WRKY family in the activation of the final steps of lupulin biosynthesis. In particular, lupulin gland-specific transcription factor HlWRKY1 that shows significant similarity to AtWRKY75, has ability to activate the set of promoters driving key genes of xanthohumol and bitter acids biosynthesis such as chalcone synthase H1, valerophenone synthase, prenyltransferase 1, 1L and 2 and O-methyltransferase-1. When combined with co-factor HlWDR1 and silencing suppressor p19, HlWRKY1 is able to enhance transient expression of gus gene driven by Omt1 and Chs_H1 promoters to significant level as compared to 35S promoter of CaMV in Nicotiana. benthamiana. Transformation of hop with dual Agrobacterium vector bearing HlWRKY1/HlWDR1 led to ectopic overexpression of these transgenes and further activation of lupulin-specific genes expression in hop leaves. It was further showed that (1) HlWRKY1 is endowed with promoter autoactivation; (2) It is regulated by post-transcriptional gene silencing (PTGS) mechanism; (3) It is stimulated by kinase co-expression. Since HlWRKY1 promotes expression of lupulin-specific HlMyb3 gene therefore it can constitute a significant component in hop lupulin regulation network. Putative involvement of HlWRKY1 in the regulation of lupulin biosynthesis may suggest the original physiological function of lupulin components in hop as flower and seed protective compounds.

Diversity, expansion, and evolutionary novelty of plant DNA-binding transcription factor families

Plant transcription factors (TFs) that interact with specific sequences via DNA-binding domains are crucial for regulating transcriptional initiation and are fundamental to plant development and environmental response. In addition, expansion of TF families has allowed functional divergence of duplicate copies, which has contributed to novel, and in some cases adaptive, traits in plants. Thus, TFs are central to the generation of the diverse plant species that we see today. Major plant agronomic traits, including those relevant to domestication, have also frequently arisen through changes in TF coding sequence or expression patterns. Here our goal is to provide an overview of plant TF evolution by first comparing the diversity of DNA-binding domains and the sizes of these domain families in plants and other eukaryotes. Because TFs are among the most highly expanded gene families in plants, the birth and death process of TFs as well as the mechanisms contributing to their retention are discussed. We also provide recent examples of how TFs have contributed to novel traits that are important in plant evolution and in agriculture.This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

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.

[The roles of MYB transcription factors on plant defense responses and its molecular mechanism.]

Transcriptional regulation of defense gene expression is a crucial part of plant defense responses in plant defense environment stresses. As one of the largest plant transcription factor families, MYB (v-myb avian myeloblastosis viral on-cogene homolog) transcription factors play an important role in plant stress tolerance. In this paper, we review the structural features, functional characterization and molecular mechanism of MYB transcription factor family, and discuss the regula-tory roles of transcription factors in plant defense responses.

Characterization of NtREL1, a novel root-specific gene from tobacco, and upstream promoter activity analysis in homologous and heterologous hosts

KEY MESSAGE

A novel root-specific gene and its upstream promoter were cloned and characterized for potential application in root-specific expression of transgenes. The root is an important plant organ subjected to many biotic and abiotic stresses, such as infection by Ralstonia solanacearum. To isolate tobacco root-specific promoters for genetic applications, microarray screening was performed to identify genes highly and specifically expressed in the root. One root-specific gene encoding an extensin-like protein (NtREL1) was isolated, and its expression pattern was further characterized by both microarray analysis and reverse transcription-polymerase chain reaction. NtREL1 was highly expressed only in roots but not in any other organ. NtREL1 expression was affected by hormone treatment (salicylic acid, abscisic acid, and ethephon) as well as low temperature, drought, and R. solanacearum infection. A full-length 849 bp cDNA containing a 657-nucleotide open reading frame was cloned by Rapid Amplification of cDNA Ends. Subsequently, a fragment of 1,574 bp upstream of NtREL1 was isolated by flanking PCR and named pNtREL1. This promoter fragment contains TATA, GATA, and CAAT-boxes, the basic elements of a promoter, and six root-specific expression elements, namely OSE1, OSE2, ROOTMOTIFTAPOX1, SURECOREATSULTR11, P1BS, and WUSATAg. A construct containing the bacterial uidA reporter gene (β-glucuronidase, GUS) driven by the pNtREL1 promoter was transformed into tobacco plants. GUS staining was only detected in the root, but not in leaves and stems. Additionally, transgenic tobacco plants containing peanut resveratrol synthase gene (AhRS) driven by the pNtREL1 promoter produced resveratrol only in the root. Thus, the pNtREL1 promoter can be used to direct root-specific expression of target genes to protect the root from stress or for biological studies.

The VviMYB80 Gene is Abnormally Expressed in Vitis vinifera L. cv. 'Zhong Shan Hong' and its Expression in Tobacco Driven by the 35S Promoter Causes Male Sterility

Anther development is a very precise and complicated process. In Arabidopsis, the AtMYB80 transcription factor regulates genes involved in pollen development and controls the timing of tapetal programmed cell death (PCD). In this study, we isolated and characterized cDNA for VviMYB80 expressed in flower buds of male-sterile Vitis vinifera L. cv. 'Zhong Shan Hong', a late-maturing cultivar derived from self-progeny of cv. 'Wink'. VviMYB80 belongs to the MYB80 subfamily and clusters with AtMYB35/TDF1 in a distinct clade. We found that in flower buds, expression of the VviMYB80 gene in cv. 'Zhong Shan Hong' sharply increased at the tetrad stage, resulting in a higher and earlier transcript level than that found in cv. 'Wink'. Expression of the VviMYB80 gene, driven by the 35S promoter, caused pleiotropic effects on the stamens, including smaller and shriveled anthers, delayed dehiscence, fewer seeds, shorter anther filaments, distorted pollen shape and a lack of cytoplasm, with the tapetum exhibiting hypertrophy in transformed tobacco. These results suggest that VviMYB80 may play an important role in stamen development and that expression of VviMYB80 driven by the 35S promoter in tobacco induces male sterility.

Analysis of Microbe-Associated Molecular Pattern-Responsive Synthetic Promoters with the Parsley Protoplast System

Plants recognize pathogens by microbe-associated molecular patterns (MAMPs) and subsequently induce an immune response. The regulation of gene expression during the immune response depends largely on cis-sequences conserved in promoters of MAMP-responsive genes. These cis-sequences can be analyzed by constructing synthetic promoters linked to a reporter gene and by testing these constructs in transient expression systems. Here, the use of the parsley (Petroselinum crispum) protoplast system for analyzing MAMP-responsive synthetic promoters is described. The synthetic promoter consists of four copies of a potential MAMP-responsive cis-sequence cloned upstream of a minimal promoter and the uidA reporter gene. The reporter plasmid contains a second reporter gene, which is constitutively expressed and hence eliminates the requirement of a second plasmid used as a transformation control. The reporter plasmid is transformed into parsley protoplasts that are elicited by the MAMP Pep25. The MAMP responsiveness is validated by comparing the reporter gene activity from MAMP-treated and untreated cells and by normalizing reporter gene activity using the constitutively expressed reporter gene.

Transcriptional regulatory networks in Arabidopsis thaliana during single and combined stresses

Differentially evolved responses to various stress conditions in plants are controlled by complex regulatory circuits of transcriptional activators, and repressors, such as transcription factors (TFs). To understand the general and condition-specific activities of the TFs and their regulatory relationships with the target genes (TGs), we have used a homogeneous stress gene expression dataset generated on ten natural ecotypes of the model plant Arabidopsis thaliana, during five single and six combined stress conditions. Knowledge-based profiles of binding sites for 25 stress-responsive TF families (187 TFs) were generated and tested for their enrichment in the regulatory regions of the associated TGs. Condition-dependent regulatory sub-networks have shed light on the differential utilization of the underlying network topology, by stress-specific regulators and multifunctional regulators. The multifunctional regulators maintain the core stress response processes while the transient regulators confer the specificity to certain conditions. Clustering patterns of transcription factor binding sites (TFBS) have reflected the combinatorial nature of transcriptional regulation, and suggested the putative role of the homotypic clusters of TFBS towards maintaining transcriptional robustness against cis-regulatory mutations to facilitate the preservation of stress response processes. The Gene Ontology enrichment analysis of the TGs reflected sequential regulation of stress response mechanisms in plants.

Operative photo assimilation associated proteome modulations are critical for iron-dependent cadmium tolerance in Oryza sativa L

Iron-dependent Cd tolerance in Oryza sativa L. cv 7029 had been explored. Photo assimilatory process such as photosynthesis and nitrogen fixation found to be functional in the presence of excess Fe during Cd stress. Cd-inducible Fe deficiency demonstrated with upregulation of iron uptake gene families such as OsIRT and OsYSL was attributed as foremost reason for retardation of photo assimilation in the course of Cd treatment. Upholding of photo assimilation during Fe supplement was associated with proteome modulations. Monitoring of proteome responses in leaf and root revealed proteins imperative for alleviation of Cd stress. Specifically, proteins that take part in photosynthesis, glucose metabolism, nitrogen fixation, and abiotic stress tolerance played key role in Fe-dependent Cd detoxification. It is concluded that Fe supply help to operate photo assimilation which enable favorable proteome responses of Cd stress tolerance.

Deciphering Cis-Regulatory Element Mediated Combinatorial Regulation in Rice under Blast Infected Condition

Combinations of cis-regulatory elements (CREs) present at the promoters facilitate the binding of several transcription factors (TFs), thereby altering the consequent gene expressions. Due to the eminent complexity of the regulatory mechanism, the combinatorics of CRE-mediated transcriptional regulation has been elusive. In this work, we have developed a new methodology that quantifies the co-occurrence tendencies of CREs present in a set of promoter sequences; these co-occurrence scores are filtered in three consecutive steps to test their statistical significance; and the significantly co-occurring CRE pairs are presented as networks. These networks of co-occurring CREs are further transformed to derive higher order of regulatory combinatorics. We have further applied this methodology on the differentially up-regulated gene-sets of rice tissues under fungal (Magnaporthe) infected conditions to demonstrate how it helps to understand the CRE-mediated combinatorial gene regulation. Our analysis includes a wide spectrum of biologically important results. The CRE pairs having a strong tendency to co-occur often exhibit very similar joint distribution patterns at the promoters of rice. We couple the network approach with experimental results of plant gene regulation and defense mechanisms and find evidences of auto and cross regulation among TF families, cross-talk among multiple hormone signaling pathways, similarities and dissimilarities in regulatory combinatorics between different tissues, etc. Our analyses have pointed a highly distributed nature of the combinatorial gene regulation facilitating an efficient alteration in response to fungal attack. All together, our proposed methodology could be an important approach in understanding the combinatorial gene regulation. It can be further applied to unravel the tissue and/or condition specific combinatorial gene regulation in other eukaryotic systems with the availability of annotated genomic sequences and suitable experimental data.

Clustering and Differential Alignment Algorithm: Identification of Early Stage Regulators in the Arabidopsis thaliana Iron Deficiency Response

Time course transcriptome datasets are commonly used to predict key gene regulators associated with stress responses and to explore gene functionality. Techniques developed to extract causal relationships between genes from high throughput time course expression data are limited by low signal levels coupled with noise and sparseness in time points. We deal with these limitations by proposing the Cluster and Differential Alignment Algorithm (CDAA). This algorithm was designed to process transcriptome data by first grouping genes based on stages of activity and then using similarities in gene expression to predict influential connections between individual genes. Regulatory relationships are assigned based on pairwise alignment scores generated using the expression patterns of two genes and some inferred delay between the regulator and the observed activity of the target. We applied the CDAA to an iron deficiency time course microarray dataset to identify regulators that influence 7 target transcription factors known to participate in the Arabidopsis thaliana iron deficiency response. The algorithm predicted that 7 regulators previously unlinked to iron homeostasis influence the expression of these known transcription factors. We validated over half of predicted influential relationships using qRT-PCR expression analysis in mutant backgrounds. One predicted regulator-target relationship was shown to be a direct binding interaction according to yeast one-hybrid (Y1H) analysis. These results serve as a proof of concept emphasizing the utility of the CDAA for identifying unknown or missing nodes in regulatory cascades, providing the fundamental knowledge needed for constructing predictive gene regulatory networks. We propose that this tool can be used successfully for similar time course datasets to extract additional information and infer reliable regulatory connections for individual genes.

A combination of gene expression ranking and co-expression network analysis increases discovery rate in large-scale mutant screens for novel Arabidopsis thaliana abiotic stress genes

As challenges to food security increase, the demand for lead genes for improving crop production is growing. However, genetic screens of plant mutants typically yield very low frequencies of desired phenotypes. Here, we present a powerful computational approach for selecting candidate genes for screening insertion mutants. We combined ranking of Arabidopsis thaliana regulatory genes according to their expression in response to multiple abiotic stresses (Multiple Stress [MST] score), with stress-responsive RNA co-expression network analysis to select candidate multiple stress regulatory (MSTR) genes. Screening of 62 T-DNA insertion mutants defective in candidate MSTR genes, for abiotic stress germination phenotypes yielded a remarkable hit rate of up to 62%; this gene discovery rate is 48-fold greater than that of other large-scale insertional mutant screens. Moreover, the MST score of these genes could be used to prioritize them for screening. To evaluate the contribution of the co-expression analysis, we screened 64 additional mutant lines of MST-scored genes that did not appear in the RNA co-expression network. The screening of these MST-scored genes yielded a gene discovery rate of 36%, which is much higher than that of classic mutant screens but not as high as when picking candidate genes from the co-expression network. The MSTR co-expression network that we created, AraSTressRegNet is publicly available at http://netbio.bgu.ac.il/arnet. This systems biology-based screening approach combining gene ranking and network analysis could be generally applicable to enhancing identification of genes regulating additional processes in plants and other organisms provided that suitable transcriptome data are available.

WsMAGO2, a duplicated MAGO NASHI protein with fertility attributes interacts with MPF2-like MADS-box proteins

WsMAGO2 a duplicated protein in Withania through interactions with MPF2-like proteins affects male fertility by producing fewer flowers and aborted non-viable pollens/seeds regulated by anther-specific GAATTTGTGA motif. The MAGO NASHIs are highly conserved genes that encode proteins known to be involved in RNA physiology and many other developmental processes including germ cell differentiation in animals. However, their structural and functional implications in plants as fertility function proteins remained fragmented. MAGO (shorter name of MAGO NASHI) proteins form heterodimers with MPF2-like MADS-box proteins which are recruited in calyx identity and male fertility in Solanaceous plants. Four MAGO genes namely WsMAGO1 and WsMAGO2 and TaMAGO1 and TaMAGO2 were isolated from Withania somnifera and Tubocapsicum anomalum, respectively. These genes have duplicated probably due to whole genome duplication event. Dysfunction of WsMAGO2 through double-stranded RNAi in Withania revealed suppression of RNA transcripts, non-viable pollens, fewer flowers and aborted non-viable seeds in the developing berry suggesting a role of this protein in many traits particularly male fertility. WsMAGO2 flaunted stronger yeast 2-hybrid interactions with MPF2-like proteins WSA206, WSB206 and TAB201 than other MAGO counterparts. The native transcripts of WsMAGO2 culminated in stamens and seed-bearing berries though other MAGO orthologs also exhibited expression albeit at lower level. Coding sequences of the two orthologs are highly conserved, but they differ substantially in their upstream promoter regions. Remarkably, WsMAGO2 promoter is enriched with many anther-specific cis-motifs common in fertility function genes promoters. Among them, disruption of GAATTTGTGA abolished YFP/GUS gene expression in anthers alluding towards its involvement in regulating expression of MAGO in anther. Our findings support a possible recruitment of WsMAGO2 in fertility trait in Withania. These genes have practical application in hybrid production through cytoplasmic male sterility maintenance for enhancement in crops yield.

Combined correlation-based network and mQTL analyses efficiently identified loci for branched-chain amino acid, serine to threonine, and proline metabolism in tomato seeds

Correlation-based network analysis (CNA) of the metabolic profiles of seeds of a tomato introgression line mapping population revealed a clique of proteinogenic amino acids: Gly, Ile, Pro, Ser, Thr, and Val. Correlations between profiles of these amino acids exhibited a statistically significant average correlation coefficient of 0.84 as compared with an average correlation coefficient of 0.39 over the 16 119 other metabolite cliques containing six metabolites. In silico removal of cliques was used to quantify their importance in determining seminal network properties, highlighting the strong effects of the amino acid clique. Quantitative trait locus analysis revealed co-localization for the six amino acids on chromosome 2, 4 and 10. Sequence analysis identified a unique set of 10 genes on chromosome 2 only, which were associated with amino acid metabolism and specifically the metabolism of Ser-Gly and their conversion into branched-chain amino acids. Metabolite profiling of a set of sublines, with introgressions on chromosome 2, identified a significant change in the abundance of the six amino acids in comparison with M82. Expression analysis of candidate genes affecting Ser metabolism matched the observation from the metabolite data, suggesting a coordinated behavior of the level of these amino acids at the genetic level. Analysis of transcription factor binding sites in the promoter regions of the identified genes suggested combinatorial response to light and the circadian clock.

Expression patterns and protein structure of a lipid transfer protein END1 from Arabidopsis

Arabidopsis END1-LIKE (AtEND1) was identified as a homolog of the barley endosperm-specific gene END1 and provides a model for the study of this class of genes and their products. The END1 is expressed in the endosperm transfer cells (ETC) of grasses. The ETC are responsible for transfer of nutrients from maternal tissues to the developing endosperm. Identification of several ETC-specific genes encoding lipid transfer proteins (LTP), including the END1, provided excellent markers for identification of ETC during seed development. To understand how AtEND1 forms complexes with lipid molecules, a three-dimensional (3D) molecular model was generated and reconciled with AtEND1 function. The spatial and temporal expression patterns of AtEND1 were examined in transgenic Arabidopsis plants transformed with an AtEND1 promoter-GUS fusion construct. The AtEND1 promoter was found to be seed and pollen specific. In contrast to ETC-specific expression of homologous genes in wheat and barley, expression of AtEND1 is less specific. It was observed in ovules and a few gametophytic tissues. A series of AtEND1 promoter deletions fused to coding sequence (CDS) of the uidA were transformed in Arabidopsis and the promoter region responsible for AtEND1 expression was identified. A 163 bp fragment of the promoter was found to be sufficient for both spatial and temporal patterns of expression reflecting that of AtEND1. Our data suggest that AtEND1 could be used as a marker gene for gametophytic tissues and developing endosperm. The role of the gene is unclear but it may be involved in fertilization and/or endosperm cellularization.

Transcription factors and glyoxylate cycle genes prominent in the transition of soybean cotyledons to the first functional leaves of the seedling

During early seedling growth, the cotyledons transition from a storage tissue to become effectively the first leaf-like structures of the plant. In this programmed developmental process, they likely undergo a massive change in gene expression to redirect their metabolism and physiological processes. To define the developmental shifts in gene expression and begin to understand the gene regulatory networks that set this transition in motion, we carried out high-throughput RNA sequencing of cotyledons from seven developmental stages of soybean seedlings. We identified 154 gene models with high expression exclusively in the early seedling stages. A significant number (about 25 %) of those genes with known annotations were involved in carbohydrate metabolism. A detailed examination of glyoxylate cycle genes revealed the upregulation of their expression in the early stages of development. A total of approximately 50 % of the highly expressed genes whose expression peaked in the mid-developmental stages encoded ribosomal family proteins. Our analysis also identified 219 gene models with high expression at late developmental stages. The majority of these genes are involved in photosynthesis, including photosystem I- and II-associated genes. Additionally, the advantage of RNA-Seq to detect genes expressed at low levels revealed approximately 460 transcription factors with notable expression in at least one stage of the developing soybean seedling. Relatively over-represented transcription factor genes encode AP2, zinc finger, NAC, WRKY, and MYB families. These transcription factor genes may lead to the transcriptional reprogramming during the transition of seedling cotyledons from storage tissue to metabolically active organs that serve as the first functional leaves of the plant.

Transcriptional responses of Arabidopsis thaliana to chewing and sucking insect herbivores

We tested the hypothesis that Arabidopsis can recognize and respond differentially to insect species at the transcriptional level using a genome wide microarray. Transcriptional reprogramming was characterized using co-expression analysis in damaged and undamaged leaves at two times in response to mechanical wounding and four insect species. In all, 2778 (10.6%) of annotated genes on the array were differentially expressed in at least one treatment. Responses differed mainly between aphid and caterpillar and sampling times. Responses to aphids and caterpillars shared only 10% of up-regulated and 8% of down-regulated genes. Responses to two caterpillars shared 21 and 12% of up- and down-regulated genes, whereas responses to the two aphids shared only 7 and 4% of up-regulated and down-regulated genes. Overlap in genes expressed between 6 and 24 h was 3-15%, and depended on the insect species. Responses in attacked and unattacked leaves differed at 6 h but converged by 24 h. Genes responding to the insects are also responsive to many stressors and included primary metabolism. Aphids down-regulated amino acid catabolism; caterpillars stimulated production of amino acids involved in glucosinolate synthesis. Co-expression analysis revealed 17 response networks. Transcription factors were a major portion of differentially expressed genes throughout and responsive genes shared most of the known or postulated binding sites. However, cis-element composition of genes down regulated by the aphid M. persicae was unique, as were those of genes down-regulated by caterpillars. As many as 20 cis-elements were over-represented in one or more treatments, including some from well-characterized classes and others as yet uncharacterized. We suggest that transcriptional changes elicited by wounding and insects are heavily influenced by transcription factors and involve both enrichment of a common set of cis-elements and a unique enrichment of a few cis-elements in responding genes.

A functional and evolutionary perspective on transcription factor binding in Arabidopsis thaliana

Understanding the mechanisms underlying gene regulation is paramount to comprehend the translation from genotype to phenotype. The two are connected by gene expression, and it is generally thought that variation in transcription factor (TF) function is an important determinant of phenotypic evolution. We analyzed publicly available genome-wide chromatin immunoprecipitation experiments for 27 TFs in Arabidopsis thaliana and constructed an experimental network containing 46,619 regulatory interactions and 15,188 target genes. We identified hub targets and highly occupied target (HOT) regions, which are enriched for genes involved in development, stimulus responses, signaling, and gene regulatory processes in the currently profiled network. We provide several lines of evidence that TF binding at plant HOT regions is functional, in contrast to that in animals, and not merely the result of accessible chromatin. HOT regions harbor specific DNA motifs, are enriched for differentially expressed genes, and are often conserved across crucifers and dicots, even though they are not under higher levels of purifying selection than non-HOT regions. Distal bound regions are under purifying selection as well and are enriched for a chromatin state showing regulation by the Polycomb repressive complex. Gene expression complexity is positively correlated with the total number of bound TFs, revealing insights in the regulatory code for genes with different expression breadths. The integration of noncanonical and canonical DNA motif information yields new hypotheses on cobinding and tethering between specific TFs involved in flowering and light regulation.

Temporal and spatial control of gene expression in horticultural crops

Biotechnology provides plant breeders an additional tool to improve various traits desired by growers and consumers of horticultural crops. It also provides genetic solutions to major problems affecting horticultural crops and can be a means for rapid improvement of a cultivar. With the availability of a number of horticultural genome sequences, it has become relatively easier to utilize these resources to identify DNA sequences for both basic and applied research. Promoters play a key role in plant gene expression and the regulation of gene expression. In recent years, rapid progress has been made on the isolation and evaluation of plant-derived promoters and their use in horticultural crops, as more and more species become amenable to genetic transformation. Our understanding of the tools and techniques of horticultural plant biotechnology has now evolved from a discovery phase to an implementation phase. The availability of a large number of promoters derived from horticultural plants opens up the field for utilization of native sequences and improving crops using precision breeding. In this review, we look at the temporal and spatial control of gene expression in horticultural crops and the usage of a variety of promoters either isolated from horticultural crops or used in horticultural crop improvement.

Nuclear factor Y-mediated H3K27me3 demethylation of the SOC1 locus orchestrates flowering responses of Arabidopsis

Nuclear factor Y (NF-Y) is a conserved heterotrimeric transcription factor complex that binds to the CCAAT motifs within the promoter region of many genes. In plants, a large number of genes code for variants of each NF-YA, B or C subunit that can assemble in a combinatorial fashion. Here, we report the discovery of an Arabidopsis NF-Y complex that exerts epigenetic control over flowering time by integrating environmental and developmental signals. We show that NF-Y interacts with CONSTANS in the photoperiod pathway and DELLAs in the gibberellin pathway, to directly regulate the transcription of SOC1, a major floral pathway integrator. This NF-Y complex binds to a unique cis-element within the SOC1 promoter to modulate trimethylated H3K27 levels, partly through a H3K27 demethylase REF6. Our findings establish NF-Y complexes as critical mediators of epigenetic marks that regulate the response to environmental or intrinsic signals in plants.