Identification of plant promoter constituents by analysis of local distribution of short sequences

Core promoterome of barley embryo

Precise localization and dissection of gene promoters are key to understanding transcriptional gene regulation and to successful bioengineering applications. The core RNA polymerase II initiation machinery is highly conserved among eukaryotes, leading to a general expectation of equivalent underlying mechanisms. Still, less is known about promoters in the plant kingdom. In this study, we employed cap analysis of gene expression (CAGE) at three embryonic developmental stages in barley to accurately map, annotate, and quantify transcription initiation events. Unsupervised discovery of de novo sequence clusters grouped promoters based on characteristic initiator and position-specific core-promoter motifs. This grouping was complemented by the annotation of transcription factor binding site (TFBS) motifs. Integration with genome-wide epigenomic data sets and gene ontology (GO) enrichment analysis further delineated the chromatin environments and functional roles of genes associated with distinct promoter categories. The TATA-box presence governs all features explored, supporting the general model of two separate genomic regulatory environments. We describe the extent and implications of alternative transcription initiation events, including those that are specific to developmental stages, which can affect the protein sequence or the presence of regions that regulate translation. The generated promoterome dataset provides a valuable genomic resource for enhancing the functional annotation of the barley genome. It also offers insights into the transcriptional regulation of individual genes and presents opportunities for the informed manipulation of promoter architecture, with the aim of enhancing traits of agronomic importance.

The landscape of transcription factor promoter activity during vegetative development in Marchantia

Transcription factors (TFs) are essential for the regulation of gene expression and cell fate determination. Characterizing the transcriptional activity of TF genes in space and time is a critical step toward understanding complex biological systems. The vegetative gametophyte meristems of bryophytes share some characteristics with the shoot apical meristems of flowering plants. However, the identity and expression profiles of TFs associated with gametophyte organization are largely unknown. With only ∼450 putative TF genes, Marchantia (Marchantia polymorpha) is an outstanding model system for plant systems biology. We have generated a near-complete collection of promoter elements derived from Marchantia TF genes. We experimentally tested reporter fusions for all the TF promoters in the collection and systematically analyzed expression patterns in Marchantia gemmae. This allowed us to build a map of expression domains in early vegetative development and identify a set of TF-derived promoters that are active in the stem-cell zone. The cell markers provide additional tools and insight into the dynamic regulation of the gametophytic meristem and its evolution. In addition, we provide an online database of expression patterns for all promoters in the collection. We expect that these promoter elements will be useful for cell-type-specific expression, synthetic biology applications, and functional genomics.

Predmoter-cross-species prediction of plant promoter and enhancer regions

Motivation

Identifying cis-regulatory elements (CREs) is crucial for analyzing gene regulatory networks. Next generation sequencing methods were developed to identify CREs but represent a considerable expenditure for targeted analysis of few genomic loci. Thus, predicting the outputs of these methods would significantly cut costs and time investment.

Results

We present Predmoter, a deep neural network that predicts base-wise Assay for Transposase Accessible Chromatin using sequencing (ATAC-seq) and histone Chromatin immunoprecipitation DNA-sequencing (ChIP-seq) read coverage for plant genomes. Predmoter uses only the DNA sequence as input. We trained our final model on 21 species for 13 of which ATAC-seq data and for 17 of which ChIP-seq data was publicly available. We evaluated our models on Arabidopsis thaliana and Oryza sativa. Our best models showed accurate predictions in peak position and pattern for ATAC- and histone ChIP-seq. Annotating putatively accessible chromatin regions provides valuable input for the identification of CREs. In conjunction with other in silico data, this can significantly reduce the search space for experimentally verifiable DNA-protein interaction pairs.

Availability and implementation

The source code for Predmoter is available at: https://github.com/weberlab-hhu/Predmoter. Predmoter takes a fasta file as input and outputs h5, and optionally bigWig and bedGraph files.

Calcium regulation of the Arabidopsis Na+/K+ transporter HKT1;1 improves seed germination under salt stress

Calcium is known to improve seed-germination rates under salt stress. We investigated the involvement of calcium ions (Ca2+) in regulating HIGH-AFFINITY K+ TRANSPORTER 1 (HKT1; 1), which encodes a Na+/K+ transporter, and its post-translational regulator TYPE 2C PROTEIN PHOSPHATASE 49 (PP2C49), in germinating Arabidopsis (Arabidopsis thaliana) seedlings. Germination rates of hkt1 mutant seeds under salt stress remained unchanged by CaCl2 treatment in wild-type Arabidopsis, whereas pp2c49 mutant seeds displayed improved salt-stress tolerance in the absence of CaCl2 supplementation. Analysis of HKT1;1 and PP2C49 promoter activity revealed that CaCl2 treatment results in radicle-focused expression of HKT1;1 and reduction of the native radicle-exclusive expression of PP2C49. Ion-content analysis indicated that CaCl2 treatment improves K+ retention in germinating wild-type seedlings under salt stress, but not in hkt1 seedlings. Transgenic seedlings designed to exclusively express HKT1;1 in the radicle during germination displayed higher germination rates under salt stress than the wild type in the absence of CaCl2 treatment. Transcriptome analysis of germinating seedlings treated with CaCl2, NaCl, or both revealed 118 upregulated and 94 downregulated genes as responsive to the combined treatment. Bioinformatics analysis of the upstream sequences of CaCl2-NaCl-treatment-responsive upregulated genes revealed the abscisic acid response element CACGTGTC, a potential CaM-binding transcription activator-binding motif, as most prominent. Our findings suggest a key role for Ca2+ in mediating salt-stress responses during germination by regulating genes that function to maintain Na+ and K+ homeostasis, which is vital for seed germination under salt stress.

An improved method for the highly specific detection of transcription start sites

Precise detection of the transcriptional start site (TSS) is a key for characterizing transcriptional regulation of genes and for annotation of newly sequenced genomes. Here, we describe the development of an improved method, designated 'TSS-seq2.' This method is an iterative improvement of TSS-seq, a previously published enzymatic cap-structure conversion method to detect TSSs in base sequences. By modifying the original procedure, including by introducing split ligation at the key cap-selection step, the yield and the accuracy of the reaction has been substantially improved. For example, TSS-seq2 can be conducted using as little as 5 ng of total RNA with an overall accuracy of 96%; this yield a less-biased and more precise detection of TSS. We then applied TSS-seq2 for TSS analysis of four plant species that had not yet been analyzed by any previous TSS method.

A comparative analysis of stably expressed genes across diverse angiosperms exposes flexibility in underlying promoter architecture

Promoters regulate both the amplitude and pattern of gene expression-key factors needed for optimization of many synthetic biology applications. Previous work in Arabidopsis found that promoters that contain a TATA-box element tend to be expressed only under specific conditions or in particular tissues, while promoters that lack any known promoter elements, thus designated as Coreless, tend to be expressed more uniformly. To test whether this trend represents a conserved promoter design rule, we identified stably expressed genes across multiple angiosperm species using publicly available RNA-seq data. Comparisons between core promoter architectures and gene expression stability revealed differences in core promoter usage in monocots and eudicots. Furthermore, when tracing the evolution of a given promoter across species, we found that core promoter type was not a strong predictor of expression pattern. Our analysis suggests that core promoter types are correlative rather than causative in promoter expression patterns and highlights the challenges in finding or building constitutive promoters that will work across diverse plant species.

Plant Promoters and Terminators for High-Precision Bioengineering

High-precision bioengineering and synthetic biology require fine-tuning gene expression at both transcriptional and posttranscriptional levels. Gene transcription is tightly regulated by promoters and terminators. Promoters determine the timing, tissues and cells, and levels of the expression of genes. Terminators mediate transcription termination of genes and affect mRNA levels posttranscriptionally, e.g., the 3'-end processing, stability, translation efficiency, and nuclear to cytoplasmic export of mRNAs. The promoter and terminator combination affects gene expression. In the present article, we review the function and features of plant core promoters, proximal and distal promoters, and terminators, and their effects on and benchmarking strategies for regulating gene expression.

