Conservation and diversity of microRNA-associated copper-regulatory networks in Populus trichocarpa

Genome-Wide Identification and Characterization of miRNAs and Natural Antisense Transcripts Show the Complexity of Gene Regulatory Networks for Secondary Metabolism in Aristolochia contorta

Aristolochia contorta Bunge is an academically and medicinally important plant species. It belongs to the magnoliids, with an uncertain phylogenetic position, and is one of the few plant species lacking a whole-genome duplication (WGD) event after the angiosperm-wide WGD. A. contorta has been an important traditional Chinese medicine material. Since it contains aristolochic acids (AAs), chemical compounds with nephrotoxity and carcinogenicity, the utilization of this plant has attracted widespread attention. Great efforts are being made to increase its bioactive compounds and reduce or completely remove toxic compounds. MicroRNAs (miRNAs) and natural antisense transcripts (NATs) are two classes of regulators potentially involved in metabolism regulation. Here, we report the identification and characterization of 223 miRNAs and 363 miRNA targets. The identified miRNAs include 51 known miRNAs belonging to 20 families and 172 novel miRNAs belonging to 107 families. A negative correlation between the expression of miRNAs and their targets was observed. In addition, we identified 441 A. contorta NATs and 560 NAT-sense transcript (ST) pairs, of which 12 NATs were targets of 13 miRNAs, forming 18 miRNA-NAT-ST modules. Various miRNAs and NATs potentially regulated secondary metabolism through the modes of miRNA-target gene-enzyme genes, NAT-STs, and NAT-miRNA-target gene-enzyme genes, suggesting the complexity of gene regulatory networks in A. contorta. The results lay a solid foundation for further manipulating the production of its bioactive and toxic compounds.

microRNAs: Key Players in Plant Response to Metal Toxicity

Environmental metal pollution is a common problem threatening sustainable and safe crop production. Heavy metals (HMs) cause toxicity by targeting key molecules and life processes in plant cells. Plants counteract excess metals in the environment by enhancing defense responses, such as metal chelation, isolation to vacuoles, regulating metal intake through transporters, and strengthening antioxidant mechanisms. In recent years, microRNAs (miRNAs), as a small non-coding RNA, have become the central regulator of a variety of abiotic stresses, including HMs. With the introduction of the latest technologies such as next-generation sequencing (NGS), more and more miRNAs have been widely recognized in several plants due to their diverse roles. Metal-regulated miRNAs and their target genes are part of a complex regulatory network. Known miRNAs coordinate plant responses to metal stress through antioxidant functions, root growth, hormone signals, transcription factors (TF), and metal transporters. This article reviews the research progress of miRNAs in the stress response of plants to the accumulation of HMs, such as Cu, Cd, Hg, Cr, and Al, and the toxicity of heavy metal ions.

In Planta Cell Wall Engineering: From Mutants to Artificial Cell Walls

To mitigate the effects of global warming and to preserve the limited fossil fuel resources, an increased exploitation of plant-based materials and fuels is required, which would be one of the most important innovations related to sustainable development. Cell walls account for the majority of plant dry biomass and so is the target of such innovations. In this review, we discuss recent advances in in planta cell wall engineering through genetic manipulations, with a focus on wild-type-based and mutant-based approaches. The long history of using a wild-type-based approach has resulted in the development of many strategies for manipulating lignin, hemicellulose and pectin to decrease cell wall recalcitrance. In addition to enzyme-encoding genes, many transcription factor genes important for changing relevant cell wall characteristics have been identified. Although mutant-based cell wall engineering is relatively new, it has become feasible due to the rapid development of genome-editing technologies and systems biology-related research; we will soon enter an age of designed artificial wood production via complex genetic manipulations of many industrially important trees and crops.

Comparative analysis identifies micro-RNA associated with nutrient homeostasis, development and stress response in Arabidopsis thaliana upon high Zn and metal hyperaccumulator Arabidopsis halleri

miRNAs have been found to be key players in mineral homeostasis, both in the control of nutrient balance and in the response to toxic trace elements. However, the effect of Zn excess on miRNAs has not been elucidated; moreover, no data are present regarding miRNAs in hyperaccumulator species, where metal homeostasis is tightly regulated. Therefore, expression levels of mature miRNAs were measured by RNA-Seq in Zn-sensitive Arabidopsis thaliana grown in control conditions and upon high Zn, in soil and in Zn-hyperaccumulator Arabidopsis halleri grown in control conditions. Differential expression of notable miRNAs and their targets was confirmed by real-time RT-PCR. The comparison in A. thaliana revealed a small subset modulated upon Zn treatment that is associated with stress response and nutrient homeostasis. On the other hand, a more consistent group of miRNAs was differentially expressed in A. halleri compared with A. thaliana, reflecting inherent differences in nutritional requirements and response to stresses and plant growth and development. Overall, these results confirm the involvement of miRNAs in Zn homeostasis and support the hypothesis of distinct regulatory pathways in hyperaccumulator species.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

Plant miR397 and its functions

MicroRNAs (miRNAs) are noncoding, small RNAs of 20-24 nucleotides (nt) and function critically at the post-transcriptional level to regulate gene expression through cleaving mRNA targets or interfering with translation of the target mRNAs. They are broadly involved in many biological processes in plants. The miR397 family in plants contains several conserved members either in 21-nt or in 22-nt that mainly target the laccase (LAC) genes functioning in lignin synthesis and are involved in the development of plants under various conditions. Recent findings showed that miR397b in Arabidopsis could also target to Casein Kinase II Subunit Beta 3 (CKB3) and mediate circadian regulation and plant flowering. This review aims to summarise recent updates on miR397 and provides the available basis for understanding the functional mechanisms of miR397 in plant growth and development regulation and in response to external adverse stimulation.

Genome-Wide Identification of Copper Stress-Regulated and Novel MicroRNAs in Mulberry Leaf

Copper (Cu) is an essential trace element for plant growth and development. It is widely involved in respiration, photosynthesis, pollen formation, and other biological processes. Therefore, low or excessive copper causes damage to plants. Mulberry is an essential perennial economic tree. At present, research on the abiotic stress responses in mulberry is mainly focused on the identification of resistant germplasm resources and cloning of resistant genes. In contrast, studies on the resistance function of microRNAs and the regulatory gene responses to stress are rare. In this study, small RNA libraries (control and copper stressed) were constructed from mulberry leaf RNA. High-throughput sequencing and screening were employed, a total of 65 known miRNAs and 78 predicted novel mature miRNAs were identified, among which 40 miRNAs were differentially expressed under copper stress. Subsequently, expression patterns were verified for 14 miRNAs by real-time fluorescence quantitative PCR (qPCR). The target genes of miRNAs were validated by 5' RLM-RACE. Our results provide the bases for further study on the molecular mechanism of copper stress regulation in mulberry.

