Plant promoter prediction with confidence estimation

Isolation and functional characterization of cold-induced gene (AmCIP) promoter from Ammopiptanthus mongolicus

AmCIP is a dehydrin-like protein which involved in abiotic stress tolerance in xerophytes evergreen woody plant A. mongolicus. AmCIP could be induced in the cotyledon and radicle during cold acclimation. To further elucidate the regulation of the upstream region of the gene, we isolated and characterized the promoter of AmCIP. Herein, a 1115 bp 5'-flanking region of AmCIP genomic DNA was isolated and cloned by genome walking from A. mongolicus and the segment sequence was identified as "PrAmCIP" promoter. Analysis of the promoter sequence revealed the presences of some basic cis-acting elements, which were related to various environmental stresses and plant hormones. GUS histochemical staining of transgene tobacco showed that PrAmCIP was induced by 4℃, 55℃, NaCl, mannitol and ABA, whereas it could hardly drive GUS gene expression under normal conditions. Furthermore, we constructed three deletion fragments and genetically transformed them into Arabidopsis thaliana. GUS histochemical staining showed that the MYCATERD1 element of the CP7 fragment (-189 ∼ -1) may be a key element in response to drought. In conclusion, we provide an inducible promoter, PrAmCIP, which can be applied to the development of transgenic plants for abiotic stresse tolerance.

Isolation of Salvia miltiorrhiza Kaurene Synthase-like (KSL) Gene Promoter and Its Regulation by Ethephon and Yeast Extract

The presented study describes the regulation of the promoter region of the Salvia miltiorrhiza kaurene synthase-like gene (SmKSL) by ethylene and yeast extract. The isolated fragment is 897 bp and is composed of a promoter (763 bp), 5'UTR (109 bp), and a short CDS (25 bp). The initial in silico analysis revealed the presence of numerous putative cis-active sites for trans-factors responding to different stress conditions. However, this study examines the influence of ethylene and yeast extract on SmKSL gene expression and tanshinone biosynthesis regulation. The results of 72h RT-PCR indicate an antagonistic interaction between ethylene, provided as ethephon (0.05, 0.10, 0.25, and 0.50 mM), and yeast extract (0.5%) on SmKSL gene expression in callus cultures of S. miltiorrhiza. A similar antagonistic effect was observed on total tanshinone concentration for up to 60 days. Ethylene provided as ethephon (0.05, 0.10, 0.25, and 0.50 mM) is a weak inducer of total tanshinone biosynthesis, increasing them only up to the maximum value of 0.67 ± 0.04 mg g^-1 DW (60-day induction with 0.50 mM ethephon). Among the tanshinones elicited by ethephon, cryptotanshinone (52.21%) dominates, followed by dihydrotanshinone (45.00%) and tanshinone IIA (3.79%). In contrast, the 0.5% yeast extract strongly increases the total tanshinone concentration up to a maximum value of 13.30 ± 1.09 mg g^-1 DW, observed after 50 days of induction. Yeast extract and ethylene appear to activate different fragments of the tanshinone biosynthesis route; hence the primary tanshinones induced by yeast extract were cryptotanshinone (81.42%), followed by dihydrotanshinone (17.06%) and tanshinone IIA (1.52%).

Biological and Molecular Components for Genetically Engineering Biosensors in Plants

Plants adapt to their changing environments by sensing and responding to physical, biological, and chemical stimuli. Due to their sessile lifestyles, plants experience a vast array of external stimuli and selectively perceive and respond to specific signals. By repurposing the logic circuitry and biological and molecular components used by plants in nature, genetically encoded plant-based biosensors (GEPBs) have been developed by directing signal recognition mechanisms into carefully assembled outcomes that are easily detected. GEPBs allow for in vivo monitoring of biological processes in plants to facilitate basic studies of plant growth and development. GEPBs are also useful for environmental monitoring, plant abiotic and biotic stress management, and accelerating design-build-test-learn cycles of plant bioengineering. With the advent of synthetic biology, biological and molecular components derived from alternate natural organisms (e.g., microbes) and/or de novo parts have been used to build GEPBs. In this review, we summarize the framework for engineering different types of GEPBs. We then highlight representative validated biological components for building plant-based biosensors, along with various applications of plant-based biosensors in basic and applied plant science research. Finally, we discuss challenges and strategies for the identification and design of biological components for plant-based biosensors.

Machine learning: its challenges and opportunities in plant system biology

Sequencing technologies are evolving at a rapid pace, enabling the generation of massive amounts of data in multiple dimensions (e.g., genomics, epigenomics, transcriptomic, metabolomics, proteomics, and single-cell omics) in plants. To provide comprehensive insights into the complexity of plant biological systems, it is important to integrate different omics datasets. Although recent advances in computational analytical pipelines have enabled efficient and high-quality exploration and exploitation of single omics data, the integration of multidimensional, heterogenous, and large datasets (i.e., multi-omics) remains a challenge. In this regard, machine learning (ML) offers promising approaches to integrate large datasets and to recognize fine-grained patterns and relationships. Nevertheless, they require rigorous optimizations to process multi-omics-derived datasets. In this review, we discuss the main concepts of machine learning as well as the key challenges and solutions related to the big data derived from plant system biology. We also provide in-depth insight into the principles of data integration using ML, as well as challenges and opportunities in different contexts including multi-omics, single-cell omics, protein function, and protein-protein interaction. KEY POINTS: • The key challenges and solutions related to the big data derived from plant system biology have been highlighted. • Different methods of data integration have been discussed. • Challenges and opportunities of the application of machine learning in plant system biology have been highlighted and discussed.

Methyl Jasmonate Activates the 2C Methyl-D-erithrytol 2,4-cyclodiphosphate Synthase Gene and Stimulates Tanshinone Accumulation in Salvia miltiorrhiza Solid Callus Cultures

The present study characterizes the 5′ regulatory region of the SmMEC gene. The isolated fragment is 1559 bp long and consists of a promoter, 5′UTR and 31 nucleotide 5′ fragments of the CDS region. In silico bioinformatic analysis found that the promoter region contains repetitions of many potential cis-active elements. Cis-active elements associated with the response to methyl jasmonate (MeJa) were identified in the SmMEC gene promoter. Co-expression studies combined with earlier transcriptomic research suggest the significant role of MeJa in SmMEC gene regulation. These findings were in line with the results of the RT-PCR test showing SmMEC gene expression induction after 72 h of MeJa treatment. Biphasic total tanshinone accumulation was observed following treatment of S. miltiorrhiza solid callus cultures with 50–500 μM methyl jasmonate, with peaks observed after 10–20 and 50–60 days. An early peak of total tanshinone concentration (0.08%) occurred after 20 days of 100 μM MeJa induction, and a second, much lower one, was observed after 50 days of 50 μM MeJa stimulation (0.04%). The dominant tanshinones were cryptotanshinone (CT) and dihydrotanshinone (DHT). To better understand the inducing effect of MeJa treatment on tanshinone biosynthesis, a search was performed for methyl jasmonate-responsive cis-active motifs in the available sequences of gene proximal promoters associated with terpenoid precursor biosynthesis. The results indicate that MeJa has the potential to induce a significant proportion of the presented genes, which is in line with available transcriptomic and RT-PCR data.

