One-step identification of conserved miRNAs, their targets, potential transcription factors and effector genes of complete secondary metabolism pathways after 454 pyrosequencing of calyx cDNAs from the Labiate Salvia sclarea L

Lavandula Species, Their Bioactive Phytochemicals, and Their Biosynthetic Regulation

Lavandula species are one of the most useful aromatic and medicinal plants and have great economic potential. The phytopharmaceutical contribution of the secondary metabolites of the species is unquestionable. Most recent studies have been focusing on the elucidation of the genetic background of secondary metabolite production in lavender species. Therefore, knowledge of not only genetic but especially epigenetic mechanisms for the regulation of secondary metabolites is necessary for the modification of those biosynthesis processes and the understanding of genotypic differences in the content and compositional variability of these products. The review discusses the genetic diversity of Lavandula species in relation to the geographic area, occurrence, and morphogenetic factors. The role of microRNAs in secondary-metabolites biosynthesis is described.

MicroRNAs in Medicinal Plants

Medicinal plant microRNAs (miRNAs) are an endogenous class of small RNA central to the posttranscriptional regulation of gene expression. Biosynthetic research has shown that the mature miRNAs in medicinal plants can be produced from either the standard messenger RNA splicing mechanism or the pre-ribosomal RNA splicing process. The medicinal plant miRNA function is separated into two levels: (1) the cross-kingdom level, which is the regulation of disease-related genes in animal cells by oral intake, and (2) the intra-kingdom level, which is the participation of metabolism, development, and stress adaptation in homologous or heterologous plants. Increasing research continues to enrich the biosynthesis and function of medicinal plant miRNAs. In this review, peer-reviewed papers on medicinal plant miRNAs published on the Web of Science were discussed, covering a total of 78 species. The feasibility of the emerging role of medicinal plant miRNAs in regulating animal gene function was critically evaluated. Staged progress in intra-kingdom miRNA research has only been found in a few medicinal plants, which may be mainly inhibited by their long growth cycle, high demand for growth environment, immature genetic transformation, and difficult RNA extraction. The present review clarifies the research significance, opportunities, and challenges of medicinal plant miRNAs in drug development and agricultural production. The discussion of the latest results furthers the understanding of medicinal plant miRNAs and helps the rational design of the corresponding miRNA/target genes functional modules.

Significance of miRNA in enhancement of flavonoid biosynthesis

Flavonoid metabolism shows very strong plasticity in plant development and coping with the changing environment. Flavonoid biosynthesis is regulated by many metabolic pathways, including transcriptional regulation, post-transcriptional control, post-translational regulationand epigenetic regulation. miRNA is a form of endogenous noncoding single-strand small molecule RNA that primarily regulates the expression of target genes horizontally after transcription through splicing and translational suppression. It also plays an important role in regulating plant growth and development, secondary metabolism and biotic and abiotic stress. miRNA can regulate the formation of flavonoids by acting on structural genes or indirectly by using an MBW transcription complex comprising MYB-bHLH-WD40. This study summarizes the biosynthesis and mechanisms of miRNA, and provides a summary of the mechanisms of miRNAs involved in production of flavonoids, in order to elucidate the biosynthesis pathway and complex regulatory network of plant flavonoids. We aim to provide new insights into improving the content of flavonoid active ingredients in plants.

Biosynthesis of Secondary Metabolites Based on the Regulation of MicroRNAs

MicroRNA (miRNA), a noncoding ribonucleic acid, is considered to be important for the progression of gene expression in plants and animals by rupture or translational repression of targeted mRNAs. Many types of miRNA regulate plant metabolism, growth, and response to biotic and abiotic factors. miRNA characterization helps to expose its function in regulating the process of post-transcriptional genetic regulation. There are a lot of factors associated with miRNA function, but the function of miRNA in the organic synthesis of by-products by natural products is not yet fully elucidated. The current review is aimed at observing and characterizing miRNAs and identifying those involved in the functioning of the biosynthesis of secondary metabolites in plants, with their use in controlled manipulation.

Insights into Health-Promoting Effects of Plant MicroRNAs: A Review