Knock-down of δ-aminolevulinic acid dehydratase via virus-induced gene silencing alters the microRNA biogenesis and causes stress-related reactions in citrus plants

The δ-aminolevulinic acid (δ-ALA) is an intermediate in the biosynthetic pathway of tetrapyrroles. Tetrapyrroles play vital roles in many biological processes such as photosynthesis, respiration, and light-sensing. ALA-dehydratase (ALAD) combines two molecules of δ-ALA to form porphobilinogen. In citrus, the silencing of ALAD caused discrete yellow spots and necrosis in leaves and stems. Additionally, it caused rapid death in developing new shoots. Herein, we hypothesize that the accumulation of δ-ALA results in severe stress and reduced meristem development. For that reason, we investigated the dynamic changes in the expression profiles of 23 microRNA (miRNA) identified through small RNA sequencing, from CTV-tALAD plants in comparison with healthy C. macrophylla and C. macrophylla infiltrated with CTV-wt. Furthermore, we reported the effect of ALAD silencing on the total phenolics, H₂O_2, and reactive oxygen species (ROS) levels, to examine the possibilities of miRNAs involving the regulation of these pathways. Our results showed that the total phenolics content, H₂O₂, and O₂^- levels were increased in CTV-tALAD plants. Moreover, 63 conserved miRNA members belonging to 23 different miRNA families were differentially expressed in CTV-tALAD plants compared to controls. The identified miRNAs are implicated in auxin biosynthesis and signaling, axillary shoot meristem formation and leaf morphology, starch metabolism, and oxidative stress. Collectively, our findings suggested that ALAD silencing initiates stress on citrus plants. As a result, CTV-tALAD plants exhibit reduced metabolic rate, growth, and development in order to cope with the stress that resulted from the accumulation of δ-ALA. This cascade of events led to leaf, stem, and meristem necrosis and failure of new shoot development.

The Critical Role of miRNAs in Regulation of Flowering Time and Flower Development

Flowering is an important biological process for plants that ensures reproductive success. The onset of flowering needs to be coordinated with an appropriate time of year, which requires tight control of gene expression acting in concert to form a regulatory network. MicroRNAs (miRNAs) are non-coding RNAs known as master modulators of gene expression at the post-transcriptional level. Many different miRNA families are involved in flowering-related processes such as the induction of floral competence, floral patterning, and the development of floral organs. This review highlights the diverse roles of miRNAs in controlling the flowering process and flower development, in combination with potential biotechnological applications for miRNAs implicated in flower regulation.

MicroRNA528, a hub regulator modulating ROS homeostasis via targeting of a diverse set of genes encoding copper-containing proteins in monocots

Plant microRNAs (miRNAs) regulate vital cellular processes, including responses to extreme temperatures with which reactive oxygen species (ROS) are often closely associated. In the present study, it was found that aberrant temperatures caused extensive changes in abundance to numerous miRNAs in banana fruit, especially the copper (Cu)-associated miRNAs. Among them, miR528 was significantly downregulated under cold stress and it was found to target genes encoding polyphenol oxidase (PPO), different from those identified in rice and maize. Expression of PPO genes was upregulated by > 100-fold in cold conditions, leading to ROS surge and subsequent peel browning of banana fruit. Extensive comparative genomic analyses revealed that the monocot-specific miR528 can potentially target a large collection of genes encoding Cu-containing proteins. Most of them are actively involved in cellular ROS metabolism, including not only ROS generating oxidases, but also ROS scavenging enzymes. It also was demonstrated that miR528 has evolved a distinct preference of target genes in different monocots, with its target site varying in position among/within gene families, implying a highly dynamic process of target gene diversification. Its broad capacity to target genes encoding Cu-containing protein implicates miR528 as a key regulator for modulating the cellular ROS homeostasis in monocots.

De novo origination of MIRNAs through generation of short inverted repeats in target genes

MIRNA (MIR) gene origin and early evolutionary processes, such as hairpin precursor sequence origination, promoter activity acquirement and the sequence of these two processes, are fundamental and fascinating subjects. Three models, including inverted gene duplication, spontaneous evolution and transposon transposition, have been proposed for de novo origination of hairpin precursor sequence. However, these models still open to discussion. In addition, de novo origination of MIR gene promoters has not been well investigated. Here, I systematically investigated the origin of evolutionarily young polyphenol oxidase gene (PPO)-targeting MIRs, including MIR1444, MIR058 and MIR12112, and a genomic region termed AasPPO-as-hp, which contained a hairpin-forming sequence. I found that MIR058 precursors and the hairpin-forming sequence of AasPPO-as-hp originated in an ancient PPO gene through forming short inverted repeats. Palindromic-like sequences and imperfect inverted repeats in the ancient PPO gene contributed to initiate the generation of short inverted repeats probably by causing errors during DNA duplication. Analysis of MIR058 and AasPPO-as-hp promoters showed that they originated in the 3'-flanking region of the ancient PPO gene. Promoter activities were gained by insertion of a CAAT-box and multiple-copper-response element (CuRE)-containing miniature inverted-repeat transposable element (MITE) in the upstream of AT-rich TATA-box-like sequence. Gain of promoter activities occurred before hairpin-forming sequence origination. Sequence comparison of MIR1444, MIR058 and MIR12112 promoters showed frequent birth and death of CuREs, indicating copper could be vital for the origination and evolution of PPO-targeting MIRs. Based on the evidence obtained, a novel model for plant MIR origination and evolution is proposed.

Growth-regulating factor 15 is required for leaf size control in Populus

Growth-regulating factors (GRFs) are involved in various developmental events, particularly leaf development. However, the functions of GRFs in woody plants remain elusive. In this study, functional characterization of GRF15 in Populus was performed. Most GRFs are preferentially expressed in young leaves. As GRF15 was expressed at the highest level and with highest ratio in Populus species with large leaves, this gene was investigated through transgenic analyses. Promotor-β-glucuronidase analysis revealed expression of GRF15 at the leaf expansion zone. Additionally, GRF15 was found to be localized in the nucleus and regulated by miR396. Leaf size and palisade cell size were significantly increased and decreased in GRF15-overexpressing and dominant repression lines, respectively. Consistently, expression of EXPA11a, a homolog of cell-expansion marker EXPA11 in Arabidopsis, was strongly upregulated and downregulated in the GRF15-overexpressing and dominant repression lines, respectively, which was further manifested by activation of EXPA11a by GRF15 in transactivation assays. Therefore, GRF15 is required for leaf size control and primarily modulates cell expansion during leaf development in Populus.