Critical assessment of computational tools for prokaryotic and eukaryotic promoter prediction

Promoters are crucial regulatory DNA regions for gene transcriptional activation. Rapid advances in next-generation sequencing technologies have accelerated the accumulation of genome sequences, providing increased training data to inform computational approaches for both prokaryotic and eukaryotic promoter prediction. However, it remains a significant challenge to accurately identify species-specific promoter sequences using computational approaches. To advance computational support for promoter prediction, in this study, we curated 58 comprehensive, up-to-date, benchmark datasets for 7 different species (i.e. Escherichia coli, Bacillus subtilis, Homo sapiens, Mus musculus, Arabidopsis thaliana, Zea mays and Drosophila melanogaster) to assist the research community to assess the relative functionality of alternative approaches and support future research on both prokaryotic and eukaryotic promoters. We revisited 106 predictors published since 2000 for promoter identification (40 for prokaryotic promoter, 61 for eukaryotic promoter, and 5 for both). We systematically evaluated their training datasets, computational methodologies, calculated features, performance and software usability. On the basis of these benchmark datasets, we benchmarked 19 predictors with functioning webservers/local tools and assessed their prediction performance. We found that deep learning and traditional machine learning-based approaches generally outperformed scoring function-based approaches. Taken together, the curated benchmark dataset repository and the benchmarking analysis in this study serve to inform the design and implementation of computational approaches for promoter prediction and facilitate more rigorous comparison of new techniques in the future.

Identification of candidate MYB transcription factors that influence CslF6 expression in barley grain

(1,3;1,4)-β-Glucan is a non-cellulosic polysaccharide required for correct barley grain fill and plant development, with industrial relevance in the brewing and the functional food sector. Barley grains contain higher levels of (1,3;1,4)-β-glucan compared to other small grain cereals and this influences their end use, having undesirable effects on brewing and distilling and beneficial effects linked to human health. HvCslF6 is the main gene contributing to (1,3;1,4)-β-glucan biosynthesis in the grain. Here, the transcriptional regulation of HvCslF6 was investigated using an in-silico analysis of transcription factor binding sites (TFBS) in its putative promoter, and functional characterization in a barley protoplast transient expression system. Based on TFBS predictions, TF classes AP2/ERF, MYB, and basic helix-loop-helix (bHLH) were over-represented within a 1,000 bp proximal HvCslF6 promoter region. Dual luciferase assays based on multiple HvCslF6 deletion constructs revealed the promoter fragment driving HvCslF6 expression. Highest HvCslF6 promoter activity was narrowed down to a 51 bp region located -331 bp to -382 bp upstream of the start codon. We combined this with TFBS predictions to identify two MYB TFs: HvMYB61 and HvMYB46/83 as putative activators of HvCslF6 expression. Gene network analyses assigned HvMYB61 to the same co-expression module as HvCslF6 and other primary cellulose synthases (HvCesA1, HvCesA2, and HvCesA6), whereas HvMYB46/83 was assigned to a different module. Based on RNA-seq expression during grain development, HvMYB61 was cloned and tested in the protoplast system. The transient over-expression of HvMYB61 in barley protoplasts suggested a positive regulatory effect on HvCslF6 expression.

Identification of Regulatory SNPs Associated with Vicine and Convicine Content of Vicia faba Based on Genotyping by Sequencing Data Using Deep Learning

Faba bean (Vicia faba) is a grain legume, which is globally grown for both human consumption as well as feed for livestock. Despite its agro-ecological importance the usage of Vicia faba is severely hampered by its anti-nutritive seed-compounds vicine and convicine (V+C). The genes responsible for a low V+C content have not yet been identified. In this study, we aim to computationally identify regulatory SNPs (rSNPs), i.e., SNPs in promoter regions of genes that are deemed to govern the V+C content of Vicia faba. For this purpose we first trained a deep learning model with the gene annotations of seven related species of the Leguminosae family. Applying our model, we predicted putative promoters in a partial genome of Vicia faba that we assembled from genotyping-by-sequencing (GBS) data. Exploiting the synteny between Medicago truncatula and Vicia faba, we identified two rSNPs which are statistically significantly associated with V+C content. In particular, the allele substitutions regarding these rSNPs result in dramatic changes of the binding sites of the transcription factors (TFs) MYB4, MYB61, and SQUA. The knowledge about TFs and their rSNPs may enhance our understanding of the regulatory programs controlling V+C content of Vicia faba and could provide new hypotheses for future breeding programs.

Barley grain (1,3;1,4)-β-glucan content: effects of transcript and sequence variation in genes encoding the corresponding synthase and endohydrolase enzymes

The composition of plant cell walls is important in determining cereal end uses. Unlike other widely consumed cereal grains barley is comparatively rich in (1,3;1,4)-β-glucan, a source of dietary fibre. Previous work showed Cellulose synthase-like genes synthesise (1,3;1,4)-β-glucan in several tissues. HvCslF6 encodes a grain (1,3;1,4)-β-glucan synthase, whereas the function of HvCslF9 is unknown. Here, the relationship between mRNA levels of HvCslF6, HvCslF9, HvGlbI (1,3;1,4)-β-glucan endohydrolase, and (1,3;1,4)-β-glucan content was studied in developing grains of four barley cultivars. HvCslF6 was differentially expressed during mid (8-15 DPA) and late (38 DPA) grain development stages while HvCslF9 transcript was only clearly detected at 8-10 DPA. A peak of HvGlbI expression was detected at 15 DPA. Differences in transcript abundance across the three genes could partially explain variation in grain (1,3;1,4)-β-glucan content in these genotypes. Remarkably narrow sequence variation was found within the HvCslF6 promoter and coding sequence and does not explain variation in (1,3;1,4)-β-glucan content. Our data emphasise the genotype-dependent accumulation of (1,3;1,4)-β-glucan during barley grain development and a role for the balance between hydrolysis and synthesis in determining (1,3;1,4)-β-glucan content, and suggests that other regulatory sequences or proteins are likely to be involved in this trait in developing grain.