A comparative analysis of stably expressed genes across diverse angiosperms exposes flexibility in underlying promoter architecture

Promoters regulate both the amplitude and pattern of gene expression-key factors needed for optimization of many synthetic biology applications. Previous work in Arabidopsis found that promoters that contain a TATA-box element tend to be expressed only under specific conditions or in particular tissues, while promoters which lack any known promoter elements, thus designated as Coreless, tend to be expressed more ubiquitously. To test whether this trend represents a conserved promoter design rule, we identified stably expressed genes across multiple angiosperm species using publicly available RNA-seq data. Comparisons between core promoter architectures and gene expression stability revealed differences in core promoter usage in monocots and eudicots. Furthermore, when tracing the evolution of a given promoter across species, we found that core promoter type was not a strong predictor of expression stability. Our analysis suggests that core promoter types are correlative rather than causative in promoter expression patterns and highlights the challenges in finding or building constitutive promoters that will work across diverse plant species.

On the Role of TATA Boxes and TATA-Binding Protein in Arabidopsis thaliana

For transcription initiation by RNA polymerase II (Pol II), all eukaryotes require assembly of basal transcription machinery on the core promoter, a region located approximately in the locus spanning a transcription start site (-50; +50 bp). Although Pol II is a complex multi-subunit enzyme conserved among all eukaryotes, it cannot initiate transcription without the participation of many other proteins. Transcription initiation on TATA-containing promoters requires the assembly of the preinitiation complex; this process is triggered by an interaction of TATA-binding protein (TBP, a component of the general transcription factor TFIID (transcription factor II D)) with a TATA box. The interaction of TBP with various TATA boxes in plants, in particular Arabidopsis thaliana, has hardly been investigated, except for a few early studies that addressed the role of a TATA box and substitutions in it in plant transcription systems. This is despite the fact that the interaction of TBP with TATA boxes and their variants can be used to regulate transcription. In this review, we examine the roles of some general transcription factors in the assembly of the basal transcription complex, as well as functions of TATA boxes of the model plant A. thaliana. We review examples showing not only the involvement of TATA boxes in the initiation of transcription machinery assembly but also their indirect participation in plant adaptation to environmental conditions in responses to light and other phenomena. Examples of an influence of the expression levels of A. thaliana TBP1 and TBP2 on morphological traits of the plants are also examined. We summarize available functional data on these two early players that trigger the assembly of transcription machinery. This information will deepen the understanding of the mechanisms underlying transcription by Pol II in plants and will help to utilize the functions of the interaction of TBP with TATA boxes in practice.

Transcriptional and post-transcriptional controls for tuning gene expression in plants

Plant biotechnologists seek to modify plants through genetic reprogramming, but our ability to precisely control gene expression in plants is still limited. Here, we review transcription and translation in the model plants Arabidopsis thaliana and Nicotiana benthamiana with an eye toward control points that may be used to predictably modify gene expression. We highlight differences in gene expression requirements between these plants and other species, and discuss the ways in which our understanding of gene expression has been used to engineer plants. This review is intended to serve as a resource for plant scientists looking to achieve precise control over gene expression.

Whole-genome sequencing of a year-round fruiting jackfruit (Artocarpus heterophyllus Lam.) reveals high levels of single nucleotide variation

Jackfruit (Artocarpus heterophyllus Lam.) is the national fruit of Bangladesh and produces fruit in the summer season only. However, jackfruit is not commercially grown in Bangladesh because of an extremely high variation in fruit quality, short seasonal fruiting (June-August) and susceptibility to abiotic stresses. Conversely, a year-round high yielding (ca. 4-fold higher than the seasonal variety) jackfruit variety, BARI Kanthal-3 developed by the Bangladesh Agricultural Research Institute (BARI) derived from a wild accession found in Ramgarh of Chattogram Hiltracts of Bangladesh, provides fruits from September to June. This study aimed to generate a draft whole-genome sequence (WGS) of BARI Kanthal-3 to obtain molecular insights including genes associated with year-round fruiting trait of this important unique variety. The estimated genome size of BARI Kanthal-3 was 1.04-gigabase-pair (Gbp) with a heterozygosity rate of 1.62%. De novo assembly yielded a scaffolded 817.7 Mb genome while a reference-guided approach, yielded 843 Mb of genome sequence. The estimated GC content was 34.10%. Variant analysis revealed that BARI Kanthal-3 included 5.7 M (35%) and 10.4 M (65%) simple and heterozygous single nucleotide polymorphisms (SNPs), and about 90% of all these polymorphisms are in inter-genic regions. Through BUSCO assessment, 97.2% of the core genes were represented in the assembly with 1.3% and 1.5% either fragmented or missing, respectively. By comparing identified orthologous gene groups in BARI Kanthal-3 with five closely and one distantly related species of 10,092 common orthogroups were found across the genomes of the six species. The phylogenetic analysis of the shared orthogroups showed that A. heterophyllus was the closest species to BARI Kanthal-3 and orthogroups related to flowering time were found to be more highly prevalent in BARI Kanthal-3 compared to the other Arctocarpus spp. The findings of this study will help better understanding the evolution, domestication, phylogenetic relationships, year-round fruiting of this highly nutritious fruit crop as well as providing a resource for molecular breeding.

Genomic distribution and context dependent functionality of novel WRKY transcription factor binding sites

Background

The WT-boxes NGACTTTN are novel microbe-associated molecular pattern (MAMP)-responsive cis-regulatory sequences. Many of them are uncommon WRKY transcription factor (TF) binding sites.

Results

To understand their functional relevance, a genomic distribution analysis of the 16 possible WT-boxes and a functional analysis of a WT-box rich promoter was done. The genomic distribution analysis shows an enrichment of specific WT-boxes within 500 bp upstream of all Arabidopsis thaliana genes. Those that harbour a T 5′ to the core sequence GACTTT can also be part of the classic WRKY binding site the W-box TTGACT/C. The MAMP-responsive gene ATEP3, a class IV chitinase, harbours seven WT-boxes within its 1000 bp upstream region. In the context of synthetic promoters, the four proximal WT-boxes confer MAMP responsivity while the three WT-boxes further upstream have no effect. Rendering the nucleotides adjacent and in the vicinity of the WT-box core sequence reveals their functional importance for gene expression. A 158 bp long ATEP3 minimal promoter harbouring the two WT-boxes CGACTTTT, confers WT-box-dependent basal and MAMP-responsive reporter gene expression. The ATEP3 gene is a proposed target of WRKY50 and WRKY70. WRKY50 negatively regulates MAMP responsivity of the two WT-boxes CGACTTTT, while WRKY70 activates gene expression in a WT-box dependent manner. Both WRKY factors bind directly to the WT-box CGACTTTT.

Conclusion

In summary, WT-boxes are enriched in promoter regions and comprise novel and uncommon WRKY binding sites required for basal and MAMP-induced gene expression. WT-boxes not being part of a W-box may be a missing link for WRKY target gene prediction when these genes do not harbour a W-box.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08877-y.

Plant_SNP_TATA_Z-Tester: A Web Service That Unequivocally Estimates the Impact of Proximal Promoter Mutations on Plant Gene Expression

Synthetic targeted optimization of plant promoters is becoming a part of progress in mainstream postgenomic agriculture along with hybridization of cultivated plants with wild congeners, as well as marker-assisted breeding. Therefore, here, for the first time, we compiled all the experimental data—on mutational effects in plant proximal promoters on gene expression—that we could find in PubMed. Some of these datasets cast doubt on both the existence and the uniqueness of the sought solution, which could unequivocally estimate effects of proximal promoter mutation on gene expression when plants are grown under various environmental conditions during their development. This means that the inverse problem under study is ill-posed. Furthermore, we found experimental data on in vitro interchangeability of plant and human TATA-binding proteins allowing the application of Tikhonov’s regularization, making this problem well-posed. Within these frameworks, we created our Web service Plant_SNP_TATA_Z-tester and then determined the limits of its applicability using those data that cast doubt on both the existence and the uniqueness of the sought solution. We confirmed that the effects (of proximal promoter mutations on gene expression) predicted by Plant_SNP_TATA_Z-tester correlate statistically significantly with all the experimental data under study. Lastly, we exemplified an application of Plant_SNP_TATA_Z-tester to agriculturally valuable mutations in plant promoters.