Identification of copper (Cu) stress-responsive grapevine microRNAs and their target genes by high-throughput sequencing

MicroRNAs (miRNAs) are a class of single-stranded non-coding small RNAs (sRNAs) that are 20-24 nucleotides (nt) in length. Extensive studies have indicated that miRNAs play important roles in plant growth, development and stress responses. With more copper (Cu) and copper containing compounds used as bactericides and fungicides in plants, Cu stress has become one of the serious environmental problems that affect plant growth and development. In order to uncover the hidden response mechanisms of Cu stress, two small RNA libraries were constructed from Cu-treated and water-treated (Control) leaves of 'Summer Black' grapevine. Following high-throughput sequencing and filtering, a total of 158 known and 98 putative novel miRNAs were identified in the two libraries. Among these, 100 known and 47 novel miRNAs were identified as differentially expressed under Cu stress. Subsequently, the expression patterns of nine Cu-responsive miRNAs were validated by quantitative real-time PCR (qRT-PCR). There existed some consistency in expression levels of Cu-responsive miRNAs between high throughput sequencing and qRT-PCR assays. The targets prediction of miRNAs indicates that miRNA may regulate some transcription factors, including AP2, SBP, NAC, MYB and ARF during Cu stress. The target genes for two known and two novel miRNAs showed specific cleavage sites at the 10th and/or 11th nucleotide from the 5'-end of the miRNA corresponding to their miRNA complementary sequences. The findings will lay the foundation for exploring the role of the regulation of miRNAs in response to Cu stress and provide valuable gene information for breeding some Cu-tolerant grapevine cultivars.

Overexpression of microRNA408 enhances photosynthesis, growth, and seed yield in diverse plants

The ability of a plant to produce grain, fruit, or forage depends ultimately on photosynthesis. There have been few attempts, however, to study microRNAs, which are a class of endogenous small RNAs post-transcriptionally programming gene expression, in relation to photosynthetic traits. We focused on miR408, one of the most conserved plant miRNAs, and overexpressed it in parallel in Arabidopsis, tobacco, and rice. The transgenic plants all exhibited increased copper content in the chloroplast, elevated abundance of plastocyanin, and an induction of photosynthetic genes. By means of gas exchange and optical spectroscopy analyses, we showed that higher expression of miR408 leads to enhanced photosynthesis through improving efficiency of irradiation utilization and the capacity for carbon dioxide fixation. Consequently, miR408 hyper-accumulating plants exhibited higher rate of vegetative growth. An enlargement of seed size was also observed in all three species overproducing miR408. Moreover, we conducted a 2-year-two-location field trial and observed miR408 overexpression in rice significantly increased yield, which was primarily attributed to an elevation in grain weight. Taken together, these results demonstrate that miR408 is a positive regulator of photosynthesis and that its genetic engineering is a promising route for enhancing photosynthetic performance and yield in diverse plants.

Characterization of Conserved and Novel microRNAs in Lilium lancifolium Thunb. by High-Throughput Sequencing

MicroRNAs (miRNAs) are among the class of noncoding small RNA molecules and play a crucial role in post-transcriptional regulation in plants. Although Lilium is one of the most popular ornamental flowers worldwide, however, there is no report on miRNAs identification. In the present study, therefore, miRNAs and their targets were identified from flower, leaf, bulblet and bulb of Lilium lancifolium Thunb. by high-throughput sequencing and bioinformatics analysis. In this study, a total of 38 conserved miRNAs belonging to 17 miRNA families and 44 novel miRNAs were identified. In total, 366 target genes for conserved miRNAs and 415 target genes for novel miRNAs were predicted. The majority of the target genes for conserved miRNAs were transcriptional factors and novel miRNAs targeted mainly protein coding genes. A total of 53 cleavage sites belonging to 6 conserved miRNAs families and 14 novel miRNAs were identified using degradome sequencing. Twenty-three miRNAs were randomly selected, then, their credibility was confirmed using northern blot or stem-loop qRT-PCR. The results from qRT-PCR analysis showed the expression pattern of 4 LL-miRNAs was opposite to their targets. Therefore, our finding provides an important basis to understand the biological functions of miRNAs in Lilium.

Roles of Non-Coding RNA in Sugarcane-Microbe Interaction

Studies have highlighted the importance of non-coding RNA regulation in plant-microbe interaction. However, the roles of sugarcane microRNAs (miRNAs) in the regulation of disease responses have not been investigated. Firstly, we screened the sRNA transcriptome of sugarcane infected with Acidovorax avenae. Conserved and novel miRNAs were identified. Additionally, small interfering RNAs (siRNAs) were aligned to differentially expressed sequences from the sugarcane transcriptome. Interestingly, many siRNAs aligned to a transcript encoding a copper-transporter gene whose expression was induced in the presence of A. avenae, while the siRNAs were repressed in the presence of A. avenae. Moreover, a long intergenic non-coding RNA was identified as a potential target or decoy of miR408. To extend the bioinformatics analysis, we carried out independent inoculations and the expression patterns of six miRNAs were validated by quantitative reverse transcription-PCR (qRT-PCR). Among these miRNAs, miR408-a copper-microRNA-was downregulated. The cleavage of a putative miR408 target, a laccase, was confirmed by a modified 5'RACE (rapid amplification of cDNA ends) assay. MiR408 was also downregulated in samples infected with other pathogens, but it was upregulated in the presence of a beneficial diazotrophic bacteria. Our results suggest that regulation by miR408 is important in sugarcane sensing whether microorganisms are either pathogenic or beneficial, triggering specific miRNA-mediated regulatory mechanisms accordingly.

Spatial regulation of monolignol biosynthesis and laccase genes control developmental and stress-related lignin in flax

BACKGROUND

Bast fibres are characterized by very thick secondary cell walls containing high amounts of cellulose and low lignin contents in contrast to the heavily lignified cell walls typically found in the xylem tissues. To improve the quality of the fiber-based products in the future, a thorough understanding of the main cell wall polymer biosynthetic pathways is required. In this study we have carried out a characterization of the genes involved in lignin biosynthesis in flax along with some of their regulation mechanisms.

RESULTS