Barley grain (1,3;1,4)-β-glucan content: effects of transcript and sequence variation in genes encoding the corresponding synthase and endohydrolase enzymes

The composition of plant cell walls is important in determining cereal end uses. Unlike other widely consumed cereal grains barley is comparatively rich in (1,3;1,4)-β-glucan, a source of dietary fibre. Previous work showed Cellulose synthase-like genes synthesise (1,3;1,4)-β-glucan in several tissues. HvCslF6 encodes a grain (1,3;1,4)-β-glucan synthase, whereas the function of HvCslF9 is unknown. Here, the relationship between mRNA levels of HvCslF6, HvCslF9, HvGlbI (1,3;1,4)-β-glucan endohydrolase, and (1,3;1,4)-β-glucan content was studied in developing grains of four barley cultivars. HvCslF6 was differentially expressed during mid (8-15 DPA) and late (38 DPA) grain development stages while HvCslF9 transcript was only clearly detected at 8-10 DPA. A peak of HvGlbI expression was detected at 15 DPA. Differences in transcript abundance across the three genes could partially explain variation in grain (1,3;1,4)-β-glucan content in these genotypes. Remarkably narrow sequence variation was found within the HvCslF6 promoter and coding sequence and does not explain variation in (1,3;1,4)-β-glucan content. Our data emphasise the genotype-dependent accumulation of (1,3;1,4)-β-glucan during barley grain development and a role for the balance between hydrolysis and synthesis in determining (1,3;1,4)-β-glucan content, and suggests that other regulatory sequences or proteins are likely to be involved in this trait in developing grain.

Machine learning approaches and their current application in plant molecular biology: A systematic review

Machine learning (ML) is a field of artificial intelligence that has rapidly emerged in molecular biology, thus allowing the exploitation of Big Data concepts in plant genomics. In this context, the main challenges are given in terms of how to analyze massive datasets and extract new knowledge in all levels of cellular systems research. In summary, ML techniques allow complex interactions to be inferred in several biological systems. Despite its potential, ML has been underused due to complex computational algorithms and definition terms. Therefore, a systematic review to disentangle ML approaches is relevant for plant scientists and has been considered in this study. We presented the main steps for ML development (from data selection to evaluation of classification/prediction models) with a respective discussion approaching functional genomics mainly in terms of pathogen effector genes in plant immunity. Additionally, we also considered how to access public source databases under an ML framework towards advancing plant molecular biology and introduced novel powerful tools, such as deep learning.

Abscisic Acid Regulates the 3-Hydroxy-3-methylglutaryl CoA Reductase Gene Promoter and Ginsenoside Production in Panax quinquefolium Hairy Root Cultures

Panax quinquefolium hairy root cultures synthesize triterpenoid saponins named ginsenosides, that have multidirectional pharmacological activity. The first rate-limiting enzyme in the process of their biosynthesis is 3-hydroxy-3-methylglutaryl CoA reductase (HMGR). In this study, a 741 bp fragment of the P. quinquefoliumHMGR gene (PqHMGR), consisting of a proximal promoter, 5′UTR (5′ untranslated region) and 5′CDS (coding DNA sequence) was isolated. In silico analysis of an isolated fragment indicated a lack of tandem repeats, miRNA binding sites, and CpG/CpNpG elements. However, the proximal promoter contained potential cis-elements involved in the response to light, salicylic, and abscisic acid (ABA) that was represented by the motif ABRE (TACGTG). The functional significance of ABA on P. quinquefolium HMGR gene expression was evaluated, carrying out quantitative RT-PCR experiments at different ABA concentrations (0.1, 0.25, 0.5, and 1 mg·L⁻¹). Additionally, the effect of abscisic acid and its time exposure on biomass and ginsenoside level in Panax quinquefolium hairy root was examined. The saponin content was determined using HPLC. The 28 day elicitation period with 1 mg·L⁻¹ ABA was the most efficient for Rg2 and Re (17.38 and 1.83 times increase, respectively) accumulation; however, the protopanaxadiol derivative content decreased in these conditions.

GCTTCA as a novel motif for regulating mesocarp-specific expression of the oil palm (Elaeis guineensis Jacq.) stearoyl-ACP desaturase gene

TAAAAT and a novel motif, GCTTCA found in the oil palm stearoyl-ACP desaturase (SAD1) promoter are involved in regulating mesocarp-specific expression. Two key fatty acid biosynthetic genes, stearoyl-ACP desaturase (SAD1), and acyl-carrier protein (ACP3) in Elaeis guineensis (oil palm) showed high level of expression during the period of oil synthesis in the mesocarp [12-19 weeks after anthesis (w.a.a.)] and kernel (12-15 w.a.a.). Both genes are expressed in spear leaves at much lower levels and the expression increased by 1.5-fold to 2.5-fold following treatments with ethylene and abscisic acid (ABA). Both SAD1 and ACP3 promoters contain phytohormone-responsive, light-responsive, abiotic factors/wounding-responsive, endosperm specificity and fruit maturation/ripening regulatory motifs. The activities of the full length and six 5' deletion fragments of the SAD1 promoter were analyzed in transiently transformed oil palm tissues by quantitative β-glucuronidase (GUS) fluorometric assay. The highest SAD1 promoter activity was observed in the mesocarp followed by kernel and the least in the leaves. GUS activity in the D3 deletion construct (- 486 to + 108) was the highest, while the D2 (- 535 to + 108) gave the lowest suggesting the presence of negative cis-acting regulatory element(s) in the deleted - 535 to - 486 (49 bp). It was found that the 49-bp region binds to the nuclear protein extract from mesocarp but not from leaves in electrophoretic mobility shift assay (EMSA). Further fine-tuned analysis of this 49-bp region using truncated DNA led to the identification of GCTTCA as a novel motif in the SAD1 promoter. Interestingly, another known fruit ripening-related motif, LECPLEACS2 (TAAAAT) was found to be required for effective binding of the novel motif to the mesocarp nuclear protein extract.

Nucleotide patterns aiding in prediction of eukaryotic promoters

Computational analysis of promoters is hindered by the complexity of their architecture. In less studied genomes with complex organization, false positive promoter predictions are common. Accurate identification of transcription start sites and core promoter regions remains an unsolved problem. In this paper, we present a comprehensive analysis of genomic features associated with promoters and show that probabilistic integrative algorithms-driven models allow accurate classification of DNA sequence into “promoters” and “non-promoters” even in absence of the full-length cDNA sequences. These models may be built upon the maps of the distributions of sequence polymorphisms, RNA sequencing reads on genomic DNA, methylated nucleotides, transcription factor binding sites, as well as relative frequencies of nucleotides and their combinations. Positional clustering of binding sites shows that the cells of Oryza sativa utilize three distinct classes of transcription factors: those that bind preferentially to the [-500,0] region (188 “promoter-specific” transcription factors), those that bind preferentially to the [0,500] region (282 “5′ UTR-specific” TFs), and 207 of the “promiscuous” transcription factors with little or no location preference with respect to TSS. For the most informative motifs, their positional preferences are conserved between dicots and monocots.