A successful hybrid deep learning model aiming at promoter identification

BACKGROUND

The zone adjacent to a transcription start site (TSS), namely, the promoter, is primarily involved in the process of DNA transcription initiation and regulation. As a result, proper promoter identification is critical for further understanding the mechanism of the networks controlling genomic regulation. A number of methodologies for the identification of promoters have been proposed. Nonetheless, due to the great heterogeneity existing in promoters, the results of these procedures are still unsatisfactory. In order to establish additional discriminative characteristics and properly recognize promoters, we developed the hybrid model for promoter identification (HMPI), a hybrid deep learning model that can characterize both the native sequences of promoters and the morphological outline of promoters at the same time. We developed the HMPI to combine a method called the PSFN (promoter sequence features network), which characterizes native promoter sequences and deduces sequence features, with a technique referred to as the DSPN (deep structural profiles network), which is specially structured to model the promoters in terms of their structural profile and to deduce their structural attributes.

RESULTS

The HMPI was applied to human, plant and Escherichia coli K-12 strain datasets, and the findings showed that the HMPI was successful at extracting the features of the promoter while greatly enhancing the promoter identification performance. In addition, after the improvements of synthetic sampling, transfer learning and label smoothing regularization, the improved HMPI models achieved good results in identifying subtypes of promoters on prokaryotic promoter datasets.

CONCLUSIONS

The results showed that the HMPI was successful at extracting the features of promoters while greatly enhancing the performance of identifying promoters on both eukaryotic and prokaryotic datasets, and the improved HMPI models are good at identifying subtypes of promoters on prokaryotic promoter datasets. The HMPI is additionally adaptable to different biological functional sequences, allowing for the addition of new features or models.

Sequence-based evaluation of promoter context for prediction of transcription start sites in Arabidopsis and rice

Genes are transcribed from transcription start sites (TSSs), and their position in a genome is strictly controlled to avoid mis-expression of undesired regions. In this study, we designed and developed a methodology for the evaluation of promoter context, which detects proximal promoter regions from − 200 to − 60 bp relative to a TSS, in Arabidopsis and rice genomes. The method positively evaluates spacer sequences and Regulatory Element Groups, but not core promoter elements like TATA boxes, and is able to predict the position of a TSS within a width of 200 bp. An important feature of the evaluation/prediction method is its independence of the core promoter elements, which was demonstrated by successful prediction of all the TATA, GA, and coreless types of promoters without notable differences in the accuracy of prediction. The positive relationship identified between the evaluation scores and gene expression levels suggests that this method is useful for the evaluation of promoter maturity.

Mechanosensory trichome cells evoke a mechanical stimuli-induced immune response in Arabidopsis thaliana

Perception of pathogen-derived ligands by corresponding host receptors is a pivotal strategy in eukaryotic innate immunity. In plants, this is complemented by circadian anticipation of infection timing, promoting basal resistance even in the absence of pathogen threat. Here, we report that trichomes, hair-like structures on the epidermis, directly sense external mechanical forces, including raindrops, to anticipate pathogen infections in Arabidopsis thaliana. Exposure of leaf surfaces to mechanical stimuli initiates the concentric propagation of intercellular calcium waves away from trichomes to induce defence-related genes. Propagating calcium waves enable effective immunity against pathogenic microbes through the CALMODULIN-BINDING TRANSCRIPTION ACTIVATOR 3 (CAMTA3) and mitogen-activated protein kinases. We propose an early layer of plant immunity in which trichomes function as mechanosensory cells that detect potential risks.

Genome-wide identification and expression profiling of the COBRA-like genes reveal likely roles in stem strength in rapeseed (Brassica napus L.)

The COBRA-like (COBL) genes play key roles in cell anisotropic expansion and the orientation of microfibrils. Mutations in these genes cause the brittle stem and induce pathogen responsive phenotypes in Arabidopsis and several crop plants. In this study, an in silico genome-wide analysis was performed to identify the COBL family members in Brassica. We identified 44, 20 and 23 COBL genes in B. napus and its diploid progenitor species B. rapa and B. oleracea, respectively. All the predicted COBL genes were phylogenetically clustered into two groups: the AtCOB group and the AtCOBL7 group. The conserved chromosome locations of COBLs in Arabidopsis and Brassica, together with clustering, indicated that the expansion of the COBL gene family in B. napus was primarily attributable to whole-genome triplication. Among the BnaCOBLs, 22 contained all the conserved motifs and derived from 9 of 12 subgroups. RNA-seq analysis was used to determine the tissue preferential expression patterns of various subgroups. BnaCOBL9, BnaCOBL35 and BnaCOBL41 were highly expressed in stem with high-breaking resistance, which implies these AtCOB subgroup members may be involved in stem development and stem breaking resistance of rapeseed. Our results of this study may help to elucidate the molecular properties of the COBRA gene family and provide informative clues for high stem-breaking resistance studies.

Disparate Dynamics of Gene Body and cis-Regulatory Element Evolution Illustrated for the Senescence-Associated Cysteine Protease Gene SAG12 of Plants

Gene regulation networks precisely orchestrate the expression of genes that are closely associated with defined physiological and developmental processes such as leaf senescence in plants. The Arabidopsis thaliana senescence-associated gene 12 (AtSAG12) encodes a cysteine protease that is (i) involved in the degradation of chloroplast proteins and (ii) almost exclusively expressed during senescence. Transcription factors, such as WRKY53 and WRKY45, bind to W-boxes in the promoter region of AtSAG12 and play key roles in its activation. Other transcription factors, such as bZIPs, might have accessory functions in their gene regulation, as several A-boxes have been identified and appear to be highly overrepresented in the promoter region compared to the whole genome distribution but are not localized within the regulatory regions driving senescence-associated expression. To address whether these two regulatory elements exhibiting these different properties are conserved in other closely related species, we constructed phylogenetic trees of the coding sequences of orthologs of AtSAG12 and screened their respective 2000 bp promoter regions for the presence of conserved cis-regulatory elements, such as bZIP and WRKY binding sites. Interestingly, the functional relevant upstream located W-boxes were absent in plant species as closely related as Arabidopsis lyrata, whereas an A-box cluster appeared to be conserved in the Arabidopsis species but disappeared in Brassica napus. Several orthologs were present in other species, possibly because of local or whole genome duplication events, but with distinct cis-regulatory sites in different locations. However, at least one gene copy in each family analyzed carried one W-box and one A-box in its promoter. These gene differences in SAG12 orthologs are discussed in the framework of cis- and trans-regulatory factors, of promoter and gene evolution, of genetic variation, and of the enhancement of the adaptability of plants to changing environmental conditions.

A 1232 bp upstream sequence of glutamine synthetase 1b from Eichhornia crassipes is a root-preferential promoter sequence

BACKGROUND

Glutamine synthetase (GS) acts as a key enzyme in plant nitrogen (N) metabolism. It is important to understand the regulation of GS expression in plant. Promoters can initiate the transcription of its downstream gene. Eichhornia crassipes is a most prominent aquatic invasive plant, which has negative effects on environment and economic development. It also can be used in the bioremediation of pollutants present in water and the production of feeding and energy fuel. So identification and characterization of GS promoter in E. crassipes can help to elucidate its regulation mechanism of GS expression and further to control its N metabolism.

RESULTS