We have first identified the members of the phenylpropanoid gene families through a combination of in silico approaches. The more specific lignin genes were further characterized by high throughput transcriptomic approaches in different organs and physiological conditions and their cell/tissue expression was localized in the stems, roots and leaves. Laccases play an important role in the polymerization of monolignols. This multigenic family was determined and a miRNA was identified to play a role in the posttranscriptional regulation by cleaving the transcripts of some specific genes shown to be expressed in lignified tissues. In situ hybridization also showed that the miRNA precursor was expressed in the young xylem cells located near the vascular cambium. The results obtained in this work also allowed us to determine that most of the genes involved in lignin biosynthesis are included in a unique co-expression cluster and that MYB transcription factors are potentially good candidates for regulating these genes.

CONCLUSIONS

Target engineering of cell walls to improve plant product quality requires good knowledge of the genes responsible for the production of the main polymers. For bast fiber plants such as flax, it is important to target the correct genes from the beginning since the difficulty to produce transgenic material does not make possible to test a large number of genes. Our work determined which of these genes could be potentially modified and showed that it was possible to target different regulatory pathways to modify lignification.

Genome-wide identification and expression analysis of SBP-like transcription factor genes in Moso Bamboo (Phyllostachys edulis)

BACKGROUND

The SQUAMOSA promoter binding protein-like (SPL) proteins are plant-specific transcription factors (TFs) that function in a variety of developmental processes including growth, flower development, and signal transduction. SPL proteins are encoded by a gene family, and these genes have been characterized in two model grass species, Zea mays and Oryza sativa. The SPL gene family has not been well studied in moso bamboo (Phyllostachys edulis), a woody grass species.

RESULTS

We identified 32 putative PeSPL genes in the P. edulis genome. Phylogenetic analysis arranged the PeSPL protein sequences in eight groups. Similarly, phylogenetic analysis of the SBP-like and SBP proteins from rice and maize clustered them into eight groups analogous to those from P. edulis. Furthermore, the deduced PeSPL proteins in each group contained very similar conserved sequence motifs. Our analyses indicate that the PeSPL genes experienced a large-scale duplication event ~15 million years ago (MYA), and that divergence between the PeSPL and OsSPL genes occurred 34 MYA. The stress-response expression profiles and tissue-specificity of the putative PeSPL gene promoter regions showed that SPL genes in moso bamboo have potential biological functions in stress resistance as well as in growth and development. We therefore examined PeSPL gene expression in response to different plant hormone and drought (polyethylene glycol-6000; PEG) treatments to mimic biotic and abiotic stresses. Expression of three (PeSPL10, -12, -17), six (PeSPL1, -10, -12, -17, -20, -31), and nine (PeSPL5, -8, -9, -14, -15, -19, -20, -31, -32) genes remained relatively stable after treating with salicylic acid (SA), gibberellic acid (GA), and PEG, respectively, while the expression patterns of other genes changed. In addition, analysis of tissue-specific expression of the moso bamboo SPL genes during development showed differences in their spatiotemporal expression patterns, and many were expressed at high levels in flowers and leaves.

CONCLUSIONS

The PeSPL genes play important roles in plant growth and development, including responses to stresses, and most of the genes are expressed in different tissues. Our study provides a comprehensive understanding of the PeSPL gene family and may enable future studies on the function and evolution of SPL genes in moso bamboo.

Characterization of the polyphenol oxidase gene family reveals a novel microRNA involved in posttranscriptional regulation of PPOs in Salvia miltiorrhiza

Salvia miltiorrhiza is a well-known material of traditional Chinese medicine. Understanding the regulatory mechanisms of phenolic acid biosynthesis and metabolism are important for S. miltiorrhiza quality improvement. We report here that S. miltiorrhiza contains 19 polyphenol oxidases (PPOs), forming the largest PPO gene family in plant species to our knowledge. Analysis of gene structures and sequence features revealed the conservation and divergence of SmPPOs. SmPPOs were differentially expressed in plant tissues and eight of them were predominantly expressed in phloem and xylem, indicating that some SmPPOs are functionally redundant, whereas the others are associated with different physiological processes. Expression patterns of eighteen SmPPOs were significantly altered under MeJA treatment, and twelve were yeast extract and Ag⁺-responsive, suggesting the majority of SmPPOs are stress-responsive. Analysis of high-throughput small RNA sequences and degradome data showed that miR1444-mediated regulation of PPOs existing in P. trichocarpa is absent from S. miltiorrhiza. Instead, a subset of SmPPOs was posttranscriptionally regulated by a novel miRNA, termed Smi-miR12112. It indicates the specificity and significance of miRNA-mediated regulation of PPOs. The results shed light on the regulation of SmPPO expression and suggest the complexity of SmPPO-associated phenolic acid biosynthesis and metabolism.

The copper microRNAs

1030 I. 1030 II. 1030 III. 1031 IV. 1031 V. 1032 VI. 1033 VII. 1034 VIII. 1034 1034 References 1034 SUMMARY: Copper (Cu) microRNAs are upregulated by Cu deficiency and mediate the post-transcriptional downregulation of transcripts that encode Cu proteins, suggesting a role directly related to Cu. However, expression and phenotypic analyses of copper microRNA mutants and over-expressors have suggested roles mainly in tolerance to abiotic stresses. To reconcile available data, a model is proposed which emphasizes the mobile nature of copper microRNA molecules in the regulation of Cu homeostasis. It is proposed that the Cu-microRNA regulatory circuits are further co-opted by plants to regulate both beneficial and pathogenic interactions with microbes. Further exploration of Cu-microRNA functions that account for the cell-to-cell mobility should give novel insight into plant microbe interactions and the integration of micronutrition and development.

Origin and evolution of MIR1444 genes in Salicaceae

miR1444s are functionally significant miRNAs targeting polyphenol oxidase (PPO) genes for cleavage. MIR1444 genes were reported only in Populus trichocarpa. Through the computational analysis of 215 RNA-seq data, four whole genome sequences of Salicaceae species and deep sequencing of six P. trichocarpa small RNA libraries, we investigated the origin and evolution history of MIR1444s. A total of 23 MIR1444s were identified. Populus and Idesia species contain two MIR1444 genes, while Salix includes only one. Populus and Idesia MIR1444b genes and Salix MIR1444s were phylogenetically separated from Populus and Idesia MIR1444a genes. Ptr-miR1444a and ptr-miR1444b showed sequence divergence. Compared with ptr-miR1444b, ptr-miR1444a started 2 nt upstream of precursor, resulting in differential regulation of PPO targets. Sequence alignments showed that MIR1444 genes exhibited extensive similarity to their PPO targets, the characteristics of MIRs originated from targets through an inverted gene duplication event. Genome sequence comparison showed that MIR1444 genes in Populus and Idesia were expanded through the Salicoid genome duplication event. A copy of MIR1444 gene was lost in Salix through DNA segment deletion during chromosome rearrangements. The results provide significant information for the origin of plant miRNAs and the mechanism of Salicaceae gene evolution and divergence.