Control of root cap maturation and cell detachment by BEARSKIN transcription factors in Arabidopsis

The root cap supports root growth by protecting the root meristem, sensing gravity and interacting with the rhizosphere through metabolite secretion and cell dispersal. Sustained root cap functions therefore rely on balanced proliferation of proximal stem cells and regulated detachment of distal mature cells. Although the gene regulatory network that governs stem cell activity in the root cap has been extensively studied in Arabidopsis, the mechanisms by which root cap cells mature and detach from the root tip are poorly understood. We performed a detailed expression analysis of three regulators of root cap differentiation, SOMBRERO, BEARSKIN1 and BEARSKIN2, and identified their downstream genes. Our results indicate that expression of BEARSKIN1 and BEARSKIN2 is associated with cell positioning on the root surface. We identified a glycosyl hydrolase 28 (GH28) family polygalacturonase (PG) gene as a direct target of BEARSKIN1. Overexpression and loss-of-function analyses demonstrated that the protein encoded by this PG gene facilitates cell detachment. We thus revealed a molecular link between the key regulators of root cap differentiation and the cellular events underlying root cap-specific functions.

Differential fine-tuning of gene expression regulation in coffee leaves by CcDREB1D promoter haplotypes under water deficit

Despite the importance of the DREB1D gene (also known as CBF4) in plant responses to water deficit and cold stress, studies analysing its regulation by transgenic approaches are lacking. In the current work, a functional study of three CcDREB1D promoter haplotypes (named HP15, HP16 and HP17) isolated from drought-tolerant and drought-sensitive clones of Coffea canephora was carried out in plants of C. arabica stably transformed by Agrobacterium tumefaciens by analysing their ability to regulate the expression of the uidA reporter gene in response to water deficit mimicked by polyethylene glycol (-2.0 MPa) and low relative humidity treatments. A deletion analysis of their corresponding 5'-upstream regions revealed increased specificity of β-glucuronidase activity in the polyethylene glycol and low relative humidity treatments, with high expression in leaf mesophyll and guard cells in full-length constructs. RT-qPCR assays also revealed that the HP16 haplotype (specific to clone tolerant to water deficit) had stronger and earlier activity compared with the HP15 and HP17 haplotypes. As most of the cis-regulatory elements involved in ABA-dependent and -independent networks, tissue specificity and light regulation are common to these haplotypes, we propose that their organization, as well as the nucleic acid polymorphisms present outside these boxes, may play a role in modulating activities of DREB1D promoters in guard cells.

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/.

Recognition of prokaryotic and eukaryotic promoters using convolutional deep learning neural networks

Accurate computational identification of promoters remains a challenge as these key DNA regulatory regions have variable structures composed of functional motifs that provide gene-specific initiation of transcription. In this paper we utilize Convolutional Neural Networks (CNN) to analyze sequence characteristics of prokaryotic and eukaryotic promoters and build their predictive models. We trained a similar CNN architecture on promoters of five distant organisms: human, mouse, plant (Arabidopsis), and two bacteria (Escherichia coli and Bacillus subtilis). We found that CNN trained on sigma70 subclass of Escherichia coli promoter gives an excellent classification of promoters and non-promoter sequences (Sn = 0.90, Sp = 0.96, CC = 0.84). The Bacillus subtilis promoters identification CNN model achieves Sn = 0.91, Sp = 0.95, and CC = 0.86. For human, mouse and Arabidopsis promoters we employed CNNs for identification of two well-known promoter classes (TATA and non-TATA promoters). CNN models nicely recognize these complex functional regions. For human promoters Sn/Sp/CC accuracy of prediction reached 0.95/0.98/0,90 on TATA and 0.90/0.98/0.89 for non-TATA promoter sequences, respectively. For Arabidopsis we observed Sn/Sp/CC 0.95/0.97/0.91 (TATA) and 0.94/0.94/0.86 (non-TATA) promoters. Thus, the developed CNN models, implemented in CNNProm program, demonstrated the ability of deep learning approach to grasp complex promoter sequence characteristics and achieve significantly higher accuracy compared to the previously developed promoter prediction programs. We also propose random substitution procedure to discover positionally conserved promoter functional elements. As the suggested approach does not require knowledge of any specific promoter features, it can be easily extended to identify promoters and other complex functional regions in sequences of many other and especially newly sequenced genomes. The CNNProm program is available to run at web server http://www.softberry.com.

Identification, Characterization, and Functional Validation of Drought-responsive MicroRNAs in Subtropical Maize Inbreds

MicroRNA-mediated gene regulation plays a crucial role in controlling drought tolerance. In the present investigation, 13 drought-associated miRNA families consisting of 65 members and regulating 42 unique target mRNAs were identified from drought-associated microarray expression data in maize and were subjected to structural and functional characterization. The largest number of members (14) was found in the zma-miR166 and zma-miR395 families, with several targets. However, zma-miR160, zma-miR390, zma-miR393, and zma-miR2275 each showed a single target. Twenty-three major drought-responsive cis-regulatory elements were found in the upstream regions of miRNAs. Many drought-related transcription factors, such as GAMYB, HD-Zip III, and NAC, were associated with the target mRNAs. Furthermore, two contrasting subtropical maize genotypes (tolerant: HKI-1532 and sensitive: V-372) were used to understand the miRNA-assisted regulation of target mRNA under drought stress. Approximately 35 and 31% of miRNAs were up-regulated in HKI-1532 and V-372, respectively. The up-regulation of target mRNAs was as high as 14.2% in HKI-1532 but was only 2.38% in V-372. The expression patterns of miRNA-target mRNA pairs were classified into four different types: Type I- up-regulation, Type II- down-regulation, Type III- neutral regulation, and Type IV- opposite regulation. HKI-1532 displayed 46 Type I, 13 Type II, and 23 Type III patterns, whereas V-372 had mostly Type IV interactions (151). A low level of negative regulations of miRNA associated with a higher level of mRNA activity in the tolerant genotype helped to maintain crucial biological functions such as ABA signaling, the auxin response pathway, the light-responsive pathway and endosperm expression under stress conditions, thereby leading to drought tolerance. Our study identified candidate miRNAs and mRNAs operating in important pathways under drought stress conditions, and these candidates will be useful in the development of drought-tolerant maize hybrids.