A 1232 bp genomic fragment upstream of EcGS1b sequence from E. crassipes (EcGS1b-P) has been cloned, analyzed and functionally characterized. TSSP-TCM software and PlantCARE analysis showed a TATA-box core element, a CAAT-box, root specific expression element, light regulation elements including chs-CMA1a, Box I, and Sp1 and other cis-acting elements in the sequence. Three 5'-deletion fragments of EcGS1b upstream sequence with 400 bp, 600 bp and 900 bp length and the 1232 bp fragment were used to drive the expression of β-glucuronidase (GUS) in tobacco. The quantitative test revealed that GUS activity decreased with the decreasing of the promoter length, which indicated that there were no negative regulated elements in the EcGS1-P. The GUS expressions of EcGS1b-P in roots were significantly higher than those in leaves and stems, indicating EcGS1b-P to be a root-preferential promoter. Real-time Quantitative Reverse Transcription-Polymerase Chain Reaction (qRT-PCR) analysis of EcGS1b gene also showed higher expression in the roots of E.crassipes than in stems and leaves.

CONCLUSIONS

EcGS1b-P is a root-preferential promoter sequence. It can specifically drive the transcription of its downstream gene in root. This study will help to elucidate the regulatory mechanisms of EcGS1b tissue-specific expression and further study its other regulatory mechanisms in order to utilize E.crassipes in remediation of eutrophic water and control its overgrowth from the point of nutrient metabolism.

Genome-wide identification and expression analysis of Arabidopsis GRAM-domain containing gene family in response to abiotic stresses and PGPR treatment

Characterization of stress-responsive genes is important to understand the genomics perspective of stress tolerance. In this purview, several gene-families are being identified and characterized in the model and non-model plant species, which has greatly enhanced the knowledge of molecular intricacies associated with stress tolerance. One such gene family is the GRAM-domain containing which have been found to be upregulated in response to plant growth-promoting rhizobacteria (PGPR) treatment followed by salinity stress. Thus, we aimed at understanding the involvement of GRAM domain-containing proteins in abiotic stress response under the influence of rhizobacteria in Arabidopsis thaliana. The study identified fourteen AtGRAM genes in A. thaliana. Further, comprehensive analyses of domain family including phylogenetic studies, domain architecture, gene structure and genomic composition analysis, promoter analysis, homology modelling, and duplication and divergence rates estimation was performed. RNA-Seq derived expression profiling of AtGRAM genes using GENVESTIGATOR in different stresses, developmental stages and hormonal treatments was performed, followed by qRT-PCR analysis under abiotic stresses in response to PGPR. Altogether, the study provided insights into the structure, organization, and evolutionary properties of AtGRAM gene family. Modulation in expression pattern in response to stresses influenced by PGPR-treatment suggests its multifaceted role in cross-talk among abiotic stresses and phytohormones. Further functional characterization of the selected candidate genes would enable understanding of the precise roles of GRAM-genes underlying stress tolerance.

Molecular basis for neofunctionalization of duplicated E3 ubiquitin ligases underlying adaptation to drought tolerance in Arabidopsis thaliana

Multigene families in plants expanded from ancestral genes via gene duplication mechanisms constitute a significant fraction of the coding genome. Although most duplicated genes are lost over time, many are retained in the genome. Clusters of tandemly arrayed genes are commonly found in the plant genome where they can promote expansion of gene families. In the present study, promoter fusion to the GUS reporter gene was used to examine the promoter architecture of duplicated E3 ligase genes that are part of group C in the Arabidopsis thaliana ATL family. Acquisition of gene expression by AtATL78, possibly generated from defective AtATL81 expression, is described. AtATL78 expression was purportedly enhanced by insertion of a TATA box within the core promoter region after a short tandem duplication that occurred during evolution of Brassicaceae lineages. This gene is associated with an adaptation to drought tolerance of A. thaliana. These findings also suggest duplicated genes could serve as a reservoir of tacit genetic information, and expression of these duplicated genes is activated upon acquisition of core promoter sequences. Remarkably, drought transcriptome profiling in response to rehydration suggests that ATL78-dependent gene expression predominantly affects genes with root-specific activities.

Epigenetic regulation of spurious transcription initiation in Arabidopsis

In plants, epigenetic regulation is critical for silencing transposons and maintaining proper gene expression. However, its impact on the genome-wide transcription initiation landscape remains elusive. By conducting a genome-wide analysis of transcription start sites (TSSs) using cap analysis of gene expression (CAGE) sequencing, we show that thousands of TSSs are exclusively activated in various epigenetic mutants of Arabidopsis thaliana and referred to as cryptic TSSs. Many have not been identified in previous studies, of which up to 65% are contributed by transposons. They possess similar genetic features to regular TSSs and their activation is strongly associated with the ectopic recruitment of RNAPII machinery. The activation of cryptic TSSs significantly alters transcription of nearby TSSs, including those of genes important for development and stress responses. Our study, therefore, sheds light on the role of epigenetic regulation in maintaining proper gene functions in plants by suppressing transcription from cryptic TSSs.

Epigenetic regulation can silence transposons and maintain gene expression. Here the authors survey Arabidopsis mutants defective in epigenetic regulation and show ectopic activation of thousands of cryptic TSSs and altered expression of nearby genes demonstrating the importance of suppressing spurious transcription.

Insights Into Genetic and Molecular Elements for Transgenic Crop Development

Climate change and the exploration of new areas of cultivation have impacted the yields of several economically important crops worldwide. Both conventional plant breeding based on planned crosses between parents with specific traits and genetic engineering to develop new biotechnological tools (NBTs) have allowed the development of elite cultivars with new features of agronomic interest. The use of these NBTs in the search for agricultural solutions has gained prominence in recent years due to their rapid generation of elite cultivars that meet the needs of crop producers, and the efficiency of these NBTs is closely related to the optimization or best use of their elements. Currently, several genetic engineering techniques are used in synthetic biotechnology to successfully improve desirable traits or remove undesirable traits in crops. However, the features, drawbacks, and advantages of each technique are still not well understood, and thus, these methods have not been fully exploited. Here, we provide a brief overview of the plant genetic engineering platforms that have been used for proof of concept and agronomic trait improvement, review the major elements and processes of synthetic biotechnology, and, finally, present the major NBTs used to improve agronomic traits in socioeconomically important crops.

Cytosolic GLUTAMINE SYNTHETASE1;1 Modulates Metabolism and Chloroplast Development in Roots

Nitrogen (N) is an essential macronutrient, and the final form of endogenous inorganic N is ammonium, which is assimilated by Gln synthetase (GS) into Gln. However, how the multiple isoforms of cytosolic GSs contribute to metabolic systems via the regulation of ammonium assimilation remains unclear. In this study, we compared the effects of two rice (Oryza sativa) cytosolic GSs, namely OsGS1;1 and OsGS1;2, on central metabolism in roots using reverse genetics, metabolomic and transcriptomic profiling, and network analyses. We observed (1) abnormal sugar and organic N accumulation and (2) significant up-regulation of genes associated with photosynthesis and chlorophyll biosynthesis in the roots of Osgs1;1 but not Osgs1;2 knockout mutants. Network analysis of the Osgs1;1 mutant suggested that metabolism of Gln was coordinated with the metabolic modules of sugar metabolism, tricarboxylic acid cycle, and carbon fixation. Transcript profiling of Osgs1;1 mutant roots revealed that expression of the rice sigma-factor (OsSIG) genes in the mutants were transiently upregulated. GOLDEN2-LIKE transcription factor-encoding genes, which are involved in chloroplast biogenesis in rice, could not compensate for the lack of OsSIGs in the Osgs1;1 mutant. Microscopic analysis revealed mature chloroplast development in Osgs1;1 roots but not in the roots of Osgs1;2, Osgs1;2-complemented lines, or the wild type. Thus, organic N assimilated by OsGS1;1 affects a broad range of metabolites and transcripts involved in maintaining metabolic homeostasis and plastid development in rice roots, whereas OsGS1;2 has a more specific role, affecting mainly amino acid homeostasis but not carbon metabolism.

The strength and pattern of natural selection on gene expression in rice