High-throughput sequencing of small RNAs and anatomical characteristics associated with leaf development in celery

MicroRNAs (miRNAs) exhibit diverse and important roles in plant growth, development, and stress responses and regulate gene expression at the post-transcriptional level. Knowledge about the diversity of miRNAs and their roles in leaf development in celery remains unknown. To elucidate the roles of miRNAs in celery leaf development, we identified leaf development-related miRNAs through high-throughput sequencing. Small RNA libraries were constructed using leaves from three stages (10, 20, and 30 cm) of celery cv.'Ventura' and then subjected to high-throughput sequencing and bioinformatics analysis. At Stage 1, Stage 2, and Stage 3 of 'Ventura', a total of 333, 329, and 344 conserved miRNAs (belonging to 35, 35, and 32 families, respectively) were identified. A total of 131 miRNAs were identified as novel in 'Ventura'. Potential miRNA target genes were predicted and annotated using the eggNOG, GO, and KEGG databases to explore gene functions. The abundance of five conserved miRNAs and their corresponding potential target genes were validated. Expression profiles of novel potential miRNAs were also detected. Anatomical characteristics of the leaf blades and petioles at three leaf stages were further analyzed. This study contributes to our understanding on the functions and molecular regulatory mechanisms of miRNAs in celery leaf development.

To bloom or not to bloom: role of microRNAs in plant flowering

During the course of their life cycles, plants undergo various morphological and physiological changes underlying juvenile-to-adult and adult-to-flowering phase transitions. To flower or not to flower is a key step of plasticity of a plant toward the start of its new life cycle. In addition to the previously revealed intrinsic genetic programs, exogenous cues, and endogenous cues, a class of small non-coding RNAs, microRNAs (miRNAs), plays a key role in plants making the decision to flower by integrating into the known flowering pathways. This review highlights the age-dependent flowering pathway with a focus on a number of timing miRNAs in determining such a key process. The contributions of other miRNAs which exist mainly outside the age pathway are also discussed. Approaches to study the flowering-determining miRNAs, their interactions, and applications are presented.

Prioritization of copper for the use in photosynthetic electron transport in developing leaves of hybrid poplar

Plastocyanin (PC) is an essential and abundant copper (Cu) protein required for photosynthesis in higher plants. Severe copper deprivation has the potential to cause a defect in photosynthetic electron transport due to a lack in PC. The Cu-microRNAs, which are up-regulated under Cu deficiency, down-regulate the expression of target Cu proteins other than PC, cytochrome-c oxidase and the ethylene receptors. It has been proposed that this mechanism saves Cu for PC maturation. We aimed to test how hybrid poplar, a species that has capacity to rapidly expand its photosynthetically active tissue, responds to variations in Cu availability over time. Measurement of chlorophyll fluorescence after Cu depletion revealed a drastic effect on photosynthesis in hybrid poplar. The decrease in photosynthetic capacity was correlated with a reduction in PC protein levels. Compared to older leaves, PC decreased more strongly in developing leaves, which also lost more photosynthetic electron transport capacity. The effect of Cu depletion on older and more developed leaves was minor and these leaves maintained much of their photosynthetic capacity. Interestingly, upon resupply of Cu to the medium a very rapid recovery of Cu levels was seen in the younger leaves with a concomitant rise in the expression and activity of PC. In contrast, the expression of those Cu proteins, which are targets of microRNAs was under the same circumstances delayed. At the same time, Cu resupply had only minor effects on the older leaves. The data suggest a model where rapid recovery of photosynthetic capacity in younger leaves is made possible by a preferred allocation of Cu to PC in younger leaves, which is supported by Cu-microRNA expression.

Novel and conserved heat-responsive microRNAs in wheat (Triticum aestivum L.)

MicroRNAs (miRNAs) are small endogenous RNAs of ~22 nucleotides that have been shown to play regulatory role by negatively affecting the expression of genes at the post-transcriptional level. Information of miRNAs on some important crops like soybean, Arabidopsis, and rice, etc. are available, but no study on heat-responsive novel miRNAs has yet been reported in wheat (Triticum aestivum L.). In the present investigation, a popular wheat cultivar HD2985 was used in small RNA library construction and Illumina HiSeq 2000 was used to perform high-throughput sequencing of the library after cluster generation; 110,896,604 and 87,743,861 reads were generated in the control (22 °C) and heat-treated (42 °C for 2 h) samples, respectively. Forty-four precursor and mature miRNAs were found in T. aestivum from miRBase v 19. The frequencies of the miRNA families varied from 2 (tae-miR1117) to 60,672 (tae-miR159b). We identify 1052 and 902 mature miRNA sequences in HD2985 control and HS-treated samples by mapping on reference draft genome of T. aestivum. Maximum identified miRNAs were located on IWGSC_CSS_3B_scaff (chromosome 3B). We could identify 53 and 46 mature miRNA in the control and HS samples and more than 516 target genes by mapping on the reference genome of Oryza sativa, Zea mays, and Sorghum bicolor. Using different pipelines and plant-specific criteria, 37 novel miRNAs were identified in the control and treated samples. Six novel miRNA were validated using qRT-PCR to be heat-responsive. A negative correlation was, however, observed between the expression of novel miRNAs and their targets. Target prediction and pathway analysis revealed their involvement in the heat stress tolerance. These novel miRNAs are new additions to miRNA database of wheat, and the regulatory network will be made use of in deciphering the mechanism of thermotolerance in wheat.

Molecular characterization of the SPL gene family in Populus trichocarpa

BACKGROUND

SPLs, a family of transcription factors specific to plants, play vital roles in plant growth and development through regulation of various physiological and biochemical processes. Although Populus trichocarpa is a model forest tree, the PtSPL gene family has not been systematically studied.

RESULTS