The pathogen-inducible promoter of defense-related LsGRP1 gene from Lilium functioning in phylogenetically distinct species of plants

A suitable promoter greatly enhances the efficiency of target gene expression of plant molecular breeding and farming; however, only very few promoters are available for economically important non-graminaceous ornamental monocots. In this study, an 868-bp upstream region of defense-related LsGRP1 of Lilium, named PLsGRP1, was cloned by genome walking and proven to exhibit promoter activity in Nicotiana benthamiana and Lilium 'Stargazer' as assayed by agroinfiltration-based β-glucuronidase (GUS) expression system. Many putative biotic stress-, abiotic stress- and physiological regulation-related cis-acting elements were found in PLsGRP1. Serial deletion analysis of PLsGRP1 performed in Nicotiana tabacum var. Wisconsin 38 accompanied with types of treatments indicated that 868-bp PLsGRP1 was highly induced upon pathogen challenges and cold stress while the 131-bp 3'-end region of PLsGRP1 could be dramatically induced by many kinds of abiotic stresses, biotic stresses and phytohormone treatments. Besides, transient GUS expression in a fern, gymnosperms, monocots and dicots revealed good promotor activity of PLsGRP1 in many phylogenetically distinct plant species. Thus, pathogen-inducible PLsGRP1 and its 131-bp 3'-end region are presumed potential as tools for plant molecular breeding and farming.

Light response and potential interacting proteins of a grape flavonoid 3'-hydroxylase gene promoter

Flavonoid 3'-hydroxylase (F3'H), a member of cytochrome P450 protein family, introduces B-ring hydroxyl group in the 3' position of the flavonoid. In this study, the cDNA sequence of a F3'H gene (VviF3'H), which contains an open reading frame of 1530 bp encoding a polypeptide of 509 amino acids, was cloned and characterized from Vitis vinifera L. cv. Cabernet Sauvignon. VviF3'H showed high homology to known F3'H genes, especially F3'Hs from the V. vinifera reference genome (Pinot Noir) and lotus. Expression profiling analysis using real-time PCR revealed that VviF3'H was ubiquitously expressed in all tested tissues including berries, leaves, flowers, roots, stems and tendrils, suggesting its important physiological role in plant growth and development. Moreover, the transcript level of VviF3'H gene in grape berries was relatively higher at early developmental stages and gradually decreased during véraison, and then increased in the mature phase. In addition, the promoter of VviF3'H was isolated by using TAIL-PCR. Yeast one-hybrid screening of the Cabernet Sauvignon cDNA library and subsequent in vivo/vitro validations revealed the interaction between VviF3'H promoter and several transcription factors, including members of HD-Zip, NAC, MYB and EIN families. A transcriptional regulation mechanism of VviF3'H expression is proposed for the first time.

Developing fiber specific promoter-reporter transgenic lines to study the effect of abiotic stresses on fiber development in cotton

Cotton is one of the most important cash crops in US agricultural industry. Environmental stresses, such as drought, high temperature and combination of both, not only reduce the overall growth of cotton plants, but also greatly decrease cotton lint yield and fiber quality. The impact of environmental stresses on fiber development is poorly understood due to technical difficulties associated with the study of developing fiber tissues and lack of genetic materials to study fiber development. To address this important question and provide the need for scientific community, we have generated transgenic cotton lines harboring cotton fiber specific promoter (CFSP)-reporter constructs from six cotton fiber specific genes (Expansin, E6, Rac13, CelA1, LTP, and Fb late), representing genes that are expressed at different stages of fiber development. Individual CFSP::GUS or CFSP::GFP construct was introduced into Coker 312 via Agrobacterium mediated transformation. Transgenic cotton lines were evaluated phenotypically and screened for the presence of selectable marker, reporter gene expression, and insertion numbers. Quantitative analysis showed that the patterns of GUS reporter gene activity during fiber development in transgenic cotton lines were similar to those of the native genes. Greenhouse drought and heat stress study showed a correlation between the decrease in promoter activities and decrease in fiber length, increase in micronaire and changes in other fiber quality traits in transgenic lines grown under stressed condition. These newly developed materials provide new molecular tools for studying the effects of abiotic stresses on fiber development and may be used in study of cotton fiber development genes and eventually in the genetic manipulation of fiber quality.

Identification and Characterization of High Temperature Stress Responsive Novel miRNAs in French Bean (Phaseolus vulgaris)

MicroRNAs are important gene regulators controlling almost all biological and metabolic functions. They elicit their regulatory response through modulation of their target gene expression. In this study, we identified eight novel microRNAs (miRNAs) belonging to four miRNA families and one miR* sequence from the French bean genome which responded to high temperature. The precursor miRNAs varied in length and showed conserved signatures of RNA polymerase II transcripts in their upstream regions. Promoter region analysis indicated the prevalence of MYB and WRKY binding sites emphasizing auto-inhibition of miRNA biogenesis. The genomic organization study revealed the presence of 150 putative regulatory motifs of which 41 are unique. Common motifs shared by miRNAs involved in more than one abiotic stresses were also identified. Further, the miRNA validation was carried out by stem-loop real-time PCR, and the results emphasize that the differential expression of miRNAs confers stress tolerance. Functional analysis revealed that most of the targets represent transcription factors. The results obtained would provide new insights to the complex regulatory mechanism employing small non-coding regulatory RNAs toward stress adaptation.

Isolation and characterization of rubisco small subunit gene promoter from common wheat (Triticum aestivum L.)

Choice of an appropriate promoter is critical to express target genes in intended tissues and developmental stages. However, promoters capable of directing gene expression in specific tissues and stages are not well characterized in monocot species. To identify such a promoter in wheat, this study isolated a partial sequence of the wheat small subunit of RuBisCO (TarbcS) promoter. In silico analysis revealed the presence of elements that are characteristic to rbcS promoters of other, mainly dicot, species. Transient expression of the TarbcS:GUS in immature wheat embryos and tobacco leaves but not in the wheat roots indicate the functionality of the TarbcS promoter fragment in directing the expression of target genes in green plant tissues.

Paired-end analysis of transcription start sites in Arabidopsis reveals plant-specific promoter signatures

Understanding plant gene promoter architecture has long been a challenge due to the lack of relevant large-scale data sets and analysis methods. Here, we present a publicly available, large-scale transcription start site (TSS) data set in plants using a high-resolution method for analysis of 5' ends of mRNA transcripts. Our data set is produced using the paired-end analysis of transcription start sites (PEAT) protocol, providing millions of TSS locations from wild-type Columbia-0 Arabidopsis thaliana whole root samples. Using this data set, we grouped TSS reads into "TSS tag clusters" and categorized clusters into three spatial initiation patterns: narrow peak, broad with peak, and weak peak. We then designed a machine learning model that predicts the presence of TSS tag clusters with outstanding sensitivity and specificity for all three initiation patterns. We used this model to analyze the transcription factor binding site content of promoters exhibiting these initiation patterns. In contrast to the canonical notions of TATA-containing and more broad "TATA-less" promoters, the model shows that, in plants, the vast majority of transcription start sites are TATA free and are defined by a large compendium of known DNA sequence binding elements. We present results on the usage of these elements and provide our Plant PEAT Peaks (3PEAT) model that predicts the presence of TSSs directly from sequence.

NPEST: a nonparametric method and a database for transcription start site prediction

In this paper we present NPEST, a novel tool for the analysis of expressed sequence tags (EST) distributions and transcription start site (TSS) prediction. This method estimates an unknown probability distribution of ESTs using a maximum likelihood (ML) approach, which is then used to predict positions of TSS. Accurate identification of TSS is an important genomics task, since the position of regulatory elements with respect to the TSS can have large effects on gene regulation, and performance of promoter motif-finding methods depends on correct identification of TSSs. Our probabilistic approach expands recognition capabilities to multiple TSS per locus that may be a useful tool to enhance the understanding of alternative splicing mechanisms. This paper presents analysis of simulated data as well as statistical analysis of promoter regions of a model dicot plant Arabidopsis thaliana. Using our statistical tool we analyzed 16520 loci and developed a database of TSS, which is now publicly available at www.glacombio.net/NPEST.

Identification of putative TAL effector targets of the citrus canker pathogens shows functional convergence underlying disease development and defense response

Background

Transcriptional activator-like (TAL) effectors, formerly known as the AvrBs3/PthA protein family, are DNA-binding effectors broadly found in Xanthomonas spp. that transactivate host genes upon injection via the bacterial type three-secretion system. Biologically relevant targets of TAL effectors, i.e. host genes whose induction is vital to establish a compatible interaction, have been reported for xanthomonads that colonize rice and pepper; however, citrus genes modulated by the TAL effectors PthA“s” and PthC“s” of the citrus canker bacteria Xanthomonas citri (Xc) and Xanthomonas aurantifolii pathotype C (XaC), respectively, are poorly characterized. Of particular interest, XaC causes canker disease in its host lemon (Citrus aurantifolia), but triggers a defense response in sweet orange.

Results

Based on, 1) the TAL effector-DNA binding code, 2) gene expression data of Xc and XaC-infiltrated sweet orange leaves, and 3) citrus hypocotyls transformed with PthA2, PthA4 or PthC1, we have identified a collection of Citrus sinensis genes potentially targeted by Xc and XaC TAL effectors. Our results suggest that similar with other strains of Xanthomonas TAL effectors, PthA2 and PthA4, and PthC1 to some extent, functionally converge. In particular, towards induction of genes involved in the auxin and gibberellin synthesis and response, cell division, and defense response. We also present evidence indicating that the TAL effectors act as transcriptional repressors and that the best scoring predicted DNA targets of PthA“s” and PthC“s” in citrus promoters predominantly overlap with or localize near to TATA boxes of core promoters, supporting the idea that TAL effectors interact with the host basal transcriptional machinery to recruit the RNA pol II and start transcription.