Plant-derived microRNAs (miRNAs) play a significant role in human health and are "dark nutrients", as opposed to traditional plant nutrients, as well as important components of food diversification. Studies have revealed that multiple plant-derived miRNA pathways affect human health. First, plant miRNAs regulate plant growth and development and accumulation of metabolites, which alters the food quality and thus indirectly interferes with the health of the host. Moreover, when absorbed in vivo, some miRNAs may target the host cell mRNAs to affect protein expression. In addition, plant miRNAs target and reshape the human gut microbiota (GM), which interferes with the physiology and metabolism of the host. Therefore, miRNAs play a significant role in the cross-kingdom communication of plants, GM, and the host and in maintaining a balance of the three. Future contributions of plant miRNAs can bring new perspectives and opportunities to better understand food nutrition and health care research, which will facilitate the right exploitation of plant resources.

Identification and profiling of conserved microRNAs in different developmental stages of crown imperial (Fritillaria imperialis L.) using high-throughput sequencing

BACKGROUND

Novel strategies for improvement of ornamental plants and their properties relay on miRNA control of differential plant gene expression modulation. Still, in response to the same abiotic stresses, some conserved miRNA families show different expression patterns in different plant species. In parallel, the use of deep sequencing technologies reveals new levels of complexity of regulatory networks in plants through identification of new miRNAs.

METHODS AND RESULTS

Fritillaria imperialis plants were collected from their natural habitats in Koohrang, Chaharmahal va Bakhtiari, Iran. Several tissues including stamen, pistil, petal, sepal, leaf, stem, bulb and fruit were collected during three developmental stages (stem elongation, flower development and seed head stages). Using RNAseq and qRT-PCR approach, this research revealed 21 conserved miRNAs, matching 15 miRNA families, in Fritilaria imperialis.

CONCLUSIONS

The expression of seven conserved miRNAs (Fim-miR156b, Fim-miR159, Fim-miR166a-5p, Fim-miR169d-5p, Fim-miR171c, Fim-miR393 and Fim-miR396e-3p) was further investigated in different tissues and three developmental stages, suggesting different roles for these miRNAs during growth and development of crown imperial. Gained knowledge from this research can open the door to find efficient ways to secure crown imperial survival, preservation and utilization and if proven useful may be applied in other plant species as well.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Phyto-miRNAs-based regulation of metabolites biosynthesis in medicinal plants

Medicinal plants, are known to produce a wide range of plant secondary metabolites (PSMs) applied as insecticides, drugs, dyes and toxins in agriculture, medicine, industry and bio-warfare plus bio-terrorism, respectively. However, production of PSMs is usually in small quantities, so we need to find novel ways to increase both quantity and quality of them. Fortunately, biotechnology suggests several options through which secondary metabolism in plants can be engineered in innovative ways to: 1) over-produce the useful metabolites, 2) down-produce the toxic metabolites, 3) produce the new metabolites. Among the ways, RNA interference (RNAi) technology which involves gene-specific regulation by small non-coding RNAs (sncRNAs) have been recently emerged as a promising tool for plant biotechnologist, not only to decipher the function of plant genes, but also for development of the plants with improved and novel traits through manipulation of both desirable and undesirable genes. Among sncRNAs, miRNAs have been recorded various regulatory roles in plants such as development, signal transduction, response to environmental stresses, metabolism. Certainly, the use of miRNAs in metabolic engineering requires identification of miRNAs involved in metabolites biosynthesis, understanding of the biosynthetic pathways, as well as the identification of key points of the pathways in which the miRNAs have their own effect. Thus, we firstly consider these three issues on metabolic engineering of medicinal plants. Our review shows, application of miRNAs can open a novel perspective to metabolic engineering of medicinal plants.

Increasing the synthesis of bioactive abietane diterpenes in Salvia sclarea hairy roots by elicited transcriptional reprogramming

Transcriptional activation of genes belonging to the plastidial MEP-derived isoprenoid pathway by elicitation with methyl jasmonate and coronatine enhanced the content of bioactive abietane diterpenes in Salvia sclarea hairy roots. We have shown that aethiopinone, an abietane diterpene synthesized in Salvia sclarea roots is cytotoxic and induces apoptosis in human melanoma cells. To develop a production platform for this compound and other abietane diterpenes, hairy root technology was combined with the elicitation of methyl jasmonate (MeJA) or the phytotoxin coronatine (Cor). Both MeJA and Cor induced a significant accumulation of aethiopinone, but prolonged exposure to MeJA irremediably caused inhibition of hairy root growth, which was unaffected by Cor treatment. Considering together the fold increase in aethiopinone content and the final hairy root biomass, the best combination was a Cor treatment for 28 days, which allowed to obtain up to 105.34 ± 2.30 mg L^-1 of this compound to be obtained, corresponding to a 24-fold increase above the basal content in untreated hairy roots. MeJA or Cor elicitation also enhanced the synthesis of other bioactive abietane-quinone diterpenes. The elicitor-dependent steering effect was due to a coordinated transcriptional activation of several biosynthetic genes belonging to the plastidial MEP-derived isoprenoid pathway. High correlations between aethiopinone content and MeJA or Cor-elicited level of gene transcripts were found for DXS2 (r ² = 0.99), DXR (r ² = 0.99), and GGPPS (r ² = 0.98), encoding enzymes acting upstream of GGPP, the common precursor of diterpenes and other plastidial-derived terpenes, as well as CPPS (r ² = 0.99), encoding the enzyme involved in the first cyclization steps leading to copalyl-diphosphate, the precursor of abietane-like diterpenes. These results point to these genes as possible targets of metabolic engineering approaches to establish a more efficient production platform for such promising anti-proliferative plant-derived compounds.