Here we report the identification of 28 full-length PtSPLs, which distribute on 14 P. trichocarpa chromosomes. Based on the phylogenetic relationships of SPLs in P. trichocarpa and Arabidopsis, plant SPLs can be classified into 6 groups. Each group contains at least a PtSPL and an AtSPL. The N-terminal zinc finger 1 (Zn1) of SBP domain in group 6 SPLs has four cysteine residues (CCCC-type), while Zn1 of SPLs in the other groups mainly contains three cysteine and one histidine residues (C2HC-type). Comparative analyses of gene structures, conserved motifs and expression patterns of PtSPLs and AtSPLs revealed the conservation of plant SPLs within a group, whereas among groups, the P. trichocarpa and Arabidopsis SPLs were significantly different. Various conserved motifs were identified in PtSPLs but not found in AtSPLs, suggesting the diversity of plant SPLs. A total of 11 pairs of intrachromosome-duplicated PtSPLs were identified, suggesting the importance of gene duplication in SPL gene expansion in P. trichocarpa. In addition, 18 of the 28 PtSPLs, belonging to G1, G2 and G5, were found to be targets of miR156. Consistently, all of the AtSPLs in these groups are regulated by miR156. It suggests the conservation of miR156-mediated posttranscriptional regulation in plants.

CONCLUSIONS

A total of 28 full-length SPLs were identified from the whole genome sequence of P. trichocarpa. Through comprehensive analyses of gene structures, phylogenetic relationships, chromosomal locations, conserved motifs, expression patterns and miR156-mediated posttranscriptional regulation, the PtSPL gene family was characterized. Our results provide useful information for evolution and biological function of plant SPLs.

Genome-wide characterization and comparative analysis of R2R3-MYB transcription factors shows the complexity of MYB-associated regulatory networks in Salvia miltiorrhiza

Background

MYB is the largest plant transcription factor gene family playing vital roles in plant growth and development. However, it has not been systematically studied in Salvia miltiorrhiza, an economically important medicinal plant.

Results

Here we report the genome-wide identification and characterization of 110 R2R3-MYBs, the largest subfamily of MYBs in S. miltiorrhiza. The MYB domain and other motifs of SmMYBs are largely conserved with Arabidopsis AtMYBs, whereas the divergence of SmMYBs and AtMYBs also exists, suggesting the conservation and diversity of plant MYBs. SmMYBs and AtMYBs may be classified into 37 subgroups, of which 31 include proteins from S. miltiorrhiza and Arabidopsis, whereas 6 are specific to a species, indicating that the majority of MYBs play conserved roles, while others may exhibit species-specialized functions. SmMYBs are differentially expressed in various tissues of S. miltiorrhiza. The expression profiles are largely consistent with known functions of their Arabidopsis counterparts. The expression of a subset of SmMYBs is regulated by microRNAs, such as miR159, miR319, miR828 and miR858. Based on functional conservation of MYBs in a subgroup, SmMYBs potentially involved in the biosynthesis of bioactive compounds were identified.

Conclusions

A total of 110 R2R3-MYBs were identified and analyzed. The results suggest the complexity of MYB-mediated regulatory networks in S. miltiorrhiza and provide a foundation for understanding the regulatory mechanism of SmMYBs.

Identification and characterization of a subset of microRNAs in wheat (Triticum aestivum L.)

MicroRNAs (miRNAs) represent a class of endogenous regulator for post-transcriptionally modulating gene expression. Elucidating complete miRNA repertoires for individual species is a long-desired goal in miRNA research. So far only 42 have been annotated for common wheat (Triticum aestivum) due to its large genome. Here, we employed miRDeep-P, a program developed previously for retrieving miRNAs from deep sequencing data in plants, to parse 14 sequenced small RNA libraries of wheat using expression sequence tags (ESTs) as the reference in lieu of a complete genome sequence. This effort identified 145 miRNAs along with the corresponding stem-looped precursors with many differentially expressed in development and associated with powdery mildew infection. Examination of the phylogenetic distribution of these miRNAs and their target genes revealed that many are conserved in monocots. The set of miRNAs identified in this study is thus useful to further miRNA research in wheat and other important cereal crop species.

Genome-wide analysis and molecular dissection of the SPL gene family in Salvia miltiorrhiza

SQUAMOSA promoter binding protein-likes (SPLs) are plant-specific transcription factors playing vital regulatory roles in plant growth and development. There is no information about SPLs in Salvia miltiorrhiza (Danshen), a significant medicinal plant widely used in Traditional Chinese medicine (TCM) for >1,700 years and an emerging model plant for TCM studies. Through genome-wide identification and subsequent molecular cloning, we identified a total 15 SmSPLs with divergent sequence features, gene structures, and motifs. Comparative analysis showed sequence conservation between SmSPLs and their Arabidopsis counterparts. A phylogenetic tree clusters SmSPLs into six groups. Many of the motifs identified commonly exist in a group/subgroup, implying their functional redundancy. Eight SmSPLs were predicted and experimentally validated to be targets of miR156/157. SmSPLs were differentially expressed in various tissues of S. milltiorrhiza. The expression of miR156/157-targeted SmSPLs was increased with the maturation of S. miltiorrhiza, whereas the expression of miR156/157 was decreased, confirming the regulatory roles of miR156/157 in SmSPLs and suggesting the functions of SmSPLs in S. miltiorrhiza development. The expression of miR156/157 was negatively correlated with miR172 during the maturation of S. miltiorrhiza. The results indicate the significance and complexity of SmSPL-, miR156-, and miR172-mediated regulation of developmental timing in S. miltiorrhiza.

Copper and iron homeostasis in plants: the challenges of oxidative stress

SIGNIFICANCE

Photosynthesis, the process that drives life on earth, relies on transition metal (e.g., Fe and Cu) containing proteins that participate in electron transfer in the chloroplast. However, the light reactions also generate high levels of reactive oxygen species (ROS), which makes metal use in plants a challenge.

RECENT ADVANCES

Sophisticated regulatory networks govern Fe and Cu homeostasis in response to metal ion availability according to cellular needs and priorities. Molecular remodeling in response to Fe or Cu limitation leads to its economy to benefit photosynthesis. Fe toxicity is prevented by ferritin, a chloroplastic Fe-storage protein in plants. Recent studies on ferritin function and regulation revealed the interplay between iron homeostasis and the redox balance in the chloroplast.

CRITICAL ISSUES

Although the connections between metal excess and ROS in the chloroplast are established at the molecular level, the mechanistic details and physiological significance remain to be defined. The causality/effect relationship between transition metals, redox signals, and responses is difficult to establish.

FUTURE DIRECTIONS