Levels of gene expression underpin organismal phenotypes^1,2, but the nature of selection that acts on gene expression and its role in adaptive evolution remain unknown^1,2. Here we assayed gene expression in rice (Oryza sativa)³, and used phenotypic selection analysis to estimate the type and strength of selection on the levels of more than 15,000 transcripts^4,5. Variation in most transcripts appears (nearly) neutral or under very weak stabilizing selection in wet paddy conditions (with median standardized selection differentials near zero), but selection is stronger under drought conditions. Overall, more transcripts are conditionally neutral (2.83%) than are antagonistically pleiotropic⁶ (0.04%), and transcripts that display lower levels of expression and stochastic noise^7-9 and higher levels of plasticity⁹ are under stronger selection. Selection strength was further weakly negatively associated with levels of cis-regulation and network connectivity⁹. Our multivariate analysis suggests that selection acts on the expression of photosynthesis genes^4,5, but that the efficacy of selection is genetically constrained under drought conditions¹⁰. Drought selected for earlier flowering^11,12 and a higher expression of OsMADS18 (Os07g0605200), which encodes a MADS-box transcription factor and is a known regulator of early flowering¹³-marking this gene as a drought-escape gene^11,12. The ability to estimate selection strengths provides insights into how selection can shape molecular traits at the core of gene action.

Insights obtained using different modules of the cotton uceA1.7 promoter

The structure of the cotton uceA1.7 promoter and its modules was analyzed; the potential of their key sequences has been confirmed in different tissues, proving to be a good candidate for the development of new biotechnological tools. Transcriptional promoters are among the primary genetic engineering elements used to control genes of interest (GOIs) associated with agronomic traits. Cotton uceA1.7 was previously characterized as a constitutive promoter with activity higher than that of the constitutive promoter from the Cauliflower mosaic virus (CaMV) 35S gene in various plant tissues. In this study, we generated Arabidopsis thaliana homozygous events stably overexpressing the gfp reporter gene driven by different modules of the uceA1.7 promoter. The expression level of the reporter gene in different plant tissues and the transcriptional stability of these modules was determined compared to its full-length promoter and the 35S promoter. The full-length uceA1.7 promoter exhibited higher activity in different plant tissues compared to the 35S promoter. Two modules of the promoter produced a low and unstable transcription level compared to the other promoters. The other two modules rich in cis-regulatory elements showed similar activity levels to full-length uceA1.7 and 35S promoters but were less stable. This result suggests the location of a minimal portion of the promoter that is required to initiate transcription properly (the core promoter). Additionally, the full-length uceA1.7 promoter containing the 5'-untranslated region (UTR) is essential for higher transcriptional stability in various plant tissues. These findings confirm the potential use of the full-length uceA1.7 promoter for the development of new biotechnological tools (NBTs) to achieve higher expression levels of GOIs in, for example, the root or flower bud for the efficient control of phytonematodes and pest-insects, respectively, in important crops.

Variation of gene expression in plants is influenced by gene architecture and structural properties of promoters

In higher eukaryotes, gene architecture and structural properties of promoters have emerged as significant factors influencing variation in number of transcripts (expression level) and specificity of gene expression in a tissue (expression breadth), which eventually shape the phenotype. In this study, transcriptome data of different tissue types at various developmental stages of A. thaliana, O. sativa, S. bicolor and Z. mays have been used to understand the relationship between properties of gene components and its expression. Our findings indicate that in plants, among all gene architecture and structural properties of promoters, compactness of genes in terms of intron content is significantly linked to gene expression level and breadth, whereas in human an exactly opposite scenario is seen. In plants, for the first time we have carried out a quantitative estimation of effect of a particular trait on expression level and breadth, by using multiple regression analysis and it confirms that intron content of primary transcript (as %) is a powerful determinant of expression breadth. Similarly, further regression analysis revealed that among structural properties of the promoters, stability is negatively linked to expression breadth, while DNase1 sensitivity strongly governs gene expression breadth in monocots and gene expression level in dicots. In addition, promoter regions of tissue specific genes are found to be enriched with TATA box and Y-patch motifs. Finally, multi copy orthologous genes in plants are found to be longer, highly regulated and tissue specific.

Extracellular ATP Shapes a Defense-Related Transcriptome Both Independently and along with Other Defense Signaling Pathways

ATP is not only an essential metabolite of cellular biochemistry but also acts as a signal in the extracellular milieu. In plants, extracellular ATP is monitored by the purinergic receptor P2K1. Recent studies have revealed that extracellular ATP acts as a damage-associated molecular pattern in plants, and its signaling through P2K1 is important for mounting an effective defense response against various pathogenic microorganisms. Biotrophic and necrotrophic pathogens attack plants using different strategies, to which plants respond accordingly with salicylate-based or jasmonate/ethylene-based defensive signaling, respectively. Interestingly, defense mediated by P2K1 is effective against pathogens of both lifestyles, raising the question of the level of interplay between extracellular ATP signaling and that of jasmonate, ethylene, and salicylate. To address this issue, we analyzed ATP-induced transcriptomes in wild-type Arabidopsis (Arabidopsis thaliana) seedlings and mutant seedlings defective in essential components in the signaling pathways of jasmonate, ethylene, and salicylate (classic defense hormones) as well as a mutant and an overexpression line of the P2K1 receptor. We found that P2K1 function is crucial for faithful ATP-induced transcriptional changes and that a subset of genes is more responsive in the P2K1 overexpression line. We also found that more than half of the ATP-responsive genes required signaling by one or more of the pathways for the classical defense hormones, with the jasmonate-based signaling being more critical than others. By contrast, the other ATP-responsive genes were unaffected by deficiencies in signaling for any of the classical defense hormones. These ATP-responsive genes were highly enriched for defense-related Gene Ontology terms. We further tested the ATP-induced genes in knockout mutants of transcription factors, demonstrating that MYCs acting downstream of the jasmonate receptor complex and calmodulin-binding transcription activators are nuclear transducers of P2K1-mediated extracellular ATP signaling.

Dempster-Shafer Theory for the Prediction of Auxin-Response Elements (AuxREs) in Plant Genomes

Auxin is a major regulator of plant growth and development; its action involves transcriptional activation. The identification of Auxin-response element (AuxRE) is one of the most important issues to understand the Auxin regulation of gene expression. Over the past few years, a large number of motif identification tools have been developed. Despite these considerable efforts provided by computational biologists, building reliable models to predict regulatory elements has still been a difficult challenge. In this context, we propose in this work a data fusion approach for the prediction of AuxRE. Our method is based on the combined use of Dempster-Shafer evidence theory and fuzzy theory. To evaluate our model, we have scanning the DORNRÖSCHEN promoter by our model. All proven AuxRE present in the promoter has been detected. At the 0.9 threshold we have no false positive. The comparison of the results of our model and some previous motifs finding tools shows that our model can predict AuxRE more successfully than the other tools and produce less false positive. The comparison of the results before and after combination shows the importance of Dempster-Shafer combination in the decrease of false positive and to improve the reliability of prediction. For an overall evaluation we have chosen to present the performance of our approach in comparison with other methods. In fact, the results indicated that the data fusion method has the highest degree of sensitivity (Sn) and Positive Predictive Value (PPV).

Identification of Leaf Promoters for Use in Transgenic Wheat

Wheat yields have plateaued in recent years and given the growing global population there is a pressing need to develop higher yielding varieties to meet future demand. Genetic manipulation of photosynthesis in elite wheat varieties offers the opportunity to significantly increase yields. However, the absence of a well-defined molecular tool-box of promoters to manipulate leaf processes in wheat hinders advancements in this area. Two promoters, one driving the expression of sedoheptulose-1,7-bisphosphatase (SBPase) and the other fructose-1,6-bisphosphate aldolase (FBPA) from Brachypodium distachyon were identified and cloned into a vector in front of the GUS reporter gene. Both promoters were shown to be functionally active in wheat in both transient assays and in stably transformed wheat plants. Analysis of the stable transformants of wheat (cv. Cadenza) showed that both promoters controlled gus expression throughout leaf development as well as in other green tissues. The availability of these promoters provides new tools for the expression of genes in transgenic wheat leaves and also paves the way for multigene manipulation of photosynthesis to improve yields.