Suppression Substractive Hybridization and NGS Reveal Differential Transcriptome Expression Profiles in Wayfaring Tree (Viburnum lantana L.) Treated with Ozone

Tropospheric ozone (O3) is a global air pollutant that causes high economic damages by decreasing plant productivity. It enters the leaves through the stomata, generates reactive oxygen species, which subsequent decrease in photosynthesis, plant growth, and biomass accumulation. In order to identify genes that are important for conferring O3 tolerance or sensitivity to plants, a suppression subtractive hybridization analysis was performed on the very sensitive woody shrub, Viburnum lantana, exposed to chronic O3 treatment (60 ppb, 5 h d(-1) for 45 consecutive days). Transcript profiling and relative expression assessment were carried out in asymptomatic leaves, after 15 days of O3 exposure. At the end of the experiment symptoms were observed on all treated leaves and plants, with an injured leaf area per plant accounting for 16.7% of the total surface. Cloned genes were sequenced by 454-pyrosequencing and transcript profiling and relative expression assessment were carried out on sequenced reads. A total of 38,800 and 12,495 high quality reads obtained in control and O3-treated libraries, respectively (average length of 319 ± 156.7 and 255 ± 107.4 bp). The Ensembl transcriptome yielded a total of 1241 unigenes with a total sequence length of 389,126 bp and an average length size of 389 bp (guanine-cytosine content = 49.9%). mRNA abundance was measured by reads per kilobase per million and 41 and 37 ensembl unigenes showed up- and down-regulation respectively. Unigenes functionally associated to photosynthesis and carbon utilization were repressed, demonstrating the deleterious effect of O3 exposure. Unigenes functionally associated to heat-shock proteins and glutathione were concurrently induced, suggesting the role of thylakoid-localized proteins and antioxidant-detoxification pathways as an effective strategy for responding to O3. Gene Ontology analysis documented a differential expression of co-regulated transcripts for several functional categories, including specific transcription factors (MYB and WRKY). This study demonstrates that a complex sequence of events takes place in the cells at intracellular and membrane level following O3 exposure and elucidates the effects of this oxidative stress on the transcriptional machinery of the non-model plant species V. lantana, with the final aim to provide the molecular supportive knowledge for the use of this plant as O3-bioindicator.

New insights of medicinal plant therapeutic activity-The miRNA transfer

MicroRNA (miRNA) has become the spotlight of the biomedical research around the world and is considered to be a major post-transcriptional gene regulator. This small, endogenous RNA (21-25 nucleotides long) plays an important role by targeting specific mRNAs in plants, animals and humans. Herbal medicine has been used for thousands of years, however little is known about its molecular mechanism of action. Since the discovery of plant miRNA in human tissue and sera after ingestion, the connection between the two kingdoms is presented under a new perspective. Forward pharmacology, such as miRNA therapeutics could be the next best step toward identifying novel therapeutic options involving medicinal plants. Besides reporting the latest findings regarding the cross-kingdom transfer of miRNA and its therapeutic application, this review can inform further investigations that could lead to a modern definition of herbal medicine.

New opportunities for the regulation of secondary metabolism in plants: focus on microRNAs

Plant cell cultures are of particular interest in industrial applications as a source of biologically active substances. It is difficult, however, to achieve stable production of secondary metabolites for many plant cell cultures using classical techniques. Novel approaches should be developed for removal of the inhibitor blocks that prevent pathway activation and shift the regulatory balance to the activation of entire biosynthetic pathways. MicroRNAs (miRNAs) are small RNAs that play important regulatory roles in various biological processes. Only recently miRNAs have been demonstrated as active in secondary metabolism regulation. In this work, we summarize recent data on the emerging approaches based on regulation of secondary metabolism by miRNAs.