Conclusions

The identification of PthA“s” and PthC“s” targets, such as the LOB (LATERAL ORGAN BOUNDARY) and CCNBS genes that we report here, is key for the understanding of the canker symptoms development during host susceptibility, or the defenses of sweet orange against the canker bacteria. We have narrowed down candidate targets to a few, which pointed out the host metabolic pathways explored by the pathogens.

Analysis of promoters of microRNAs from a Glycine max degradome library

OBJECTIVE

MicroRNAs (miRNAs) are genome-encoded, small non-coding RNAs that play important functions in development, biotic and abiotic stress responses, and other processes. Our aim was to explore the regulation of miRNA expression.

METHODS

We used bioinformatics methods to predict the core promoters of 440 miRNAs identified from a soybean (Glycine max) degradome library and to analyze cis-acting elements for 369 miRNAs.

RESULTS

The prediction results showed that 83.86% of the 440 miRNAs contained promoters in their upstream sequences, and 8.64% (38 loci) in their downstream sequences. The distributions of two core promoter elements, TATA-boxes and transcription start sites (TSSs), were similar. The cis-acting elements were examined to provide clues to the function and regulation of spatiotemporal expression of the miRNAs. Analyses of miRNA cis-elements and targets indicated a potential auxin response factor (ARF)- and gibberellin response factor (GARF)-mediated negative feedback loop for miRNA expression.

CONCLUSIONS

The features of miRNAs from a Glycine max degradome library obtained here provide insights into the transcription regulation and functions of miRNAs in soybean.

Rule-based knowledge acquisition method for promoter prediction in human and Drosophila species

The rapid and reliable identification of promoter regions is important when the number of genomes to be sequenced is increasing very speedily. Various methods have been developed but few methods investigate the effectiveness of sequence-based features in promoter prediction. This study proposes a knowledge acquisition method (named PromHD) based on if-then rules for promoter prediction in human and Drosophila species. PromHD utilizes an effective feature-mining algorithm and a reference feature set of 167 DNA sequence descriptors (DNASDs), comprising three descriptors of physicochemical properties (absorption maxima, molecular weight, and molar absorption coefficient), 128 top-ranked descriptors of 4-mer motifs, and 36 global sequence descriptors. PromHD identifies two feature subsets with 99 and 74 DNASDs and yields test accuracies of 96.4% and 97.5% in human and Drosophila species, respectively. Based on the 99- and 74-dimensional feature vectors, PromHD generates several if-then rules by using the decision tree mechanism for promoter prediction. The top-ranked informative rules with high certainty grades reveal that the global sequence descriptor, the length of nucleotide A at the first position of the sequence, and two physicochemical properties, absorption maxima and molecular weight, are effective in distinguishing promoters from non-promoters in human and Drosophila species, respectively.

Efficient regeneration potential is closely related to auxin exposure time and catalase metabolism during the somatic embryogenesis of immature embryos in Triticum aestivum L