Features of Ppd-B1 expression regulation and their impact on the flowering time of wheat near-isogenic lines

BACKGROUND

Photoperiod insensitive Ppd-1a alleles determine early flowering of wheat. Increased expression of homoeologous Ppd-D1a and Ppd-A1a result from deletions in the promoter region, and elevated expression of Ppd-B1a is determined by an increased copy number.

RESULTS

In this study, using bread wheat cultivars Sonora and PSL2, which contrast in flowering time, and near-isogenic lines resulting from their cross, "Ppd-m" and "Ppd-w" with Ppd-B1a introgressed from Sonora, we investigated the putative factors that influence Ppd-B1a expression. By analyzing the Ppd-B1a three distinct copies, we identified an indel and the two SNPs, which distinguished the investigated allele from other alleles with a copy number variation. We studied the expression of the Ppd-A1, Ppd-B1a, and Ppd-D1 genes along with genes that are involved in light perception (PhyA, PhyB, PhyC) and the flowering initiation (Vrn-1, TaFT1) and discussed their interactions. Expression of Ppd-B1a in the "Ppd-m" line, which flowered four days earlier than "Ppd-w", was significantly higher. We found PhyC to be up-regulated in lines with Ppd-B1a alleles. Expression of PhyC was higher in "Ppd-m". Microsatellite genotyping demonstrated that in the line "Ppd-m", there is an introgression in the pericentromeric region of chromosome 5B from the early flowering parental Sonora, while the "Ppd-w" does not have this introgression. FHY3/FAR1 is known to be located in this region. Expression of the transcription factor FHY3/FAR1 was higher in the "Ppd-m" line than in "Ppd-w", suggesting that FHY3/FAR1 is important for the wheat flowering time and may cause earlier flowering of "Ppd-m" as compared to "Ppd-w".

CONCLUSIONS

We propose that there is a positive bidirectional regulation of Ppd-B1a and PhyC with an FHY3/FAR1 contribution. The bidirectional regulation can be proposed for Ppd-A1a and Ppd-D1a. Using in silico analysis, we demonstrated that the specificity of the Ppd-B1 regulation compared to that of homoeologous genes involves not only a copy number variation but also distinct regulatory elements.

Eukaryotic and prokaryotic promoter databases as valuable tools in exploring the regulation of gene transcription: a comprehensive overview

The complete exploration of the regulation of gene expression remains one of the top-priority goals for researchers. As the regulation is mainly controlled at the level of transcription by promoters, study on promoters and findings are of great importance. This review summarizes forty selected databases that centralize experimental and theoretical knowledge regarding the organization of promoters, interacting transcription factors (TFs) and microRNAs (miRNAs) in many eukaryotic and prokaryotic species. The presented databases offer researchers valuable support in elucidating the regulation of gene transcription.

Genome-wide analysis of cis-regulatory element structure and discovery of motif-driven gene co-expression networks in grapevine

Coordinated transcriptional and metabolic reprogramming ensures a plant’s continued growth and survival under adverse environmental conditions. Transcription factors (TFs) act to modulate gene expression through complex cis-regulatory element (CRE) interactions. Genome-wide analysis of known plant CREs was performed for all currently predicted protein-coding gene promoters in grapevine (Vitis vinifera L.). Many CREs such as abscisic acid (ABA)-responsive, drought-responsive, auxin-responsive, and evening elements, exhibit bona fide CRE properties such as strong position bias towards the transcription start site (TSS) and over-representation when compared with random promoters. Genes containing these CREs are enriched in a large repertoire of plant biological pathways. Large-scale transcriptome analyses also show that these CREs are highly implicated in grapevine development and stress response. Numerous CRE-driven modules in condition-specific gene co-expression networks (GCNs) were identified and many of these modules were highly enriched for plant biological functions. Several modules corroborate known roles of CREs in drought response, pathogen defense, cell wall metabolism, and fruit ripening, whereas others reveal novel functions in plants. Comparisons with Arabidopsis suggest a general conservation in promoter architecture, gene expression dynamics, and GCN structure across species. Systems analyses of CREs provide insights into the grapevine cis-regulatory code and establish a foundation for future genomic studies in grapevine.

Prediction of bipartite transcriptional regulatory elements using transcriptome data of Arabidopsis

In our previous study, a methodology was established to predict transcriptional regulatory elements in promoter sequences using transcriptome data based on a frequency comparison of octamers. Some transcription factors, including the NAC family, cannot be covered by this method because their binding sequences have non-specific spacers in the middle of the two binding sites. In order to remove this blind spot in promoter prediction, we have extended our analysis by including bipartite octamers that are composed of ‘4 bases—a spacer with a flexible length—4 bases’. 8,044 pre-selected bipartite octamers, which had an overrepresentation of specific spacer lengths in promoter sequences and sequences related to core elements removed, were subjected to frequency comparison analysis. Prediction of ER stress-responsive elements in the BiP/BiPL promoter and an ANAC017 target sequence resulted in precise detection of true positives, judged by functional analyses of a reported article and our own in vitro protein–DNA binding assays. These results demonstrate that incorporation of bipartite octamers with continuous ones improves promoter prediction significantly.

Auxin regulates functional gene groups in a fold-change-specific manner in Arabidopsis thaliana roots

Auxin plays a pivotal role in virtually every aspect of plant morphogenesis. It simultaneously orchestrates a diverse variety of processes such as cell wall biogenesis, transition through the cell cycle, or metabolism of a wide range of chemical substances. The coordination principles for such a complex orchestration are poorly understood at the systems level. Here, we perform an RNA-seq experiment to study the transcriptional response to auxin treatment  within gene groups of different biological processes, molecular functions, or cell components in a quantitative fold-change-specific manner. We find for Arabidopsis thaliana roots treated with auxin for 6 h that (i) there are functional groups within which genes respond to auxin with a surprisingly similar fold changes and that (ii) these fold changes vary from one group to another. These findings make it tempting to conjecture the existence of some transcriptional logic orchestrating the coordinated expression of genes within functional groups in a fold-change-specific manner. To obtain some initial insight about this coordinated expression, we performed a motif enrichment analysis and found cis-regulatory elements TBX1-3, SBX, REG, and TCP/site2 as the candidates conferring fold-change-specific responses to auxin in Arabidopsis thaliana.

TSSPlant: a new tool for prediction of plant Pol II promoters

Our current knowledge of eukaryotic promoters indicates their complex architecture that is often composed of numerous functional motifs. Most of known promoters include multiple and in some cases mutually exclusive transcription start sites (TSSs). Moreover, TSS selection depends on cell/tissue, development stage and environmental conditions. Such complex promoter structures make their computational identification notoriously difficult. Here, we present TSSPlant, a novel tool that predicts both TATA and TATA-less promoters in sequences of a wide spectrum of plant genomes. The tool was developed by using large promoter collections from ppdb and PlantProm DB. It utilizes eighteen significant compositional and signal features of plant promoter sequences selected in this study, that feed the artificial neural network-based model trained by the backpropagation algorithm. TSSPlant achieves significantly higher accuracy compared to the next best promoter prediction program for both TATA promoters (MCC≃0.84 and F1-score≃0.91 versus MCC≃0.51 and F1-score≃0.71) and TATA-less promoters (MCC≃0.80, F1-score≃0.89 versus MCC≃0.29 and F1-score≃0.50). TSSPlant is available to download as a standalone program at http://www.cbrc.kaust.edu.sa/download/.

Identification of Arabidopsis genic and non-genic promoters by paired-end sequencing of TSS tags