Mining NGS transcriptomes for miRNAs and dissecting their role in regulating growth, development, and secondary metabolites production in different organs of a medicinal herb, Picrorhiza kurroa

This study is the first endeavor on mining of miRNAs and analyzing their involvement in development and secondary metabolism of an endangered medicinal herb Picrorhiza kurroa (P. kurroa ). miRNAs are ubiquitous non-coding RNA species that target complementary sequences of mRNA and result in either translational repression or target degradation in eukaryotes. The role of miRNAs has not been investigated in P. kurroa which is a medicinal herb of industrial value due to the presence of secondary metabolites, picroside-I and picroside-II. Computational identification of miRNAs was done in 6 transcriptomes of P. kurroa generated from root, shoot, and stolon organs varying for growth, development, and culture conditions. All available plant miRNA entries were retrieved from miRBase and used as backend datasets to computationally identify conserved miRNAs in transcriptome data sets. Total 18 conserved miRNAs were detected in P. kurroa followed by target prediction and functional annotation which suggested their possible role in controlling various biological processes. Validation of miRNA and expression analysis by qRT-PCR and 5' RACE revealed that miRNA-4995 has a regulatory role in terpenoid biosynthesis ultimately affecting the production of picroside-I. miR-5532 and miR-5368 had negligible expression in field-grown samples as compared to in vitro-cultured samples suggesting their role in regulating P. kurroa growth in culture conditions. The study has thus identified novel functions for existing miRNAs which can be further validated for their potential regulatory role.

Exploring triacylglycerol biosynthetic pathway in developing seeds of Chia (Salvia hispanica L.): a transcriptomic approach

Chia (Salvia hispanica L.), a member of the mint family (Lamiaceae), is a rediscovered crop with great importance in health and nutrition and is also the highest known terrestrial plant source of heart-healthy omega-3 fatty acid, alpha linolenic acid (ALA). At present, there is no public genomic information or database available for this crop, hindering research on its genetic improvement through genomics-assisted breeding programs. The first comprehensive analysis of the global transcriptome profile of developing Salvia hispanica L. seeds, with special reference to lipid biosynthesis is presented in this study. RNA from five different stages of seed development was extracted and sequenced separately using the Illumina GAIIx platform. De novo assembly of processed reads in the pooled transcriptome using Trinity yielded 76,014 transcripts. The total transcript length was 66,944,462 bases (66.9 Mb), with an average length of approximately 880 bases. In the molecular functions category of Gene Ontology (GO) terms, ATP binding and nucleotide binding were found to be the most abundant and in the biological processes category, the metabolic process and the regulation of transcription-DNA-dependent and oxidation-reduction process were abundant. From the EuKaryotic Orthologous Groups of proteins (KOG) classification, the major category was “Metabolism” (31.97%), of which the most prominent class was ‘carbohydrate metabolism and transport’ (5.81% of total KOG classifications) followed by ‘secondary metabolite biosynthesis transport and catabolism’ (5.34%) and ‘lipid metabolism’ (4.57%). A majority of the candidate genes involved in lipid biosynthesis and oil accumulation were identified. Furthermore, 5596 simple sequence repeats (SSRs) were identified. The transcriptome data was further validated through confirmative PCR and qRT-PCR for select lipid genes. Our study provides insight into the complex transcriptome and will contribute to further genome-wide research and understanding of chia. The identified novel UniGenes will facilitate gene discovery and creation of genomic resource for this crop.

Temporal transcriptome changes induced by methyl jasmonate in Salvia sclarea