Regeneration of cultured tissue is a prerequisite of Agrobacterium- and biolistic-mediated plant transformation. In this study, an efficient protocol for improving wheat (Triticum aestivum L.) immature embryo regeneration was developed. Based on the statistical analysis of embryogenic callus induction efficiency, green spot differentiation efficiency, and plant regeneration efficiency from five wheat accessions, improved culture conditions were found to be more effective for embryogenic callus production than the traditional conditions. Using semi-quantitative reverse transcription polymerase chain reaction, a candidate gene, designated as TaCAT1, which encodes a catalase was identified to have a significant correlation with high-regeneration trait of wheat immature embryos. Three amino acid substitutions were found in TaCAT1 protein between high- and low-regeneration wheat accessions. Hydrogen peroxide content in the cultured calli increased from day 5 to 15, and then decreased sharply on day 20, followed by a second peak on day 25 during regeneration stage. Furthermore, a 3,500-bp 5' flanking region upstream of the first codon ATG of TaCAT1 was isolated using inverse polymerase chain reaction. In silico, analysis revealed that the TaCAT1 promoter contained two regulatory motifs associated with responses to auxin.

The tonoplast intrinsic aquaporin (TIP) subfamily of Eucalyptus grandis: Characterization of EgTIP2, a root-specific and osmotic stress-responsive gene

Aquaporins have important roles in various physiological processes in plants, including growth, development and adaptation to stress. In this study, a gene encoding a root-specific tonoplast intrinsic aquaporin (TIP) from Eucalyptus grandis (named EgTIP2) was investigated. The root-specific expression of EgTIP2 was validated over a panel of five eucalyptus organ/tissues. In eucalyptus roots, EgTIP2 expression was significantly induced by osmotic stress imposed by PEG treatment. Histochemical analysis of transgenic tobacco lines (Nicotiana tabacum SR1) harboring an EgTIP2 promoter:GUS reporter cassette revealed major GUS staining in the vasculature and in root tips. Consistent with its osmotic-stress inducible expression in eucalyptus, EgTIP2 promoter activity was up-regulated by mannitol treatment, but was down-regulated by abscisic acid. Taken together, these results suggest that EgTIP2 might be involved in eucalyptus response to drought. Additional searches in the eucalyptus genome revealed the presence of four additional putative TIP coding genes, which could be individually assigned to the classical TIP1-5 groups.

A Novel Moderate Constitutive Promoter Derived from Poplar (Populus tomentosa Carrière)

A novel sequence that functions as a promoter element for moderate constitutive expression of transgenes, designated as the PtMCP promoter, was isolated from the woody perennial Populus tomentosa. The PtMCP promoter was fused to the GUS reporter gene to characterize its expression pattern in different species. In stable Arabidopsis transformants, transcripts of the GUS reporter gene could be detected by RT-PCR in the root, stem, leaf, flower and silique. Further histochemical and fluorometric GUS activity assays demonstrated that the promoter could direct transgene expression in all tissues and organs, including roots, stems, rosette leaves, cauline leaves and flowers of seedlings and maturing plants. Its constitutive expression pattern was similar to that of the CaMV35S promoter, but the level of GUS activity was significantly lower than in CaMV35S promoter::GUS plants. We also characterized the promoter through transient expression in transgenic tobacco and observed similar expression patterns. Histochemical GUS staining and quantitative analysis detected GUS activity in all tissues and organs of tobacco, including roots, stems, leaves, flower buds and flowers, but GUS activity in PtMCP promoter::GUS plants was significantly lower than in CaMV35S promoter::GUS plants. Our results suggested that the PtMCP promoter from poplar is a constitutive promoter with moderate activity and that its function is presumably conserved in different species. Therefore, the PtMCP promoter may provide a practical choice to direct moderate level constitutive expression of transgenes and could be a valuable new tool in plant genetic engineering.

Intraspecific sequence comparisons reveal similar rates of non-collinear gene insertion in the B and D genomes of bread wheat

BACKGROUND

Polyploidization is considered one of the main mechanisms of plant genome evolution. The presence of multiple copies of the same gene reduces selection pressure and permits sub-functionalization and neo-functionalization leading to plant diversification, adaptation and speciation. In bread wheat, polyploidization and the prevalence of transposable elements resulted in massive gene duplication and movement. As a result, the number of genes which are non-collinear to genomes of related species seems markedly increased in wheat.

RESULTS

We used new-generation sequencing (NGS) to generate sequence of a Mb-sized region from wheat chromosome arm 3DS. Sequence assembly of 24 BAC clones resulted in two scaffolds of 1,264,820 and 333,768 bases. The sequence was annotated and compared to the homoeologous region on wheat chromosome 3B and orthologous loci of Brachypodium distachyon and rice. Among 39 coding sequences in the 3DS scaffolds, 32 have a homoeolog on chromosome 3B. In contrast, only fifteen and fourteen orthologs were identified in the corresponding regions in rice and Brachypodium, respectively. Interestingly, five pseudogenes were identified among the non-collinear coding sequences at the 3B locus, while none was found at the 3DS locus.

CONCLUSION

Direct comparison of two Mb-sized regions of the B and D genomes of bread wheat revealed similar rates of non-collinear gene insertion in both genomes with a majority of gene duplications occurring before their divergence. Relatively low proportion of pseudogenes was identified among non-collinear coding sequences. Our data suggest that the pseudogenes did not originate from insertion of non-functional copies, but were formed later during the evolution of hexaploid wheat. Some evidence was found for gene erosion along the B genome locus.

RNase MRP RNA and RNase P activity in plants are associated with a Pop1p containing complex

RNase P processes the 5'-end of tRNAs. An essential catalytic RNA has been demonstrated in Bacteria, Archaea and the nuclei of most eukaryotes; an organism-specific number of proteins complement the holoenzyme. Nuclear RNase P from yeast and humans is well understood and contains an RNA, similar to the sister enzyme RNase MRP. In contrast, no protein subunits have yet been identified in the plant enzymes, and the presence of a nucleic acid in RNase P is still enigmatic. We have thus set out to identify and characterize the subunits of these enzymes in two plant model systems. Expression of the two known Arabidopsis MRP RNA genes in vivo was verified. The first wheat MRP RNA sequences are presented, leading to improved structure models for plant MRP RNAs. A novel mRNA encoding the central RNase P/MRP protein Pop1p was identified in Arabidopsis, suggesting the expression of distinct protein variants from this gene in vivo. Pop1p-specific antibodies precipitate RNase P activity and MRP RNAs from wheat extracts. Our results provide evidence that in plants, Pop1p is associated with MRP RNAs and with the catalytic subunit of RNase P, either separately or in a single large complex.

Native polyubiquitin promoter of rice provides increased constitutive expression in stable transgenic rice plants