Information about transcription start sites (TSSs) provides baseline data for the analysis of promoter architecture. In this paper we used paired- and single-end deep sequencing to analyze Arabidopsis TSS tags from several libraries prepared from roots, shoots, flowers and etiolated seedlings. The clustering of approximately 33 million mapped TSS tags led to the identification of 324 461 promoters that covered 79.7% (21 672/27 206) of protein-coding genes in the Arabidopsis genome. In addition we identified intragenic, antisense and orphan promoters that were not associated with any gene models. Of these, intragenic promoters exhibited unique characteristics regarding dinucleotide sequences at TSSs and core promoter element composition, suggesting that these promoters use different mechanisms of transcriptional initiation. An analysis of base composition with regard to promoter position revealed a low GC content throughout the promoter region and several local strand biases that were evident for TATA-type promoters, but not for Coreless-type promoters. Most observed strand biases coincided with strand biases of single nucleotide polymorphism rate. Our analysis also revealed that transcription of a gene is supported by an average of 2.7 genic promoters, among which one specific promoter, designated as a top promoter, substantially determines the expression level of the gene.

Plant Promoter Database (PPDB)

ppdb ( http://ppdb.agr.gifu-u.ac.jp ) is a web-based plant promoter database that provides promoter information of each gene in genomes of Arabidopsis, rice, poplar, and Physcomitrella patens. In this database, recognition of a promoter structure is achieved by annotating genome sequences with our sequence lists of bioinformatically identified octamers for core promoter structure (TATA boxes, Initiators, Y Patches, GA and CA Elements) and regulatory element groups (REGs), together with information of transcription start sites (TSSs) that have been experimentally identified. Our promoter elements are octamer sequences that show strongly biased localization profiles in the promoter region, extracted by the local distribution of short sequence (LDSS) analysis. In addition, REGs are linked with the information of the PLACE database and also with their physiological roles that are predicted using large-scale gene expression data.

Functional Characterization of CYP94-Genes and Identification of a Novel Jasmonate Catabolite in Flowers

Over the past decades much research focused on the biosynthesis of the plant hormone jasmonyl-isoleucine (JA-Ile). While many details about its biosynthetic pathway as well about its physiological function are established nowadays, knowledge about its catabolic fate is still scarce. Only recently, the hormonal inactivation mechanisms became a stronger research focus. Two major pathways have been proposed to inactivate JA-Ile: i) The cleavage of the jasmonyl-residue from the isoleucine moiety, a reaction that is catalyzed by specific amido-hydrolases, or ii), the sequential oxidation of the ω-end of the pentenyl side-chain. This reaction is catalyzed by specific members of the cytochrome P450 (CYP) subfamily CYP94: CYP94B1, CYP94B3 and CYP94C1. In the present study, we further investigated the oxidative fate of JA-Ile by expanding the analysis on Arabidopsis thaliana mutants, lacking only one (cyp94b1, cyp94b2, cyp94b3, cyp94c1), two (cyp94b1xcyp94b2, cyp94b1xcyp94b3, cyp94b2xcyp94b3), three (cyp94b1xcyp94b2xcyp94b3) or even four (cyp94b1xcyp94b2xcyp94b3xcyp94c1) CYP94 functionalities. The results obtained in the present study show that CYP94B1, CYP94B2, CYP94B3 and CYP94C1 are responsible for catalyzing the sequential ω-oxidation of JA-Ile in a semi-redundant manner. While CYP94B-enzymes preferentially hydroxylate JA-Ile to 12-hydroxy-JA-Ile, CYP94C1 catalyzes primarily the subsequent oxidation, yielding 12-carboxy-JA-Ile. In addition, data obtained from investigating the triple and quadruple mutants let us hypothesize that a direct oxidation of unconjugated JA to 12-hydroxy-JA is possible in planta. Using a non-targeted metabolite fingerprinting analysis, we identified unconjugated 12-carboxy-JA as novel jasmonate derivative in floral tissues. Using the same approach, we could show that deletion of CYP94-genes might not only affect JA-homeostasis but also other signaling pathways. Deletion of CYP94B1, for example, led to accumulation of metabolites that may be characteristic for plant stress responses like systemic acquired resistance. Evaluation of the in vivo function of the different CYP94-enzymes on the JA-sensitivity demonstrated that particularly CYP94B-enzymes might play an essential role for JA-response, whereas CYP94C1 might only be of minor importance.

LDSS-P: an advanced algorithm to extract functional short motifs associated with coordinated gene expression

Identifying functional elements in promoter sequences is a major goal in computational and experimental genome biology. Here, we describe an algorithm, Local Distribution of Short Sequences for Prokaryotes (LDSS-P), to identify conserved short motifs located at specific positions in the promoters of co-expressed prokaryotic genes. As a test case, we applied this algorithm to a symbiotic nitrogen-fixing bacterium, Sinorhizobium meliloti. The LDSS-P profiles that overlap with the 5′ section of the extracytoplasmic function RNA polymerase sigma factor RpoE2 consensus sequences displayed a sharp peak between -34 and -32 from TSS positions. The corresponding genes overlap significantly with RpoE2 targets identified from previous experiments. We further identified several groups of genes that are co-regulated with characterized marker genes. Our data indicate that in S. meliloti, and possibly in other Rhizobiaceae species, the master cell cycle regulator CtrA may recognize an expanded motif (AACCAT), which is positionally shifted from the previously reported CtrA consensus sequence in Caulobacter crescentus. Bacterial one-hybrid experiments showed that base substitution in the expanded motif either increase or decrease the binding by CtrA. These results show the effectiveness of LDSS-P as a method to delineate functional promoter elements.

A novel enrichment strategy reveals unprecedented number of novel transcription start sites at single base resolution in a model prokaryote and the gut microbiome

Background

The initiating nucleotide found at the 5’ end of primary transcripts has a distinctive triphosphorylated end that distinguishes these transcripts from all other RNA species. Recognizing this distinction is key to deconvoluting the primary transcriptome from the plethora of processed transcripts that confound analysis of the transcriptome. The currently available methods do not use targeted enrichment for the 5′end of primary transcripts, but rather attempt to deplete non-targeted RNA.

Results

We developed a method, Cappable-seq, for directly enriching for the 5' end of primary transcripts and enabling determination of transcription start sites at single base resolution. This is achieved by enzymatically modifying the 5′ triphosphorylated end of RNA with a selectable tag. We first applied Cappable-seq to E. coli, achieving up to 50 fold enrichment of primary transcripts and identifying an unprecedented 16539 transcription start sites (TSS) genome-wide at single base resolution. We also applied Cappable-seq to a mouse cecum sample and identified TSS in a microbiome.

Conclusions

Cappable-seq allows for the first time the capture of the 5′ end of primary transcripts. This enables a unique robust TSS determination in bacteria and microbiomes.  In addition to and beyond TSS determination, Cappable-seq depletes ribosomal RNA and reduces the complexity of the transcriptome to a single quantifiable tag per transcript enabling digital profiling of gene expression in any microbiome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2539-z) contains supplementary material, which is available to authorized users.

Small Genetic Circuits and MicroRNAs: Big Players in Polymerase II Transcriptional Control in Plants

RNA Polymerase II (Pol II) regulatory cascades involving transcription factors (TFs) and their targets orchestrate the genetic circuitry of every eukaryotic organism. In order to understand how these cascades function, they can be dissected into small genetic networks, each containing just a few Pol II transcribed genes, that generate specific signal-processing outcomes. Small RNA regulatory circuits involve direct regulation of a small RNA by a TF and/or direct regulation of a TF by a small RNA and have been shown to play unique roles in many organisms. Here, we will focus on small RNA regulatory circuits containing Pol II transcribed microRNAs (miRNAs). While the role of miRNA-containing regulatory circuits as modular building blocks for the function of complex networks has long been on the forefront of studies in the animal kingdom, plant studies are poised to take a lead role in this area because of their advantages in probing transcriptional and posttranscriptional control of Pol II genes. The relative simplicity of tissue- and cell-type organization, miRNA targeting, and genomic structure make the Arabidopsis thaliana plant model uniquely amenable for small RNA regulatory circuit studies in a multicellular organism. In this Review, we cover analysis, tools, and validation methods for probing the component interactions in miRNA-containing regulatory circuits. We then review the important roles that plant miRNAs are playing in these circuits and summarize methods for the identification of small genetic circuits that strongly influence plant function. We conclude by noting areas of opportunity where new plant studies are imminently needed.