Salvia sclarea is a traditional medicinal and aromatic plant that grows in Europe and produces various economically important compounds, including phenylpropanoid derivatives and terpenoids. Methyl jasmonate (MeJA) is commonly used to elicit plant stress responses. However, how MeJA enhances production of secondary metabolites in S. sclarea is not well understood. We performed a genome-wide analysis of temporal gene expression in S. sclarea leaves and roots. The transcriptome profiles 0, 10 and 26 h after MeJA treatment were analyzed by Illumina RNA-Seq. A total of 16,142 isogenes (average length 866bp; N50 1035bp) were obtained by de novo assembly of 35,757,567 raw sequencing reads. When these sequencing reads were mapped onto the assembled Unigenes, 3236, 2792 and 798 Unigenes were found to be expressed differentially between 0 and 10h, 0 and 26 h, and 10 and 26h, respectively. These included many secondary metabolite biosynthesis, stress and defense-related genes. A qRT-PCR analysis confirmed the expression profiles of selected differentially expressed genes (DEGs) revealed by RNA-Seq data, and also extended our analysis of differential gene expression to 73 h. Our investigations revealed temporal differences in the responses of S. sclarea to MeJA treatment. MeJA treatment induced the expression of a large number of genes involved in phenylpropanoid biosynthesis, especially between 0 and 10h, and 0 and 26 h. Additionally, many genes encoding transcription factors, cytochrome P450s, glycosyltransferases, methyltransferases and transporters were shown to respond to MeJA elicitation. DEGs related to structural molecule activity and cell death showed a significant temporal variation. A chromatographic analysis of metabolites at 26h, 73h and six days after MeJA treatment indicated that these transcriptomic changes precede MeJA-induced changes in secondary metabolite content. This study sheds light on the molecular mechanisms of MeJA elicitation and is helpful in understanding how exogenous MeJA treatment mediates extensive plant transcriptome reprogramming/remodeling. Our results can be utilized to characterize genes related to secondary metabolism and their regulation, and in breeding S. sclarea for desirable chemotypes.

Transcriptome analysis in cotton boll weevil (Anthonomus grandis) and RNA interference in insect pests

Cotton plants are subjected to the attack of several insect pests. In Brazil, the cotton boll weevil, Anthonomus grandis, is the most important cotton pest. The use of insecticidal proteins and gene silencing by interference RNA (RNAi) as techniques for insect control are promising strategies, which has been applied in the last few years. For this insect, there are not much available molecular information on databases. Using 454-pyrosequencing methodology, the transcriptome of all developmental stages of the insect pest, A. grandis, was analyzed. The A. grandis transcriptome analysis resulted in more than 500.000 reads and a data set of high quality 20,841 contigs. After sequence assembly and annotation, around 10,600 contigs had at least one BLAST hit against NCBI non-redundant protein database and 65.7% was similar to Tribolium castaneum sequences. A comparison of A. grandis, Drosophila melanogaster and Bombyx mori protein families’ data showed higher similarity to dipteran than to lepidopteran sequences. Several contigs of genes encoding proteins involved in RNAi mechanism were found. PAZ Domains sequences extracted from the transcriptome showed high similarity and conservation for the most important functional and structural motifs when compared to PAZ Domains from 5 species. Two SID-like contigs were phylogenetically analyzed and grouped with T. castaneum SID-like proteins. No RdRP gene was found. A contig matching chitin synthase 1 was mined from the transcriptome. dsRNA microinjection of a chitin synthase gene to A. grandis female adults resulted in normal oviposition of unviable eggs and malformed alive larvae that were unable to develop in artificial diet. This is the first study that characterizes the transcriptome of the coleopteran, A. grandis. A new and representative transcriptome database for this insect pest is now available. All data support the state of the art of RNAi mechanism in insects.

Trans-acting small interfering RNA4: key to nutraceutical synthesis in grape development?

The facility and versatility of microRNAs (miRNAs) to evolve and change likely underlies how they have become dominant constituents of eukaryotic genomes. In this opinion article I propose that trans-acting small interfering RNA gene 4 (TAS4) evolution may be important for biosynthesis of polyphenolics, arbuscular symbiosis, and bacterial pathogen etiologies. Expression-based and phylogenetic evidence shows that TAS4 targets two novel grape (Vitis vinifera L.) MYB transcription factors (VvMYBA6, VvMYBA7) that spawn phased small interfering RNAs (siRNAs) which probably function in nutraceutical bioflavonoid biosynthesis and fruit development. Characterization of the molecular mechanisms of TAS4 control of plant development and integration into biotic and abiotic stress- and nutrient-signaling regulatory networks has applicability to molecular breeding and the development of strategies for engineering healthier foods.

Combining gene expression and genetic analyses to identify candidate genes involved in cold responses in pea