Integrated approaches have led to a comprehensive understanding of Cu homeostasis in plants. However, the biological functions of several major families of Cu proteins remain unclear. The cellular priorities for Fe use under deficiency remain largely to be determined. A number of transcription factors that function to regulate Cu and Fe homeostasis under deficiency have been characterized, but we have not identified regulators that mediate responses to excess. Importantly, details of metal sensing mechanisms and cross talk to ROS-sensing mechanisms are so far poorly documented in plants.

Genome-wide identification and analysis of the SBP-box family genes in apple (Malus × domestica Borkh.)

SQUAMOSA promoter binding protein (SBP)-box genes encode a family of plant-specific transcription factors and play many crucial roles in plant development. In this study, 27 SBP-box gene family members were identified in the apple (Malus × domestica Borkh.) genome, 15 of which were suggested to be putative targets of MdmiR156. Plant SBPs were classified into eight groups according to the phylogenetic analysis of SBP-domain proteins. Gene structure, gene chromosomal location and synteny analyses of MdSBP genes within the apple genome demonstrated that tandem and segmental duplications, as well as whole genome duplications, have likely contributed to the expansion and evolution of the SBP-box gene family in apple. Additionally, synteny analysis between apple and Arabidopsis indicated that several paired homologs of MdSBP and AtSPL genes were located in syntenic genomic regions. Tissue-specific expression analysis of MdSBP genes in apple demonstrated their diversified spatiotemporal expression patterns. Most MdmiR156-targeted MdSBP genes, which had relatively high transcript levels in stems, leaves, apical buds and some floral organs, exhibited a more differential expression pattern than most MdmiR156-nontargeted MdSBP genes. Finally, expression analysis of MdSBP genes in leaves upon various plant hormone treatments showed that many MdSBP genes were responsive to different plant hormones, indicating that MdSBP genes may be involved in responses to hormone signaling during stress or in apple development.

Ptr-miR397a is a negative regulator of laccase genes affecting lignin content in Populus trichocarpa

Laccases, as early as 1959, were proposed to catalyze the oxidative polymerization of monolignols. Genetic evidence in support of this hypothesis has been elusive due to functional redundancy of laccase genes. An Arabidopsis double mutant demonstrated the involvement of laccases in lignin biosynthesis. We previously identified a subset of laccase genes to be targets of a microRNA (miRNA) ptr-miR397a in Populus trichocarpa. To elucidate the roles of ptr-miR397a and its targets, we characterized the laccase gene family and identified 49 laccase gene models, of which 29 were predicted to be targets of ptr-miR397a. We overexpressed Ptr-MIR397a in transgenic P. trichocarpa. In each of all nine transgenic lines tested, 17 PtrLACs were down-regulated as analyzed by RNA-seq. Transgenic lines with severe reduction in the expression of these laccase genes resulted in an ∼40% decrease in the total laccase activity. Overexpression of Ptr-MIR397a in these transgenic lines also reduced lignin content, whereas levels of all monolignol biosynthetic gene transcripts remained unchanged. A hierarchical genetic regulatory network (GRN) built by a bottom-up graphic Gaussian model algorithm provides additional support for a role of ptr-miR397a as a negative regulator of laccases for lignin biosynthesis. Full transcriptome-based differential gene expression in the overexpressed transgenics and protein domain analyses implicate previously unidentified transcription factors and their targets in an extended hierarchical GRN including ptr-miR397a and laccases that coregulate lignin biosynthesis in wood formation. Ptr-miR397a, laccases, and other regulatory components of this network may provide additional strategies for genetic manipulation of lignin content.

Identification of drought-responsive and novel Populus trichocarpa microRNAs by high-throughput sequencing and their targets using degradome analysis

BACKGROUND

MicroRNAs (miRNAs) are endogenous small RNAs (sRNAs) with a wide range of regulatory functions in plant development and stress responses. Although miRNAs associated with plant drought stress tolerance have been studied, the use of high-throughput sequencing can provide a much deeper understanding of miRNAs. Drought is a common stress that limits the growth of plants. To obtain more insight into the role of miRNAs in drought stress, Illumina sequencing of Populus trichocarpa sRNAs was implemented.

RESULTS

Two sRNA libraries were constructed by sequencing data of control and drought stress treatments of poplar leaves. In total, 207 P. trichocarpa conserved miRNAs were detected from the two sRNA libraries. In addition, 274 potential candidate miRNAs were found; among them, 65 candidates with star sequences were chosen as novel miRNAs. The expression of nine conserved miRNA and three novel miRNAs showed notable changes in response to drought stress. This was also confirmed by quantitative real time polymerase chain reaction experiments. To confirm the targets of miRNAs experimentally, two degradome libraries from the two treatments were constructed. According to degradome sequencing results, 53 and 19 genes were identified as targets of conserved and new miRNAs, respectively. Functional analysis of these miRNA targets indicated that they are involved in important activities such as the regulation of transcription factors, the stress response, and lipid metabolism.

CONCLUSIONS

We discovered five upregulated miRNAs and seven downregulated miRNAs in response to drought stress. A total of 72 related target genes were detected by degradome sequencing. These findings reveal important information about the regulation mechanism of miRNAs in P. trichocarpa and promote the understanding of miRNA functions during the drought response.

Evolution of variety-specific regulatory schema for expression of osa-miR408 in indica rice varieties under drought stress

Evolution of differential regulatory mechanisms can lead to quite distinct physiological attributes. In the present study, we have identified one such regulatory schema that regulates osa-miR408 and responds differentially in drought-sensitive and -tolerant indica rice varieties. A comparison of the drought stress response in drought-sensitive (Pusa Basmati 1 and IR64) and drought-tolerant (Nagina 22 and Vandana) indica rice varieties revealed that, during drought stress, levels of miR408 transcript decrease significantly in sensitive cultivars, whereas they remain elevated in the tolerant cultivars. The trend is reflected in young seedlings, as well as in flag leaf and spikelets of adult plants (heading stage). Members of the plastocyanin-like protein family targeted by miR408 also show the inverse expression profile and thus accumulate at a lower level in tolerant cultivars during drought. Interestingly, some members of this family are implicated in maintaining the cellular redox state and spikelet fertility in Arabidopsis. An investigation of miR408 loci (including promoter) in all four cultivars did not reveal any significant sequence variation indicating an involvement of the upstream regulatory schema. Indeed, a similar variety-specific stress response was found in the Oryza sativa squamosa promoter-binding-like 9 transcription factor that regulates miR408 expression. We further demonstrate that drought-mediated induction of miR408 in Nagina 22 is regulated by [Ca(2+)]cyt levels. However, [Ca(2+)]cyt does not appear to regulate miR408 levels in Pusa Basmati 1, suggesting a variety-specific evolution of regulatory schema in rice.