A cis-regulatory sequence from a short intergenic region gives rise to a strong microbe-associated molecular pattern-responsive synthetic promoter

The high gene density in Arabidopsis thaliana leaves only relatively short intergenic regions for potential cis-regulatory sequences. To learn more about the regulation of genes harbouring only very short upstream intergenic regions, this study investigates a recently identified novel microbe-associated molecular pattern (MAMP)-responsive cis-sequence located within the 101 bp long intergenic region upstream of the At1g13990 gene. It is shown that the cis-regulatory sequence is sufficient for MAMP-responsive reporter gene activity in the context of its native promoter. The 3' UTR of the upstream gene has a quantitative effect on gene expression. In context of a synthetic promoter, the cis-sequence is shown to achieve a strong increase in reporter gene activity as a monomer, dimer and tetramer. Mutation analysis of the cis-sequence determined the specific nucleotides required for gene expression activation. In transgenic A. thaliana the synthetic promoter harbouring a tetramer of the cis-sequence not only drives strong pathogen-responsive reporter gene expression but also shows a high background activity. The results of this study contribute to our understanding how genes with very short upstream intergenic regions are regulated and how these regions can serve as a source for MAMP-responsive cis-sequences for synthetic promoter design.

A One Precursor One siRNA Model for Pol IV-Dependent siRNA Biogenesis

RNA-directed DNA methylation in Arabidopsis thaliana is driven by the plant-specific RNA Polymerase IV (Pol IV). It has been assumed that a Pol IV transcript can give rise to multiple 24-nt small interfering RNAs (siRNAs) that target DNA methylation. Here, we demonstrate that Pol IV-dependent RNAs (P4RNAs) from wild-type Arabidopsis are surprisingly short in length (30 to 40 nt) and mirror 24-nt siRNAs in distribution, abundance, strand bias, and 5'-adenine preference. P4RNAs exhibit transcription start sites similar to Pol II products and are featured with 5'-monophosphates and 3'-misincorporated nucleotides. The 3'-misincorporation preferentially occurs at methylated cytosines on the template DNA strand, suggesting a co-transcriptional feedback to siRNA biogenesis by DNA methylation to reinforce silencing locally. These results highlight an unusual mechanism of Pol IV transcription and suggest a "one precursor, one siRNA" model for the biogenesis of 24-nt siRNAs in Arabidopsis.

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.

NanoCAGE-XL and CapFilter: an approach to genome wide identification of high confidence transcription start sites

BACKGROUND

Identifying the transcription start sites (TSS) of genes is essential for characterizing promoter regions. Several protocols have been developed to capture the 5' end of transcripts via Cap Analysis of Gene Expression (CAGE) or linker-ligation strategies such as Paired-End Analysis of Transcription Start Sites (PEAT), but often require large amounts of tissue. More recently, nanoCAGE was developed for sequencing on the Illumina GAIIx to overcome these difficulties.

RESULTS

Here we present the first publicly available adaptation of nanoCAGE for sequencing on recent ultra-high throughput platforms such as Illumina HiSeq-2000, and CapFilter, a computational pipeline that greatly increases confidence in TSS identification. We report excellent gene coverage, reproducibility, and precision in transcription start site discovery for samples from Arabidopsis thaliana roots.

CONCLUSION

nanoCAGE-XL together with CapFilter allows for genome wide identification of high confidence transcription start sites in large eukaryotic genomes.

Network Analyses Reveal Shifts in Transcript Profiles and Metabolites That Accompany the Expression of SUN and an Elongated Tomato Fruit

SUN controls elongated tomato (Solanum lycopersicum) shape early in fruit development through changes in cell number along the different axes of growth. The gene encodes a member of the IQ domain family characterized by a calmodulin binding motif. To gain insights into the role of SUN in regulating organ shape, we characterized genome-wide transcriptional changes and metabolite and hormone accumulation after pollination and fertilization in wild-type and SUN fruit tissues. Pericarp, seed/placenta, and columella tissues were collected at 4, 7, and 10 d post anthesis. Pairwise comparisons between SUN and the wild type identified 3,154 significant differentially expressed genes that cluster in distinct gene regulatory networks. Gene regulatory networks that were enriched for cell division, calcium/transport, lipid/hormone, cell wall, secondary metabolism, and patterning processes contributed to profound shifts in gene expression in the different fruit tissues as a consequence of high expression of SUN. Promoter motif searches identified putative cis-elements recognized by known transcription factors and motifs related to mitotic-specific activator sequences. Hormone levels did not change dramatically, but some metabolite levels were significantly altered, namely participants in glycolysis and the tricarboxylic acid cycle. Also, hormone and primary metabolite networks shifted in SUN compared with wild-type fruit. Our findings imply that SUN indirectly leads to changes in gene expression, most strongly those involved in cell division, cell wall, and patterning-related processes. When evaluating global coregulation in SUN fruit, the main node represented genes involved in calcium-regulated processes, suggesting that SUN and its calmodulin binding domain impact fruit shape through calcium signaling.

Hormonal Regulation and Expression Profiles of Wheat Genes Involved during Phytic Acid Biosynthesis Pathway

Phytic acid (PA) biosynthesis pathway genes were reported from multiple crop species. PA accumulation was enhanced during grain filling and at that time, hormones like Abscisic acid (ABA) and Gibberellic acid (GA₃) interplay to control the process of seed development. Regulation of wheat PA pathway genes has not yet been reported in seeds. In an attempt to find the clues for the regulation by hormones, the promoter region of wheat PA pathway genes was analyzed for the presence of cis-elements. Multiple cis-elements of those known to be involved for ABA, GA₃, salicylic acid (SA), and cAMP sensing were identified in the promoters of PA pathway genes. Eight genes (TaIMP, TaITPK1-4, TaPLC1, TaIPK2 and TaIPK1) involved in the wheat PA biosynthesis pathway were selected for the expression studies. The temporal expression response was studied in seeds treated with ABA and GA₃ using quantitative real time PCR. Our results suggested that exogenous application of ABA induces few PA pathway genes in wheat grains. Comparison of expression profiles for PA pathway for GA₃ and ABA suggested the antagonistic regulation of certain genes. Additionally, to reveal stress responses of wheat PA pathway genes, expression was also studied in the presence of SA and cAMP. Results suggested SA specific differential expression of few genes, whereas, overall repression of genes was observed in cAMP treated samples. This study is an effort to understand the regulation of PA biosynthesis genes in wheat.

AtmiRNET: a web-based resource for reconstructing regulatory networks of Arabidopsis microRNAs

Compared with animal microRNAs (miRNAs), our limited knowledge of how miRNAs involve in significant biological processes in plants is still unclear. AtmiRNET is a novel resource geared toward plant scientists for reconstructing regulatory networks of Arabidopsis miRNAs. By means of highlighted miRNA studies in target recognition, functional enrichment of target genes, promoter identification and detection of cis- and trans-elements, AtmiRNET allows users to explore mechanisms of transcriptional regulation and miRNA functions in Arabidopsis thaliana, which are rarely investigated so far. High-throughput next-generation sequencing datasets from transcriptional start sites (TSSs)-relevant experiments as well as five core promoter elements were collected to establish the support vector machine-based prediction model for Arabidopsis miRNA TSSs. Then, high-confidence transcription factors participate in transcriptional regulation of Arabidopsis miRNAs are provided based on statistical approach. Furthermore, both experimentally verified and putative miRNA-target interactions, whose validity was supported by the correlations between the expression levels of miRNAs and their targets, are elucidated for functional enrichment analysis. The inferred regulatory networks give users an intuitive insight into the pivotal roles of Arabidopsis miRNAs through the crosstalk between miRNA transcriptional regulation (upstream) and miRNA-mediate (downstream) gene circuits. The valuable information that is visually oriented in AtmiRNET recruits the scant understanding of plant miRNAs and will be useful (e.g. ABA-miR167c-auxin signaling pathway) for further research.

Database URL:http://AtmiRNET.itps.ncku.edu.tw/