The rice Ubiquitin1 (Ubi1) promoter was tested to evaluate its capacity to express the heterologous gusA gene encoding β-glucuronidase in transgenic rice tissue relative to the commonly used Ubi1 corn promoter and the rice gibberellic acid insensitive (GAI) gene promoter element. Experimental results showed increased expression of gusA gene in rice tissue when driven by the native Ubi1 promoter when compared to the use of corn Ubi1 promoter. Results further indicated that the cis-regulatory elements present in the native promoter element might have been responsible for high expression. However, the gusA gene expression level when driven by the rice GAI promoter was notably lower than both Ubi1 promoters. The present study, thus, for the first time helped to demonstrate that the native Ubi1 promoter is a promising genetic element in transgenic approaches for constitutive expression of any gene in rice tissue.

AtHaspin phosphorylates histone H3 at threonine 3 during mitosis and contributes to embryonic patterning in Arabidopsis

Post-translational histone modifications regulate many aspects of chromosome activity. Threonine 3 of histone H3 is highly conserved, but the significance of its phosphorylation is unclear, and the identity of the corresponding kinase in plants is unknown. Therefore, we characterized the candidate kinase in Arabidopsis thaliana, called AtHaspin. Recombinant AtHaspin in vitro phosphorylates histone H3 at threonine 3. Reduction of H3 threonine 3 phosphorylation level and reduced chromatin condensation in interphase nuclei by AtHaspin RNAi supports the proposition that this kinase is involved in histone H3 phosphorylation in vivo in mitotic cells. In addition, we provide a developmental function for a Haspin kinase. At the whole plant level, altered expression of the kinase induced pleiotropic phenotypes with defects in floral organs and vascular tissue. It reduced fertility and modified adventitious shoot apical meristems that then gave rise to plants with multi-rosettes and multi-shoots. Haspin mutant embryos frequently showed alteration in division plane orientation that could be traced back to the earliest divisions of embryo development, thus Haspin contributes to embryonic patterning.

Specific expression of LATERAL SUPPRESSOR is controlled by an evolutionarily conserved 3' enhancer

Aerial plant architecture is largely based on the activity of axillary meristems (AMs), initiated in the axils of leaves. The Arabidopsis gene LATERAL SUPPRESSOR (LAS), which is expressed in well-defined domains at the adaxial boundary of leaf primordia, is a key regulator of AM formation. The precise definition of organ boundaries is an essential step for the formation of new organs in general and for meristem initiation; however, mechanisms leading to these specific patterns are not well understood. To increase understanding of how the highly specific transcript accumulation in organ boundary regions is established, we investigated the LAS promoter. Analysis of deletion constructs revealed that an essential enhancer necessary for complementation is situated about 3.2 kb downstream of the LAS open reading frame. This enhancer is sufficient to confer promoter specificity as upstream sequences in LAS could be replaced by non-specific promoters, such as the 35S minimal promoter. Further promoter swapping experiments using the PISTILLATA or the full 35S promoter demonstrated that the LAS 3' enhancer also has suppressor functions, largely overwriting the activity of different 5' promoters. Phylogenetic analyses suggest that LAS function and regulation are evolutionarily highly conserved. Homologous elements in downstream regulatory sequences were found in all LAS orthologs, including grasses. Transcomplementation experiments demonstrated the functional conservation of non-coding sequences between Solanum lycopersicum (tomato) and Arabidopsis. In summary, our results show that a highly conserved enhancer/suppressor element is the main regulatory module conferring the boundary-specific expression of LAS.

Prediction of plant promoters based on hexamers and random triplet pair analysis

Background

With an increasing number of plant genome sequences, it has become important to develop a robust computational method for detecting plant promoters. Although a wide variety of programs are currently available, prediction accuracy of these still requires further improvement. The limitations of these methods can be addressed by selecting appropriate features for distinguishing promoters and non-promoters.

Methods

In this study, we proposed two feature selection approaches based on hexamer sequences: the Frequency Distribution Analyzed Feature Selection Algorithm (FDAFSA) and the Random Triplet Pair Feature Selecting Genetic Algorithm (RTPFSGA). In FDAFSA, adjacent triplet-pairs (hexamer sequences) were selected based on the difference in the frequency of hexamers between promoters and non-promoters. In RTPFSGA, random triplet-pairs (RTPs) were selected by exploiting a genetic algorithm that distinguishes frequencies of non-adjacent triplet pairs between promoters and non-promoters. Then, a support vector machine (SVM), a nonlinear machine-learning algorithm, was used to classify promoters and non-promoters by combining these two feature selection approaches. We referred to this novel algorithm as PromoBot.

Results

Promoter sequences were collected from the PlantProm database. Non-promoter sequences were collected from plant mRNA, rRNA, and tRNA of PlantGDB and plant miRNA of miRBase. Then, in order to validate the proposed algorithm, we applied a 5-fold cross validation test. Training data sets were used to select features based on FDAFSA and RTPFSGA, and these features were used to train the SVM. We achieved 89% sensitivity and 86% specificity.

Conclusions

We compared our PromoBot algorithm to five other algorithms. It was found that the sensitivity and specificity of PromoBot performed well (or even better) with the algorithms tested. These results show that the two proposed feature selection methods based on hexamer frequencies and random triplet-pair could be successfully incorporated into a supervised machine learning method in promoter classification problem. As such, we expect that PromoBot can be used to help identify new plant promoters. Source codes and analysis results of this work could be provided upon request.

A soybean β-expansin gene GmEXPB2 intrinsically involved in root system architecture responses to abiotic stresses

Root system architecture responds plastically to some abiotic stresses, including phosphorus (P), iron (Fe) and water deficiency, but its response mechanism is still unclear. We cloned and characterized a vegetative β-expansin gene, GmEXPB2, from a Pi starvation-induced soybean cDNA library. Transient expression of 35S::GmEXPB2-GFP in onion epidermal cells verified that GmEXPB2 is a secretory protein located on the cell wall. GmEXPB2 was found to be primarily expressed in roots, and was highly induced by Pi starvation, and the induction pattern was confirmed by GUS staining in transgenic soybean hairy roots. Results from intact soybean composite plants either over-expressing GmEXPB2 or containing knockdown constructs, showed that GmEXPB2 is involved in hairy root elongation, and subsequently affects plant growth and P uptake, especially at low P levels. The results from a heterogeneous transformation system indicated that over-expressing GmEXPB2 in Arabidopsis increased root cell division and elongation, and enhanced plant growth and P uptake at both low and high P levels. Furthermore, we found that, in addition to Pi starvation, GmEXPB2 was also induced by Fe and mild water deficiencies. Taken together, our results suggest that GmEXPB2 is a critical root β-expansin gene that is intrinsically involved in root system architecture responses to some abiotic stresses, including P, Fe and water deficiency. In the case of Pi starvation responses, GmEXPB2 may enhance both P efficiency and P responsiveness by regulating adaptive changes of the root system architecture. This finding has great agricultural potential for improving crop P uptake on both low-P and P-fertilized soils.