Cold stress affects plant growth and development. In order to better understand the responses to cold (chilling or freezing tolerance), we used two contrasted pea lines. Following a chilling period, the Champagne line becomes tolerant to frost whereas the Terese line remains sensitive. Four suppression subtractive hybridisation libraries were obtained using mRNAs isolated from pea genotypes Champagne and Terese. Using quantitative polymerase chain reaction (qPCR) performed on 159 genes, 43 and 54 genes were identified as differentially expressed at the initial time point and during the time course study, respectively. Molecular markers were developed from the differentially expressed genes and were genotyped on a population of 164 RILs derived from a cross between Champagne and Terese. We identified 5 candidate genes colocalizing with 3 different frost damage quantitative trait loci (QTL) intervals and a protein quantity locus (PQL) rich region previously reported. This investigation revealed the role of constitutive differences between both genotypes in the cold responses, in particular with genes related to glycine degradation pathway that could confer to Champagne a better frost tolerance. We showed that freezing tolerance involves a decrease of expression of genes related to photosynthesis and the expression of a gene involved in the production of cysteine and methionine that could act as cryoprotectant molecules. Although it remains to be confirmed, this study could also reveal the involvement of the jasmonate pathway in the cold responses, since we observed that two genes related to this pathway were mapped in a frost damage QTL interval and in a PQL rich region interval, respectively.

Cross-genera transferable e-microsatellite markers for 12 genera of the Lamiaceae family

BACKGROUND

The Lamiaceae family contains many high-valued medicinal, aromatic and ornamental plant species. Several members of the genera in this family are under heavy pressure of collection for commercial use. DNA markers such as microsatellites could be used to identify commercially important genotypes and to select high-yielding ones for development of new varieties.

RESULTS

A total of 12,432 expressed sequence tags (ESTs) from Salvia fruticosa, S. miltiorrhiza, S. sclarea and Stenogyne rugosa were analyzed. A total of 6216 ESTs were found to be unique according the redundancy test used. Results of this study indicated that the use of redundant ESTs in comparison to non-redundant ESTs was advantageous in terms of higher cross-genera transferability of the markers. A total of 75 EST-microsatellite primer pairs were tested using two different polymerase chain reaction amplification profiles and 52 were found to be cross-genera transferable. Cross-genera transferability of the e-microsatellite primer pairs varied from one species to 12 species tested. It was noted that cross-genera transferability of e-microsatellite primer pairs decreased as the evolutionary distance between the sources and target species increased.

CONCLUSION

This study indicated that EST resources from Salvia spp. and Stenogyne rugosa could be successfully used to identify cross-genera transferable e-microsatellite markers for uncharacterized genomes of the genera in the Lamiaceae family. These e-microsatellite markers could allow one to perform comparative analyses of population structure and genomic studies, and facilitate comparative linkage mapping in the genera studied. E-microsatellite primer pairs reported in this manuscript are equivalent to a total of 135 e-microsatellite primer pairs since many e-microsatellite primer pairs show cross-genera transferability.

Characterization of microRNAs expressed during secondary wall biosynthesis in Acacia mangium

MicroRNAs (miRNAs) play critical regulatory roles by acting as sequence specific guide during secondary wall formation in woody and non-woody species. Although thousands of plant miRNAs have been sequenced, there is no comprehensive view of miRNA mediated gene regulatory network to provide profound biological insights into the regulation of xylem development. Herein, we report the involvement of six highly conserved amg-miRNA families (amg-miR166, amg-miR172, amg-miR168, amg-miR159, amg-miR394, and amg-miR156) as the potential regulatory sequences of secondary cell wall biosynthesis. Within this highly conserved amg-miRNA family, only amg-miR166 exhibited strong differences in expression between phloem and xylem tissue. The functional characterization of amg-miR166 targets in various tissues revealed three groups of HD-ZIP III: ATHB8, ATHB15, and REVOLUTA which play pivotal roles in xylem development. Although these three groups vary in their functions, -psRNA target analysis indicated that miRNA target sequences of the nine different members of HD-ZIP III are always conserved. We found that precursor structures of amg-miR166 undergo exhaustive sequence variation even within members of the same family. Gene expression analysis showed three key lignin pathway genes: C4H, CAD, and CCoAOMT were upregulated in compression wood where a cascade of miRNAs was downregulated. This study offers a comprehensive analysis on the involvement of highly conserved miRNAs implicated in the secondary wall formation of woody plants.

Extracellular localization of the diterpene sclareol in clary sage (Salvia sclarea L., Lamiaceae)

Sclareol is a high-value natural product obtained by solid/liquid extraction of clary sage (Salvia sclarea L.) inflorescences. Because processes of excretion and accumulation of this labdane diterpene are unknown, the aim of this work was to gain knowledge on its sites of accumulation in planta. Samples were collected in natura or during different steps of the industrial process of extraction (steam distillation and solid/liquid extraction). Samples were then analysed with a combination of complementary analytical techniques (gas chromatography coupled to a mass spectrometer, polarized light microscopy, environmental scanning electron microscopy, two-photon fluorescence microscopy, second harmonic generation microscopy). According to the literature, it is hypothesized that sclareol is localized in oil pockets of secretory trichomes. This study demonstrates that this is not the case and that sclareol accumulates in a crystalline epicuticular form, mostly on calyces.