Two microRNAs linked to nodule infection and nitrogen-fixing ability in the legume Lotus japonicus

Legumes overcome nitrogen shortage by developing root nodules in which symbiotic bacteria fix atmospheric nitrogen in exchange for host-derived carbohydrates and mineral nutrients. Nodule development involves the distinct processes of nodule organogenesis, bacterial infection, and the onset of nitrogen fixation. These entail profound, dynamic gene expression changes, notably contributed to by microRNAs (miRNAs). Here, we used deep-sequencing, candidate-based expression studies and a selection of Lotus japonicus mutants uncoupling different symbiosis stages to identify miRNAs involved in symbiotic nitrogen fixation. Induction of a noncanonical miR171 isoform, which targets the key nodulation transcription factor Nodulation Signaling Pathway2, correlates with bacterial infection in nodules. A second candidate, miR397, is systemically induced in the presence of active, nitrogen-fixing nodules but not in that of noninfected or inactive nodule organs. It is involved in nitrogen fixation-related copper homeostasis and targets a member of the laccase copper protein family. These findings thus identify two miRNAs specifically responding to symbiotic infection and nodule function in legumes.

Expression-based functional investigation of the organ-specific microRNAs in Arabidopsis

MicroRNAs (miRNAs) play a pivotal role in plant development. The expression patterns of the miRNA genes significantly influence their regulatory activities. By utilizing small RNA (sRNA) high-throughput sequencing (HTS) data, the miRNA expression patterns were investigated in four organs (flowers, leaves, roots and seedlings) of Arabidopsis. Based on a set of criteria, dozens of organ-specific miRNAs were discovered. A dominant portion of the organ-specific miRNAs identified from the ARGONAUTE 4-enriched sRNA HTS libraries were highly expressed in flowers. Additionally, the expression of the precursors of the organ-specific miRNAs was analyzed. Degradome sequencing data-based approach was employed to identify the targets of the organ-specific miRNAs. The miRNA–target interactions were used for network construction. Subnetwork analysis unraveled some novel regulatory cascades, such as the feedback regulation mediated by miR161, the potential self-regulation of the genes miR172, miR396, miR398 and miR860, and the miR863-guided cleavage of the SERRATE transcript. Our bioinformatics survey expanded the organ-specific miRNA–target list in Arabidopsis, and could deepen the biological view of the miRNA expression and their regulatory roles.

Rosette iron deficiency transcript and microRNA profiling reveals links between copper and iron homeostasis in Arabidopsis thaliana

Iron (Fe) is an essential plant micronutrient, and its deficiency limits plant growth and development on alkaline soils. Under Fe deficiency, plant responses include up-regulation of genes involved in Fe uptake from the soil. However, little is known about shoot responses to Fe deficiency. Using microarrays to probe gene expression in Kas-1 and Tsu-1 ecotypes of Arabidopsis thaliana, and comparison with existing Col-0 data, revealed conserved rosette gene expression responses to Fe deficiency. Fe-regulated genes included known metal homeostasis-related genes, and a number of genes of unknown function. Several genes responded to Fe deficiency in both roots and rosettes. Fe deficiency led to up-regulation of Cu,Zn superoxide dismutase (SOD) genes CSD1 and CSD2, and down-regulation of FeSOD genes FSD1 and FSD2. Eight microRNAs were found to respond to Fe deficiency. Three of these (miR397a, miR398a, and miR398b/c) are known to regulate transcripts of Cu-containing proteins, and were down-regulated by Fe deficiency, suggesting that they could be involved in plant adaptation to Fe limitation. Indeed, Fe deficiency led to accumulation of Cu in rosettes, prior to any detectable decrease in Fe concentration. ccs1 mutants that lack functional Cu,ZnSOD proteins were prone to greater oxidative stress under Fe deficiency, indicating that increased Cu concentration under Fe limitation has an important role in oxidative stress prevention. The present results show that Cu accumulation, microRNA regulation, and associated differential expression of Fe and CuSOD genes are coordinated responses to Fe limitation.

High-throughput sequencing and characterization of the small RNA transcriptome reveal features of novel and conserved microRNAs in Panax ginseng

microRNAs (miRNAs) play vital regulatory roles in many organisms through direct cleavage of transcripts, translational repression, or chromatin modification. Identification of miRNAs has been carried out in various plant species. However, no information is available for miRNAs from Panax ginseng, an economically significant medicinal plant species. Using the next generation high-throughput sequencing technology, we obtained 13,326,328 small RNA reads from the roots, stems, leaves and flowers of P. ginseng. Analysis of these small RNAs revealed the existence of a large, diverse and highly complicated small RNA population in P. ginseng. We identified 73 conserved miRNAs, which could be grouped into 33 families, and 28 non-conserved ones belonging to 9 families. Characterization of P. ginseng miRNA precursors revealed many features, such as production of two miRNAs from distinct regions of a precursor, clusters of two precursors in a transcript, and generation of miRNAs from both sense and antisense transcripts. It suggests the complexity of miRNA production in P. gingseng. Using a computational approach, we predicted for the conserved and non-conserved miRNA families 99 and 31 target genes, respectively, of which eight were experimentally validated. Among all predicted targets, only about 20% are conserved among various plant species, whereas the others appear to be non-conserved, indicating the diversity of miRNA functions. Consistently, many miRNAs exhibited tissue-specific expression patterns. Moreover, we identified five dehydration- and ten heat-responsive miRNAs and found the existence of a crosstalk among some of the stress-responsive miRNAs. Our results provide the first clue to the elucidation of miRNA functions in P. ginseng.

Analyzing the microRNA Transcriptome in Plants Using Deep Sequencing Data

MicroRNAs (miRNAs) are 20- to 24-nucleotide endogenous small RNA molecules emerging as an important class of sequence-specific, trans-acting regulators for modulating gene expression at the post-transcription level. There has been a surge of interest in the past decade in identifying miRNAs and profiling their expression pattern using various experimental approaches. In particular, ultra-deep sampling of specifically prepared low-molecular-weight RNA libraries based on next-generation sequencing technologies has been used successfully in diverse species. The challenge now is to effectively deconvolute the complex sequencing data to provide comprehensive and reliable information on the miRNAs, miRNA precursors, and expression profile of miRNA genes. Here we review the recently developed computational tools and their applications in profiling the miRNA transcriptomes, with an emphasis on the model plant Arabidopsis thaliana. Highlighted is also progress and insight into miRNA biology derived from analyzing available deep sequencing data.