Overexpression of Arabidopsis miR157b induces bushy architecture and delayed phase transition in Torenia fournieri

miR156/157 is a small RNA molecule that is highly conserved among various plant species. Overexpression of miR156/157 has been reported to induce bushy architecture and delayed phase transition in several plant species. To investigate the effect of miR157 overexpression in a horticultural plant, and to explore the applicability of miRNA to molecular breeding, we introduced Arabidopsis MIR157b (AtMIR157b) into torenia (Torenia fournieri). The resulting 35S:AtMIR157b plants showed a high degree of branching along with small leaves, which resembled miR156/157-overexpressing plants of other species. We also isolated torenia SBP-box genes with target miR156/157 sequences and confirmed that their expression was selectively downregulated in 35S:AtMIR157b plants. The reduced accumulation of mRNA was probably due to sequence specificity. Moreover, expression of torenia homologs of the SBP-box protein-regulated genes TfLFY and TfMIR172 was also reduced by AtmiR157 overexpression. These findings suggest that the molecular mechanisms of miR156/157 regulation are conserved between Arabidopsis and torenia. The bushy architecture and small leaves of 35S:AtMIR157b torenia plants could be applied in molecular breeding of various horticultural plants as well as for increasing biomass and crop production.

Discovery and functional characterization of two diterpene synthases for sclareol biosynthesis in Salvia sclarea (L.) and their relevance for perfume manufacture

BACKGROUND

Sclareol is a diterpene natural product of high value for the fragrance industry. Its labdane carbon skeleton and its two hydroxyl groups also make it a valued starting material for semisynthesis of numerous commercial substances, including production of Ambrox® and related ambergris substitutes used in the formulation of high end perfumes. Most of the commercially-produced sclareol is derived from cultivated clary sage (Salvia sclarea) and extraction of the plant material. In clary sage, sclareol mainly accumulates in essential oil-producing trichomes that densely cover flower calices. Manool also is a minor diterpene of this species and the main diterpene of related Salvia species.

RESULTS

Based on previous general knowledge of diterpene biosynthesis in angiosperms, and based on mining of our recently published transcriptome database obtained by deep 454-sequencing of cDNA from clary sage calices, we cloned and functionally characterized two new diterpene synthase (diTPS) enzymes for the complete biosynthesis of sclareol in clary sage. A class II diTPS (SsLPPS) produced labda-13-en-8-ol diphosphate as major product from geranylgeranyl diphosphate (GGPP) with some minor quantities of its non-hydroxylated analogue, (9 S, 10 S)-copalyl diphosphate. A class I diTPS (SsSS) then transformed these intermediates into sclareol and manool, respectively. The production of sclareol was reconstructed in vitro by combining the two recombinant diTPS enzymes with the GGPP starting substrate and in vivo by co-expression of the two proteins in yeast (Saccharomyces cerevisiae). Tobacco-based transient expression assays of green fluorescent protein-fusion constructs revealed that both enzymes possess an N-terminal signal sequence that actively targets SsLPPS and SsSS to the chloroplast, a major site of GGPP and diterpene production in plants.

CONCLUSIONS

SsLPPS and SsSS are two monofunctional diTPSs which, together, produce the diterpenoid specialized metabolite sclareol in a two-step process. They represent two of the first characterized hydroxylating diTPSs in angiosperms and generate the dihydroxylated labdane sclareol without requirement for additional enzymatic oxidation by activities such as cytochrome P450 monoxygenases. Yeast-based production of sclareol by co-expresssion of SsLPPS and SsSS was efficient enough to warrant the development and use of such technology for the biotechnological production of scareol and other oxygenated diterpenes.

Bioinformatics opportunities for identification and study of medicinal plants

Plants have been used as a source of medicine since historic times and several commercially important drugs are of plant-based origin. The traditional approach towards discovery of plant-based drugs often times involves significant amount of time and expenditure. These labor-intensive approaches have struggled to keep pace with the rapid development of high-throughput technologies. In the era of high volume, high-throughput data generation across the biosciences, bioinformatics plays a crucial role. This has generally been the case in the context of drug designing and discovery. However, there has been limited attention to date to the potential application of bioinformatics approaches that can leverage plant-based knowledge. Here, we review bioinformatics studies that have contributed to medicinal plants research. In particular, we highlight areas in medicinal plant research where the application of bioinformatics methodologies may result in quicker and potentially cost-effective leads toward finding plant-based remedies.

Amphilectane diterpenes from Salvia sclarea: biosynthetic considerations

Salviatrienes A and B, two new diterpenes belonging to the amphilectane/elisabethane family, have been isolated from an extract of clary sage (Salvia sclarea). These molecules are the first representatives of this family to be described from the plant kingdom. This study has led to consideration of the possible enzymatic machinery and biosynthesis pathways within S. sclarea.

Comprehensive functional analyses of expressed sequence tags in common wheat (Triticum aestivum)

About 1 million expressed sequence tag (EST) sequences comprising 125.3 Mb nucleotides were accreted from 51 cDNA libraries constructed from a variety of tissues and organs under a range of conditions, including abiotic stresses and pathogen challenges in common wheat (Triticum aestivum). Expressed sequence tags were assembled with stringent parameters after processing with inbuild scripts, resulting in 37,138 contigs and 215,199 singlets. In the assembled sequences, 10.6% presented no matches with existing sequences in public databases. Functional characterization of wheat unigenes by gene ontology annotation, mining transcription factors, full-length cDNA, and miRNA targeting sites were carried out. A bioinformatics strategy was developed to discover single-nucleotide polymorphisms (SNPs) within our large EST resource and reported the SNPs between and within (homoeologous) cultivars. Digital gene expression was performed to find the tissue-specific gene expression, and correspondence analysis was executed to identify common and specific gene expression by selecting four biotic stress-related libraries. The assembly and associated information cater a framework for future investigation in functional genomics.

Deep sequencing of small RNAs in plants: applied bioinformatics

Small RNAs, including microRNA and short-interfering RNAs, play important roles in plants. In recent years, developments in sequencing technology have enabled the large-scale discovery of sRNAs in various cells, tissues and developmental stages and in response to various stresses. This review describes the bioinformatics challenges to analysing these large datasets of short-RNA sequences and some of the solutions to those challenges.

Identification of lignin genes and regulatory sequences involved in secondary cell wall formation in Acacia auriculiformis and Acacia mangium via de novo transcriptome sequencing

Background

Acacia auriculiformis × Acacia mangium hybrids are commercially important trees for the timber and pulp industry in Southeast Asia. Increasing pulp yield while reducing pulping costs are major objectives of tree breeding programs. The general monolignol biosynthesis and secondary cell wall formation pathways are well-characterized but genes in these pathways are poorly characterized in Acacia hybrids. RNA-seq on short-read platforms is a rapid approach for obtaining comprehensive transcriptomic data and to discover informative sequence variants.

Results

We sequenced transcriptomes of A. auriculiformis and A. mangium from non-normalized cDNA libraries synthesized from pooled young stem and inner bark tissues using paired-end libraries and a single lane of an Illumina GAII machine. De novo assembly produced a total of 42,217 and 35,759 contigs with an average length of 496 bp and 498 bp for A. auriculiformis and A. mangium respectively. The assemblies of A. auriculiformis and A. mangium had a total length of 21,022,649 bp and 17,838,260 bp, respectively, with the largest contig 15,262 bp long. We detected all ten monolignol biosynthetic genes using Blastx and further analysis revealed 18 lignin isoforms for each species. We also identified five contigs homologous to R2R3-MYB proteins in other plant species that are involved in transcriptional regulation of secondary cell wall formation and lignin deposition. We searched the contigs against public microRNA database and predicted the stem-loop structures of six highly conserved microRNA families (miR319, miR396, miR160, miR172, miR162 and miR168) and one legume-specific family (miR2086). Three microRNA target genes were predicted to be involved in wood formation and flavonoid biosynthesis. By using the assemblies as a reference, we discovered 16,648 and 9,335 high quality putative Single Nucleotide Polymorphisms (SNPs) in the transcriptomes of A. auriculiformis and A. mangium, respectively, thus yielding useful markers for population genetics studies and marker-assisted selection.

Conclusion

We have produced the first comprehensive transcriptome-wide analysis in A. auriculiformis and A. mangium using de novo assembly techniques. Our high quality and comprehensive assemblies allowed the identification of many genes in the lignin biosynthesis and secondary cell wall formation in Acacia hybrids. Our results demonstrated that Next Generation Sequencing is a cost-effective method for gene discovery, identification of regulatory sequences, and informative markers in a non-model plant.