Antiquity of microRNAs and their targets in land plants

'Movers and shakers' in the regulation of fruit ripening: a cross-dissection of climacteric versus non-climacteric fruit

Fruit ripening is a complex and highly coordinated developmental process involving the expression of many ripening-related genes under the control of a network of signalling pathways. The hormonal control of climacteric fruit ripening, especially ethylene perception and signalling transduction in tomato has been well characterized. Additionally, great strides have been made in understanding some of the major regulatory switches (transcription factors such as RIPENING-INHIBITOR and other transcriptional regulators such as COLOURLESS NON-RIPENING, TOMATO AGAMOUS-LIKE1 and ETHYLENE RESPONSE FACTORs), that are involved in tomato fruit ripening. In contrast, the regulatory network related to non-climacteric fruit ripening remains poorly understood. However, some of the most recent breakthrough research data have provided several lines of evidences for abscisic acid- and sucrose-mediated ripening of strawberry, a non-climacteric fruit model. In this review, we discuss the most recent research findings concerning the hormonal regulation of fleshy fruit ripening and their cross-talk and the future challenges taking tomato as a climacteric fruit model and strawberry as a non-climacteric fruit model. We also highlight the possible contribution of epigenetic changes including the role of plant microRNAs, which is opening new avenues and great possibilities in the fields of fruit-ripening research and postharvest biology.

Genome-wide comparison of microRNAs and their targeted transcripts among leaf, flower and fruit of sweet orange

BACKGROUND

In plants, microRNAs (miRNAs) regulate gene expression mainly at the post-transcriptional level. Previous studies have demonstrated that miRNA-mediated gene silencing pathways play vital roles in plant development. Here, we used a high-throughput sequencing approach to characterize the miRNAs and their targeted transcripts in the leaf, flower and fruit of sweet orange.

RESULTS

A total of 183 known miRNAs and 38 novel miRNAs were identified. An in-house script was used to identify all potential secondary siRNAs derived from miRNA-targeted transcripts using sRNA and degradome sequencing data. Genome mapping revealed that these miRNAs were evenly distributed across the genome with several small clusters, and 69 pre-miRNAs were co-localized with simple sequence repeats (SSRs). Noticeably, the loop size of pre-miR396c was influenced by the repeat number of CUU unit. The expression pattern of miRNAs among different tissues and developmental stages were further investigated by both qRT-PCR and RNA gel blotting. Interestingly, Csi-miR164 was highly expressed in fruit ripening stage, and was validated to target a NAC transcription factor. This study depicts a global picture of miRNAs and their target genes in the genome of sweet orange, and focused on the comparison among leaf, flower and fruit tissues.

CONCLUSIONS

This study provides a global view of miRNAs and their target genes in different tissue of sweet orange, and focused on the identification of miRNA involved in the regulation of fruit ripening. The results of this study lay a foundation for unraveling key regulators of orange fruit development and ripening on post-transcriptional level.

Composition and expression of conserved microRNA genes in diploid cotton (Gossypium) species

MicroRNAs are ubiquitous in plant genomes but vary greatly in their abundance within and conservation among plant lineages. To gain insight into the evolutionary birth/death dynamics of microRNA families, we sequenced small RNA and 5′-end PARE libraries generated from two closely related species of Gossypium. Here, we demonstrate that 33 microRNA families, with similar copy numbers and average evolutionary rates, are conserved in the two congeneric cottons. Analysis of the presence/absence of these microRNA families in other land plants sheds light on their depth of phylogenetic origin and lineage-specific loss/gain. Conserved microRNA families in Gossypium exhibit a striking interspecific asymmetry in expression, potentially connected to relative proximity to neighboring transposable elements. A complex correlated expression pattern of microRNA target genes with their controlling microRNAs indicates that possible functional divergence of conserved microRNA families can also exist even within a single plant genus.

Recovery of dicer-like 1-late flowering phenotype by miR172 expressed by the noncanonical DCL4-dependent biogenesis pathway

MicroRNAs (miRNAs) act as down-regulators of gene expression, and play a dominant role in eukaryote development. In Arabidopsis thaliana, DICER-LIKE 1 (DCL1) is the main processor in miRNA biogenesis, and dcl1 mutants show various developmental defects at the early stage of embryogenesis or at gamete formation. However, miRNAs responsible for the respective developmental stages of the dcl1 defects have not been identified. Here, we developed a DCL1-independent miRNA expression system using the unique DCL4-dependent miRNA, miR839. By replacing the mature sequence in the miR839 precursor sequence with that of miR172, one of the most widely conserved miRNAs in angiosperms, we succeeded in expressing miR172 from a chimeric miR839 precursor in dcl1-7 plants and observed the repression of miR172 target gene expression. In parallel, the DCL4-dependent miR172 expression rescued the late flowering phenotype of dcl1-7 by acceleration of flowering. We established the DCL1-independent miRNA expression system, and revealed that the reduction of miR172 expression is responsible for the dcl1-7 late flowering phenotype.

Expression patterns of conserved microRNAs in the male gametophyte of loblolly pine (Pinus taeda)

MicroRNAs (miRNAs) are small RNAs that regulate genes involved in various aspects of plant development, but their presence and expression patterns in the male gametophytes of gymnosperms have not yet been established. Therefore, this study identified and compared the expression patterns of conserved miRNAs from two stages of the male gametophyte of loblolly pine (Pinus taeda), which are the mature (ungerminated) and germinated pollen. Microarray was used to identify conserved miRNAs that varied in expression between these two stages of the loblolly pine male gametophyte. Forty-seven conserved miRNAs showed significantly different expression levels between mature and germinated loblolly pine pollen. In particular, miRNAs representing 14 and 8 families were up- and down-regulated in germinated loblolly pine pollen, respectively. qRT-PCR was used to validate their expression patterns using representative miRNAs. Target genes and proteins were identified using psRNATarget program. Predicted targets of the 22 miRNA families belong mostly to classes of genes involved in defense/stress response, metabolism, regulation, and signaling. qRT-PCR was also used to validate the expression patterns of representative target genes. This study shows that conserved miRNAs are expressed in mature and germinated loblolly pine pollen. Many of these miRNAs are differentially expressed, which indicates that the two stages of the male gametophyte examined are regulated at the miRNA level. This study also expands our knowledge of the male gametophytes of seed plants by providing insights on some similarities and differences in the types and expression patterns of conserved miRNAs between loblolly pine with those of rice and Arabidopsis.

Down-regulation of AUXIN RESPONSE FACTORS 6 and 8 by microRNA 167 leads to floral development defects and female sterility in tomato

Auxin regulates the expression of diverse genes that affect plant growth and development. This regulation requires AUXIN RESPONSE FACTORS (ARFs) that bind to the promoter regions of these genes. ARF6 and ARF8 in Arabidopsis thaliana are required to promote inflorescence stem elongation and late stages of petal, stamen, and gynoecium development. All seed plants studied thus far have ARF6 and ARF8 orthologues as well as the microRNA miR167, which targets ARF6 and ARF8. Whether these genes have broadly conserved roles in flower development is not known. To address this question, the effects of down-regulation of ARF6 and ARF8 were investigated through transgenic expression of Arabidopsis MIR167a in tomato, which diverged from Arabidopsis before the radiation of dicotyledonous plants approximately 90-112 million years ago. The transgenic tomato plants overexpressing MIR167a exhibited reductions in leaf size and internode length as well as shortened petals, stamens, and styles. More significantly, the transgenic plants were female-sterile as a result of failure of wild-type pollen to germinate on the stigma surface and/or to grow through the style. RNA-Seq analysis identified many genes with significantly altered expression patterns, including those encoding products with functions in 'transcription regulation', 'cell wall' and 'lipid metabolism' categories. Putative orthologues of a subset of these genes were also differentially expressed in Arabidopsis arf6 arf8 mutant flowers. These results thus suggest that ARF6 and ARF8 have conserved roles in controlling growth and development of vegetative and flower organs in dicots.

Expression of microRNAs and their targets regulates floral development in tobacco (Nicotiana tabacum)

MicroRNAs (miRNAs) are an extensive class of endogenous posttranscriptional gene regulators that function to mediate gene expression by cleaving target mRNAs or by preventing protein translation. Because of their importance in mediating gene regulation, identifying and elucidating the function of miRNAs have been the primary focus of many researchers. Now that many miRNAs have been identified and assessed for their functionality, the next step is to create expression profiles for miRNAs, so that gene expression studies can be further enhanced with knowledge of the basal expression levels of miRNAs and their targets. In a previous study, 259 putative miRNAs were identified in tobacco, in which 11 of them were confirmed. The primary goal of this study was to further expand on that study and create an expression profile for nine miRNAs and their targets in a tissue-specific manner in tobacco. We chose to study miRNAs that largely play a role in floral development and nutrient stress response. The results of our study show that all tested miRNAs and their targets were expressed in a differential manner. The results of our study also show that out of the tested miRNAs and their targets, miR159, miR157, miR167, miR172, and superoxide dismutase were expressed the highest, suggesting that these genes may play a vital role in the growth and development of tobacco. Corrected expression of miRNAs and their targets regulates floral development.

Identification and functional analysis of novel and conserved microRNAs in tomato

MicroRNAs are ~22 nt non-coding endogenous RNAs which play important regulation roles in various species. By using homology-based computational research, 14 novel and conserved tomato miRNAs belonging to ten families were identified from EST, GSS, and nucleotide sequences. Real-time PCR analysis of these miRNAs demonstrated their expression in tomato afterwards. Meanwhile, a total of 36 potential targets were predicted for the ten miRNAs using psRNATarget. The target genes were mainly involved in metabolism, transmembrane transport, stress response, and transcription regulation. According to our experiment, miR398 was down-regulated on different levels under biotic and abiotic stresses, suggesting that miR398 might be involved in tomato stress regulatory network. Our results supplement the findings of tomato miRNAs and also suggest crucial regulatory functions of miRNAs in stress responses.

Genome-wide identification of microRNAs in Medicago truncatula by high-throughput sequencing

MicroRNAs (miRNAs) are small, endogenous RNAs that play important regulatory roles in development and stress response in plants by negatively regulating gene expression post-transcriptionally. Medicago truncatula has been used as a model plant to study functional genomics of legume plants. It has also been widely used to functionally study miRNAs. Identification of miRNAs at the whole-genome level is essential for functional characterization of miRNAs in plants. High-throughput sequencing is a powerful technology to identify miRNAs. In this chapter, the methods used for construction of a small RNA library and high-throughput sequencing involving total RNA isolation, small RNA purification, adapter ligation, reverse transcription, PCR amplification, and Solexa sequencing are described. Bioinformatics and analysis of differential expression of miRNAs including primary disposal, miRNA identification, target prediction, and expression analysis are also discussed. These methodologies associated with identification and functional characterization of miRNAs may provide useful tools for readers to study miRNAs in plants in general and Medicago truncatula in particular.

miRNAs trigger widespread epigenetically activated siRNAs from transposons in Arabidopsis

In plants, post-transcriptional gene silencing (PTGS) is mediated by DICER-LIKE1 (DCL1)-dependent miRNAs, that also trigger 21-nt secondary siRNA via RNA DEPENDENT RNA POLYMERASE6 (RDR6), DCL4, and ARGONAUTE1 (AGO1)^1–3, while transcriptional gene silencing (TGS) of transposons is mediated by 24-nt heterochromatic (het)siRNA RDR2, DCL3 and AGO4⁴. Transposons can also give rise to abundant 21-nt “epigenetically activated” small interfering RNAs (easiRNAs) in DECREASE IN DNA METHYLATION1 (ddm1) and DNA METHYLTRANSFERASE1 (met1) mutants, as well as in the vegetative nucleus of pollen grains⁵, and in dedifferentiated plant cell cultures⁶. Here we show that easiRNAs resemble secondary siRNAs, in that thousands of transposon transcripts are specifically targeted by more than fifty miRNAs for cleavage and processing by RDR6. Loss of RDR6, DCL4 or DCL1 in a ddm1 background results in loss of 21-nt easiRNA, and severe infertility, but 24-nt hetsiRNA are partially restored, supporting an antagonistic relationship between PTGS and TGS. Thus miRNA-directed easiRNA biogenesis is a latent mechanism that specifically targets transposon transcripts, but only when they are epigenetically reactivated during reprogramming of the germline. This ancient recognition mechanism may have been retained both by transposons to evade long-term heterochromatic silencing, and by their hosts for genome defence.

Identification and characterization of the microRNA transcriptome of a moth orchid Phalaenopsis aphrodite

Orchids display unique phenotypes, functional characteristics and ecological adaptations that are not found in model plants. In this study, we aimed to characterize the microRNA (miRNA) transcriptome and identify species- and tissue-specific miRNAs in Phalaenopsis aphrodite. After data filtering and cleanup, a total of 59,387,374 reads, representing 1,649,996 unique reads, were obtained from four P. aphrodite small RNA libraries. A systematic bioinformatics analysis pipeline was developed that can be used for miRNA and precursor mining, and target gene prediction in non-model plants. A total of 3,251 unique reads for 181 known plant miRNAs (belonging to 88 miRNA families), 23 new miRNAs and 91 precursors were identified. All the miRNA star sequences (miRNA*), the complementary strands of miRNA that from miRNA/miRNA* duplexes, of the predicted new miRNAs were detected in our small RNA libraries, providing additional evidence for their existence as new miRNAs in P. aphrodite. Furthermore, 240 potential miRNA-targets that appear to be involved in many different biological activities and molecular functions, especially transcription factors, were identified, suggesting that miRNAs can impact multiple processes in P. aphrodite. We also verified the cleavage sites for six targets using RNA ligase-mediated rapid amplification of 5' ends assay. The results provide valuable information about the composition, expression and function of miRNA in P. aphrodite, and will aid functional genomics studies of orchids.

Evolutionary history of plant microRNAs

microRNAs (miRNAs) are short noncoding regulatory genes that perform important roles in plant development and physiology. With the increasing power of next generation sequencing technologies and the development of bioinformatic tools, there has been a dramatic increase in the number of studies surveying the miRNAomes of plant species, which has led to an explosion in the number of described miRNAs. Unfortunately, very many of these new discoveries have been incompletely annotated and thus fail to discriminate genuine miRNAs from small interfering RNAs (siRNAs), fragments of longer RNAs, and random sequence. We review the published repertoire of plant miRNAs, discriminating those that have been correctly annotated. We use these data to explore prevailing hypotheses on the tempo and mode of miRNA evolution within the plant kingdom.

The involvement of InMIR167 in the regulation of expression of its target gene InARF8, and their participation in the vegetative and generative development of Ipomoea nil plants

The plant hormone auxin plays a critical role in regulating plant growth and development. Recent advances have been made that having improved our understanding of auxin response pathways, primarily by characterizing the genes encoding auxin response factors (ARFs) in Arabidopsis. In addition, the expression of some ARFs is regulated by microRNAs (miRNAs). In Arabidopsis thaliana, ARF6 and ARF8 are targeted by miR167, whereas ARF10, ARF16 and ARF17 are targeted by miR160. Nevertheless, little is known about any possible interactions between miRNAs and the auxin signaling pathway during plant development. In this study, we isolated the miR167 target gene InARF8 cDNA from the cotyledons of the short day plant (SDP) Ipomoea nil (named also Pharbitis nil). Additionally, the In-miR167 precursor was identified from the I. nil EST database and analyses of InARF8 mRNA, In-pre-miR167 and mature miR167 accumulation in the plant's vegetative and generative organs were performed. The identified cDNA of InARF8 contains a miR167 complementary sequence and shows significant similarity to ARF8 cDNAs of other plant species. The predicted amino acid sequence of InARF8 includes all of the characteristic domains for ARF family transcription factors (B3 DNA-binding domain, AUX/IAA-CTD and a glutamine-rich region). Quantitative RT-PCR reactions and in situ hybridization indicated that InARF8 was expressed primarily in the shoot apices, leaf primordia and hypocotyls of I. nil seedlings, as well as in flower pistils and petals. The InARF8 transcript level increased consistently during the entire period of pistil development, whereas in the stamens, the greatest transcriptional activity occurred only during the intensive elongation phase. Additionally, an expression analysis of both the precursor In-pre-miR167 molecules identified and mature miRNA was performed. We observed that, in most of the organs examined, the InARF8 expression pattern was opposite to that of MIR167, indicating that the gene's activity was regulated by mRNA cleavage. Our findings suggested that InARF8 and InMIR167 participated in the development of young tissues, especially the shoot apices and flower elements. The main function of MIR167 appears to be to regulate InARF8 organ localization.

Characterization of regulatory mechanism of Poncirus trifoliata microRNAs on their target genes with an integrated strategy of newly developed PPM-RACE and RLM-RACE

MicroRNAs (miRNAs) play an important role in post-transcriptional gene regulation that involved various biological and metabolic processes. Many extensive studies have been done in model plant species, to discover miRNAs' regulating expression of their target genes and analyze their functions. But, the function of Poncirus trifoliata miRNAs has not been properly investigated. In this study, we employed the RNA ligase-mediated 5' rapid amplification of cDNA ends (RLM-RACE) and the newly developed method called poly (A) polymerase-mediated 3' rapid amplification of cDNA ends (PPM-RACE), which mapped the cleavage site of target mRNAs and detected expression patterns of cleaved fragments that could in turn indicate the regulatory functions of the miRNAs on their target genes. Furthermore, the spatiotemporal expression levels of target genes were analyzed by qRT-PCR, with exhibiting different expression trends from their corresponding miRNAs, thus indicating the cleavage mode of miRNAs on their target genes. The expression patterns of miRNAs, their target mRNAs and cleaved target mRNAs in different organs of juvenile and adult trifoliate orange were studied. The results showed that the expression of miRNAs and their target mRNAs was in a trade-off trend. When the miRNA expression was high, its corresponding target mRNA expression was low, while the cleaved target mRNA expression was high; when the miRNA expression was low, its target mRNA expression was high, while the expression of cleaved target mRNAs follows that of the miRNA. The validation of the cleavage site of target mRNAs and the detection of expression patterns of cleaved fragments can further broaden the knowledge of small RNA-mediated regulation in P. trifoliate.

Differential expression of miRNAs in Brassica napus root following infection with Plasmodiophora brassicae

Canola (oilseed rape, Brassica napus L.) is susceptible to infection by the biotrophic protist Plasmodiophora brassicae, the causal agent of clubroot. To understand the roles of microRNAs (miRNAs) during the post-transcriptional regulation of disease initiation and progression, we have characterized the changes in miRNA expression profiles in canola roots during clubroot disease development and have compared these to uninfected roots. Two different stages of clubroot development were targeted in this miRNA profiling study: an early time of 10-dpi for disease initiation and a later 20-dpi, by which time the pathogen had colonized the roots (as evident by visible gall formation and histological observations). P. brassicae responsive miRNAs were identified and validated by qRT-PCR of miRNAs and the subsequent validation of the target mRNAs through starBase degradome analysis, and through 5' RLM-RACE. This study identifies putative miRNA-regulated genes with roles during clubroot disease initiation and development. Putative target genes identified in this study included: transcription factors (TFs), hormone-related genes, as well as genes associated with plant stress response regulation such as cytokinin, auxin/ethylene response elements. The results of our study may assist in elucidating the role of miRNAs in post-transcriptional regulation of target genes during disease development and may contribute to the development of strategies to engineer durable resistance to this important phytopathogen.

Analysis of biochemical variations and microRNA expression in wild ( Ipomoea campanulata ) and cultivated ( Jacquemontia pentantha ) species exposed to in vivo water stress

The current study analyses few important biochemical parameters and microRNA expression in two closely related species (wild but tolerant Ipomoea campanulata L. and cultivated but sensitive Jacquemontia pentantha Jacq.G.Don) exposed to water deficit conditions naturally occurring in the field. Under soil water deficit, both the species showed reduction in their leaf area and SLA as compared to well-watered condition. A greater decrease in chlorophyll was noticed in J. pentantha (~50 %) as compared to I. campanulata (20 %) under stress. By contrast, anthocyanin and MDA accumulation was greater in J. pentantha as compared to I. campanulata. Multiple isoforms of superoxide dismutases (SODs) with differing activities were observed under stress in these two plant species. CuZnSOD isoforms showed comparatively higher induction (~10-40 %) in I. campanulata than J. pentantha. MicroRNAs, miR398, miR319, miR395 miR172, and miR408 showed opposing expression under water deficit in these two plant species. Expression of miR156, miR168, miR171, miR172, miR393, miR319, miR396, miR397 and miR408 from either I. campanulata or J. pentantha or both demonstrated opposite pattern of expression to that of drought stressed Arabidopsis. The better tolerance of the wild species (I. campanulata) to water deficit could be attributed to lesser variations in chlorophyll and anthocyanin levels; and relatively higher levels of SODs than J. pentantha. miRNA expression was different in I. campanulata than J. pentantha.

Computational identification of microRNAs in the strawberry (Fragaria x ananassa) genome sequence and validation of their precise sequences by miR-RACE

In plants, microRNAs (miRNAs) play significant roles in post-transcriptional gene regulation and have been found to control many genes involved in different biological and metabolic processes. Extensive studies were carried out to discover miRNAs and analyze their functions in model plant species, such as in Arabidopsis and rice that have been reported. In this research, we used bioinformatics to predict microRNAs in an important strawberry rootstock cultivar to discover and validate precise sequences of microRNAs in strawberry. By adopting a range of filtering criteria, we obtained 59 potential miRNAs belonging to 40 miRNA families from the Fragaria vesca genome. Using two specific 5' and 3' miRNA RACE PCR reactions and a sequence-directed cloning method, we accurately determined 34 precise sequences of candidate miRNAs, while six other sequences exhibited some minor divergence in their termini nucleotides, and 19 miRNAs that could not be cloned owing to expression abundance may be too low or these mirRNAs predicted could not be existing in strawberry. Potential target genes were further predicted for the miRNAs above. The expression of the 16 miRNAs unreported and having exact sequences and their targets by experiment could be detected in different tissues of strawberry ranging from roots, stems, leaves, flowers and fruits by qRT-PCR and some of them showed differential expression in various tissues. The functional analysis of 16 miRNAs and their targets was carried out. Finally, we conclude that there are 34 mirRNAs in strawberry and their targets play vital roles not only in growth and development, but also in diverse physiological processes. These results show that regulatory miRNAs exist in agronomically important strawberry and might have an important function in strawberry growth and development.

Molecular identification and analysis of Cd-responsive microRNAs in rice

Cadmium (Cd) is a non-essential heavy metal with high toxicity to plants. MicroRNAs (miRNAs) are a class of small non-coding RNAs that play important roles in plant abiotic stress responses. To investigate whether miRNAs function in Cd stress response, miRNA expression profiles in rice ( Oryza sativa ) under Cd stress were monitored using microarray assays. A total of 12 Cd-responsive novel miRNAs predicted previously were identified, of which 4 were further validated experimentally. A total of 44 target genes were predicted for the Cd-responsive miRNAs, many of which appeared to regulate gene networks mediating environmental stresses. Several target genes were validated to show a reciprocal regulation by miRNAs. A transgenic approach was also used to determine the role of miRNAs in rice response to Cd stress. Overexpression of miR192 retarded seed germination and seedling growth under Cd stress. These results implied the role of novel miRNAs in the involvement of Cd tolerance of rice.

Phylogenetic analysis reveals conservation and diversification of micro RNA166 genes among diverse plant species

Similar to the majority of the microRNAs, mature miR166s are derived from multiple members of MIR166 genes (precursors) and regulate various aspects of plant development by negatively regulating their target genes (Class III HD-ZIP). The evolutionary conservation or functional diversification of miRNA166 family members remains elusive. Here, we show the phylogenetic relationships among MIR166 precursor and mature sequences from three diverse model plant species. Despite strong conservation, some mature miR166 sequences, such as ppt-miR166m, have undergone sequence variation. Critical sequence variation in ppt-miR166m has led to functional diversification, as it targets non-HD-ZIPIII gene transcript (s). MIR166 precursor sequences have diverged in a lineage specific manner, and both precursors and mature osa-miR166i/j are highly conserved. Interestingly, polycistronic MIR166s were present in Physcomitrella and Oryza but not in Arabidopsis. The nature of cis-regulatory motifs on the upstream promoter sequences of MIR166 genes indicates their possible contribution to the functional variation observed among miR166 species.

Peculiar evolutionary history of miR390-guided TAS3-like genes in land plants

PCR-based approach was used as a phylogenetic profiling tool to probe genomic DNA samples from representatives of evolutionary distant moss taxa, namely, classes Bryopsida, Tetraphidopsida, Polytrichopsida, Andreaeopsida, and Sphagnopsida. We found relatives of all Physcomitrella patens miR390 and TAS3-like loci in these plant taxa excluding Sphagnopsida. Importantly, cloning and sequencing of Marchantia polymorpha genomic DNA showed miR390 and TAS3-like sequences which were also found among genomic reads of M. polymorpha at NCBI database. Our data suggest that the ancient plant miR390-dependent TAS molecular machinery firstly evolved to target AP2-like mRNAs in Marchantiophyta and only then both ARF- and AP2-specific mRNAs in mosses. The presented analysis shows that moss TAS3 families may undergone losses of tasiAP2 sites during evolution toward ferns and seed plants. These data confirm that miR390-guided genes coding for ARF- and AP2-specific ta-siRNAs have been gradually changed during land plant evolution.

Construction of microRNA functional families by a mixture model of position weight matrices

MicroRNAs (miRNAs) are small regulatory molecules that repress the translational processes of their target genes by binding to their 3′ untranslated regions (3′ UTRs). Because the target genes are predominantly determined by their sequence complementarity to the miRNA seed regions (nucleotides 2–7) which are evolutionarily conserved, it is inferred that the target relationships and functions of the miRNA family members are conserved across many species. Therefore, detecting the relevant miRNA families with confidence would help to clarify the conserved miRNA functions, and elucidate miRNA-mediated biological processes. We present a mixture model of position weight matrices for constructing miRNA functional families. This model systematically finds not only evolutionarily conserved miRNA family members but also functionally related miRNAs, as it simultaneously generates position weight matrices representing the conserved sequences. Using mammalian miRNA sequences, in our experiments, we identified potential miRNA groups characterized by similar sequence patterns that have common functions. We validated our results using score measures and by the analysis of the conserved targets. Our method would provide a way to comprehensively identify conserved miRNA functions.

Deep sequencing discovery of novel and conserved microRNAs in wild type and a white-flesh mutant strawberry

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression by base pairing to mRNA target sequences, and play crucial roles in plant development and stress responses. The knowledge on post-transcriptional regulation by miRNAs in strawberry is rather limited so far. In order to understand the role of miRNA in the molecular control during strawberry fruit development, small RNA libraries were constructed from fruits at the turning stage of strawberry cultivar 'Sachinoka' and its white-flesh mutant by using the Solexa platform. One hundred and twenty conserved miRNAs belonging to 27 miRNA families and 33 putative novel strawberry miRNAs were identified in both libraries. Their target genes were predicted using the Fragaria vesca genome. Nine of all miRNAs showed significant expression differences between two types of samples. Four miRNAs were up-regulated and five were down-regulated in white-flesh mutant. The sequencing results were partially validated by quantitative RT-PCR. Among them, the expression of miR399a shows the biggest change between the two samples. The prediction of its target gene showed that miR399 may play an important role in phosphate homeostasis of strawberry fruits. Furthermore, we deduce that the expression of miR399 has negative correlation with the content of sugars.

Comprehensive analyses of microRNA gene evolution in paleopolyploid soybean genome

miRNA genes are thought to undergo quick birth and death processes in genomes and the emergence of a MIRNA-like hairpin provides the base for functional miRNA gene formation. However, the factors affecting the formation of an active miRNA gene from a MIRNA-like hairpin within a genome remain unclear. We performed a genome-wide investigation of MIRNA-like hairpin accumulation, expression, structural changes and relationships with annotated genomic features in the paleopolyploid soybean genome. Our results showed that adjacent gene and transposable element content, rates of genetic recombination at location of emergence, along with its own gene structure divergence greatly affected miRNA gene evolution. Further investigation suggested that miRNA genes from different duplication sources followed distinct evolutionary trajectories and that the accumulation of MIRNA-like hairpins might be an important factor in causing long terminal repeat retrotransposons to lose activity during genome evolution.

Computational identification of conserved microRNAs and their putative targets in the Hypericum perforatum L. flower transcriptome

MicroRNAs (miRNAs) have recently emerged as important regulators of gene expression in plants. Many miRNA families and their targets have been extensively studied in model species and major crops. We have characterized mature miRNAs along with their precursors and potential targets in Hypericum to generate a comprehensive list of conserved miRNA families and to investigate the regulatory role of selected miRNAs in biological processes that occur in the flower. St. John's wort (Hypericum perforatum L., 2n = 4x = 32), a medicinal plant that produces pharmaceutically important metabolites with therapeutic activities, was chosen because it is regarded as an attractive model system for the study of apomixis. A computational in silico prediction of structure, in combination with an in vitro validation, allowed us to identify 7 pre-miRNAs, including miR156, miR166, miR390, miR394, miR396, and miR414. We demonstrated that H. perforatum flowers share highly conserved miRNAs and that these miRNAs potentially target dozens of genes with a wide range of molecular functions, including metabolism, response to stress, flower development, and plant reproduction. Our analysis paves the way toward identifying flower-specific miRNAs that may differentiate the sexual and apomictic reproductive pathways.

Genome-wide identification, molecular cloning, expression profiling and posttranscriptional regulation analysis of the Argonaute gene family in Salvia miltiorrhiza, an emerging model medicinal plant

Background

Argonaute (AGO) is the core component of RNA-induced silencing complex. The AGO gene family has been analyzed in various plant species; however, there is no report about AGOs in the well-known Traditional Chinese Medicine (TCM) plant, Salvia miltiorrhiza.

Results

Through a genome-wide analysis, we identified ten SmAGO genes in S. miltiorrhiza. Full-length cDNAs of all SmAGOs were subsequently cloned and sequenced. These SmAGOs were characterized using a comprehensive approach. Sequence features, gene structures and conserved domains were analyzed by the comparison of SmAGOs and AtAGOs. Phylogenetic relationships among AGO proteins from S. miltiorrhiza, Arabidopsis and rice were revealed. The expression levels of SmAGO genes in various tissues of S. miltiorrhiza were investigated. The results implied that some SmAGOs, such as SmAGO1, SmAGO2, SmAGO3, SmAGO7 and SmAGO10, probably played similar roles as their counterparts in Arabidopsis; whereas the others could be more species-specialized. It suggests the conservation and diversity of AGOs in plants. Additionally, we identified a total of 24 hairpin structures, representing six miRNA gene families, to be miRNA precursors. Using the modified 5′-RACE method, we confirmed that SmAGO1 and SmAGO2 were targeted by S. miltiorrhiza miR168a/b and miR403, respectively. It suggests the conservation of AGO1-miR168 and AGO2-miR403 regulatory modules in S. miltiorrhiza and Arabidopsis.

Conclusions

This is the first attempt to explore SmAGOs and miRNAs in S. miltiorrhiza. The results provide useful information for further elucidation of gene silencing pathways in S. miltiorrhiza.

Abundances of microRNAs in human cells can be estimated as a function of the abundances of YRHB and RHHK tetranucleotides in these microRNAs as an ill-posed inverse problem solution

Mature microRNAs (miRNAs) are small endogenous non-coding RNAs 18–25 nt in length. They program the RNA Induced Silencing Complex (RISC) to make it inhibit either messenger RNAs or promoter DNAs. We have found that the mean abundance of miRNAs in Arabidopsis is correlated with the abundance of DRYD tetranucleotides near the 3′-end and the abundance of WRHB tetranucleotides in the center of the miRNA sequence. Based on this correlation, we have estimated miRNA abundances in seven organs of this plant, namely: inflorescences, stems, siliques, seedlings, roots, cauline, and rosette leaves. We have also found that the mean affinity of miRNAs for two proteins in the Argonaute family (Ago2 and Ago3) in man is correlated with the abundance of YRHB tetranucleotides near the 3′-end and that the preference of miRNAs for Ago2 is correlated with the abundance of RHHK tetranucleotides in the center of the miRNA sequence. This allowed us to obtain statistically significant estimates of miRNA abundances in human embryonic kidney cells, HEK293T. These findings in relation to two taxonomically distant entities (man and Arabidopsis) fit one another like pieces of a jigsaw puzzle, which allowed us to heuristically generalize them and state that the miRNA abundance in the human brain may be determined by the abundance of YRHB and RHHK tetranucleotides in these miRNAs.

Deep sequencing discovery of novel and conserved microRNAs in strawberry (Fragaria×ananassa)

In plants, microRNAs (miRNAs) play a critical role in post-transcriptional gene regulation and have been shown to control many genes involved in various biological and metabolic processes. There have been extensive studies to discover miRNAs and analyze their functions in model plant species, such as Arabidopsis and rice. Deep sequencing technologies have facilitated identification of species-specific or lowly expressed as well as conserved or highly expressed miRNAs in plants. In this research, we used Solexa sequencing to discover new miRNAs in cultivated strawberry (Fragaria×ananassa). A total of 23,282 ,309 reads representing 22,500 ,402 distinct sequences were obtained from a short RNA library generated from small RNAs extracted from strawberry fruit tissues. On the basis of sequence similarity and hairpin structure prediction, we found that 156,639 reads representing 153 sequences have good matches to known miRNAs. We also identified 37 novel miRNA candidates. These sequences had not been previously described in other plant species. Potential target genes were predicted for the majority of and novel miRNAs. These results show that regulatory miRNAs exist in the agriculturally important cultivated strawberry and may play an important role in its growth, development and response to disease.

The hot pepper (Capsicum annuum) microRNA transcriptome reveals novel and conserved targets: a foundation for understanding MicroRNA functional roles in hot pepper

MicroRNAs (miRNAs) are a class of non-coding RNAs approximately 21 nt in length which play important roles in regulating gene expression in plants. Although many miRNA studies have focused on a few model plants, miRNAs and their target genes remain largely unknown in hot pepper (Capsicum annuum), one of the most important crops cultivated worldwide. Here, we employed high-throughput sequencing technology to identify miRNAs in pepper extensively from 10 different libraries, including leaf, stem, root, flower, and six developmental stage fruits. Based on a bioinformatics pipeline, we successfully identified 29 and 35 families of conserved and novel miRNAs, respectively. Northern blot analysis was used to validate further the expression of representative miRNAs and to analyze their tissue-specific or developmental stage-specific expression patterns. Moreover, we computationally predicted miRNA targets, many of which were experimentally confirmed using 5′ rapid amplification of cDNA ends analysis. One of the validated novel targets of miR-396 was a domain rearranged methyltransferase, the major de novo methylation enzyme, involved in RNA-directed DNA methylation in plants. This work provides the first reliable draft of the pepper miRNA transcriptome. It offers an expanded picture of pepper miRNAs in relation to other plants, providing a basis for understanding the functional roles of miRNAs in pepper.

Evolutionary conservation of microRNA regulatory programs in plant flower development

MicroRNAs (miRNAs) are post-transcriptional regulators of growth and development in both plants and animals. Flowering is critical for the reproduction of angiosperms. Flower development entails the transition from vegetative growth to reproductive growth, floral organ initiation, and the development of floral organs. These developmental processes are genetically regulated by miRNAs, which participate in complex genetic networks of flower development. A survey of the literature shows that miRNAs, their specific targets, and the regulatory programs in which they participate are conserved throughout the plant kingdom. This review summarizes the role of miRNAs and their targets in the regulation of gene expression during the floral developmental phase, which includes the floral transition stage, followed by floral patterning, and then the development of floral organs. The conservation patterns observed in each component of the miRNA regulatory system suggest that these miRNAs play important roles in the evolution of flower development.

Detection and evolutionary analysis of soybean miRNAs responsive to soybean mosaic virus

MicroRNAs (miRNA) are a class of non-coding RNAs that have important gene regulatory roles in various organisms. However, the miRNAs involved in soybean's response to soybean mosaic virus (SMV) are unknown. To identify novel miRNAs and biotic-stress regulated small RNAs that are involved in soybean's response to SMV, two small RNA libraries were constructed from mock-inoculated and SMV-infected soybean leaves and sequenced. This led to the discovery of 179 miRNAs, representing 52 families, among which five miRNAs belonging to three families were novel miRNAs in soybean. A large proportion (71.5 %) of miRNAs arose from segmental duplication, similar to the process that drives the evolution of protein-coding genes. In addition, we predicted 346 potential targets of these identified miRNAs, and verified 12 targets by modified 5'-RACE analysis. Finally, three miRNAs (miR160, miR393 and miR1510) that are involved in plant resistance were observed to respond to SMV infection. The interaction between miRNAs and resistance-related genes provides a novel mechanism for pathogens to evade host recognition.

Morphogenesis of simple leaves: regulation of leaf size and shape

Plants produce new organs throughout their life span. Leaves first initiate as rod-like structures protruding from the shoot apical meristem, while they need to pass through different developmental stages to become the flat organ specialized in photosynthesis. Leaf morphogenesis is an active process regulated by many genes and pathways that can generate organs with a wide variety of sizes and shapes. Important differences in leaf architecture can be seen among different species, but also in single individuals. A key aspect of leaf morphogenesis is the precise control of cell proliferation. Modification or manipulation of this process may lead to leaves with different sizes and shapes, and changes in the organ margins and curvature. Many genes required for leaf development have been identified in Arabidopsis thaliana, and the mechanisms underlying leaf morphogenesis are starting to be unraveled at the molecular level.

Genome-wide identification and characterization of microRNA genes and their targets in flax (Linum usitatissimum): Characterization of flax miRNA genes

MicroRNAs (miRNAs) are small (20-24 nucleotide long) endogenous regulatory RNAs that play important roles in plant growth and development. They regulate gene expression at the post-transcriptional level by translational repression or target degradation and gene silencing. In this study, we identified 116 conserved miRNAs belonging to 23 families from the flax (Linum usitatissimum L.) genome using a computational approach. The precursor miRNAs varied in length; while most of the mature miRNAs were 21 nucleotide long, intergenic and showed conserved signatures of RNA polymerase II transcripts in their upstream regions. Promoter region analysis of the flax miRNA genes indicated prevalence of MYB transcription factor binding sites. Four miRNA gene clusters containing members of three phylogenetic groups were identified. Further, 142 target genes were predicted for these miRNAs and most of these represent transcriptional regulators. The miRNA encoding genes were expressed in diverse tissues as determined by digital expression analysis as well as real-time PCR. The expression of fourteen miRNAs and nine target genes was independently validated using the quantitative reverse transcription PCR (qRT-PCR). This study suggests that a large number of conserved plant miRNAs are also found in flax and these may play important roles in growth and development of flax.

Genomic organization, phylogenetic comparison and differential expression of the SBP-box family genes in grape

BACKGROUND

The SBP-box gene family is specific to plants and encodes a class of zinc finger-containing transcription factors with a broad range of functions. Although SBP-box genes have been identified in numerous plants including green algae, moss, silver birch, snapdragon, Arabidopsis, rice and maize, there is little information concerning SBP-box genes, or the corresponding miR156/157, function in grapevine.

METHODOLOGY/PRINCIPAL FINDINGS

Eighteen SBP-box gene family members were identified in Vitis vinifera, twelve of which bore sequences that were complementary to miRNA156/157. Phylogenetic reconstruction demonstrated that plant SBP-domain proteins could be classified into seven subgroups, with the V. vinifera SBP-domain proteins being more closely related to SBP-domain proteins from dicotyledonous angiosperms than those from monocotyledonous angiosperms. In addition, synteny analysis between grape and Arabidopsis demonstrated that homologs of several grape SBP genes were found in corresponding syntenic blocks of Arabidopsis. Expression analysis of the grape SBP-box genes in various organs and at different stages of fruit development in V. quinquangularis 'Shang-24' revealed distinct spatiotemporal patterns. While the majority of the grape SBP-box genes lacking a miR156/157 target site were expressed ubiquitously and constitutively, most genes bearing a miR156/157 target site exhibited distinct expression patterns, possibly due to the inhibitory role of the microRNA. Furthermore, microarray data mining and quantitative real-time RT-PCR analysis identified several grape SBP-box genes that are potentially involved in the defense against biotic and abiotic stresses.

CONCLUSION

The results presented here provide a further understanding of SBP-box gene function in plants, and yields additional insights into the mechanism of stress management in grape, which may have important implications for the future success of this crop.

Effects of drought on the microtranscriptome of field-grown sugarcane plants

Sugarcane (Saccharum spp.) is the most promising crop for renewable energy. Among the diverse stresses that affect plant productivity, drought stress frequently causes losses in sugarcane fields. Although several studies have addressed plant responses to drought using controlled environments, plant responses under field conditions are largely unknown. Recently, microRNA (miRNA)-mediated post-transcriptional regulation has been described as an important and decisive component in vegetal development and stress resistance modulation. The role of miRNAs in sugarcane responses to drought under field conditions is currently not known. Two sugarcane cultivars differing in drought tolerance were grown in the field with and without irrigation (rainfed) for 7 months. By using small RNA deep sequencing, we were able to identify 18 miRNA families comprising 30 mature miRNA sequences. Among these families, we found 13 mature miRNAs that were differentially expressed in drought-stressed plants. Seven miRNAs were differentially expressed in both cultivars. The target genes for many of the differentially expressed mature miRNAs were predicted, and some of them were validated by quantitative reverse transcription PCR. Among the targets, we found transcription factors, transporters, proteins associated with senescence, and proteins involved with flower development. All of these data increase our understanding of the role of miRNAs in the complex regulation of drought stress in field-grown sugarcane, providing valuable tools to develop new sugarcane cultivars tolerant to drought stress.

MicroRNAs and their putative targets in Brassica napus seed maturation

Background

MicroRNAs (miRNAs) are 20–21 nucleotide RNA molecules that suppress the transcription of target genes and may also inhibit translation. Despite the thousands of miRNAs identified and validated in numerous plant species, only small numbers have been identified from the oilseed crop plant Brassica napus (canola) – especially in seeds.

Results

Using next-generation sequencing technologies, we performed a comprehensive analysis of miRNAs during seed maturation at 9 time points from 10 days after flowering (DAF) to 50 DAF using whole seeds and included separate analyses of radicle, hypocotyl, cotyledon, embryo, endosperm and seed coat tissues at 4 selected time points. We identified more than 500 conserved miRNA or variant unique sequences with >300 sequence reads and also found 10 novel miRNAs. Only 27 of the conserved miRNA sequences had been previously identified in B. napus (miRBase Release 18). More than 180 MIRNA loci were identified/annotated using the B. rapa genome as a surrogate for the B.napus A genome. Numerous miRNAs were expressed in a stage- or tissue-specific manner suggesting that they have specific functions related to the fine tuning of transcript abundance during seed development. miRNA targets in B. napus were predicted and their expression patterns profiled using microarray analyses. Global correlation analysis of the expression patterns of miRNAs and their targets revealed complex miRNA-target gene regulatory networks during seed development. The miR156 family was the most abundant and the majority of the family members were primarily expressed in the embryo.

Conclusions

Large numbers of miRNAs with diverse expression patterns, multiple-targeting and co-targeting of many miRNAs, and complex relationships between expression of miRNAs and targets were identified in this study. Several key miRNA-target expression patterns were identified and new roles of miRNAs in regulating seed development are suggested. miR156, miR159, miR172, miR167, miR158 and miR166 are the major contributors to the network controlling seed development and maturation through their pivotal roles in plant development. miR156 may regulate the developmental transition to germination.

Role of RNA interference (RNAi) in the Moss Physcomitrella patens

RNA interference (RNAi) is a mechanism that regulates genes by either transcriptional (TGS) or posttranscriptional gene silencing (PTGS), required for genome maintenance and proper development of an organism. Small non-coding RNAs are the key players in RNAi and have been intensively studied in eukaryotes. In plants, several classes of small RNAs with specific sizes and dedicated functions have evolved. The major classes of small RNAs include microRNAs (miRNAs) and small interfering RNAs (siRNAs), which differ in their biogenesis. miRNAs are synthesized from a short hairpin structure while siRNAs are derived from long double-stranded RNAs (dsRNA). Both miRNA and siRNAs control the expression of cognate target RNAs by binding to reverse complementary sequences mediating cleavage or translational inhibition of the target RNA. They also act on the DNA and cause epigenetic changes such as DNA methylation and histone modifications. In the last years, the analysis of plant RNAi pathways was extended to the bryophyte Physcomitrella patens, a non-flowering, non-vascular ancient land plant that diverged from the lineage of seed plants approximately 450 million years ago. Based on a number of characteristic features and its phylogenetic key position in land plant evolution P. patens emerged as a plant model species to address basic as well as applied topics in plant biology. Here we summarize the current knowledge on the role of RNAi in P. patens that shows functional overlap with RNAi pathways from seed plants, and also unique features specific to this species.

Heterochronic genes in plant evolution and development

Evolution of morphology includes evolutionary shifts of developmental processes in space or in time. Heterochronic evolution is defined as a temporal shift. The concept of heterochrony has been very rewarding to investigators of both animal and plant developmental evolution, because it has strong explanatory power when trying to understand morphological diversity. While for animals, extensive literature on heterochrony developed along with the field of evolution of development, in plants the concept has been applied less often and is less elaborately developed. Yet novel genetic findings highlight heterochrony as a developmental and evolutionary process in plants. Similar to what has been found for the worm Caenorhabditis, a heterochronic gene pathway controlling developmental timing has been elucidated in flowering plants. Two antagonistic microRNA's miR156 and miR172 target two gene families of transcription factors, SQUAMOSA PROMOTOR BINDING PROTEIN-LIKE and APETALA2-like, respectively. Here, we propose that this finding now allows the molecular investigation of cases of heterochronic evolution in plants. We illustrate this point by examining microRNA expression patterns in the Antirrhinum majus incomposita and choripetala heterochronic mutants. Some of the more beautiful putative cases of heterochronic evolution can be found outside flowering plants, but little is known about the extent of conservation of this flowering plant pathway in other land plants. We show that the expression of an APETALA2-like gene decreases with age in a fern species. This contributes to the idea that ferns share some heterochronic gene functions with flowering plants.

Bioinformatics prediction of miRNAs in the Prunus persica genome with validation of their precise sequences by miR-RACE

We predicted 262 potential MicroRNAs (miRNAs) belonging to 70 miRNA families from the peach (Prunus persica) genome and two specific 5' and 3' miRNA rapid amplification of cDNA ends (miR-RACE) PCR reactions and sequence-directed cloning were employed to accurately validate 61 unique P. persica miRNAs (Ppe-miRNAs) sequences belonging to 61 families comprising 97 Ppe-miRNAs. Validation of the termini nucleotides in particular can define the real sequences of the Ppe-miRNAs on peach genome. Comparison between predicted and validated Ppe-miRNAs through alignment revealed that 43 unique orthologous sequences were identical, while the remaining 18 exhibited some divergences at their termini nucleotides. Quantitative real-time polymerase chain reaction (qRT-PCR) was further employed to analyze the expression of all the 61 miRNAs and 10 putative targets of 8 randomly selected Ppe-miRNAs in peach leaves, flowers and fruits at different stages of development, where both the miRNAs and the putative target genes showed tissue-specific expression.

Asymmetric purine-pyrimidine distribution in cellular small RNA population of papaya

Background

The small RNAs (sRNA) are a regulatory class of RNA mainly represented by the 21 and 24-nucleotide size classes. The cellular sRNAs are processed by RNase III family enzyme dicer (Dicer like in plant) from a self-complementary hairpin loop or other type of RNA duplexes. The papaya genome has been sequenced, but its microRNAs and other regulatory RNAs are yet to be analyzed.

Results

We analyzed the genomic features of the papaya sRNA population from three sRNA deep sequencing libraries made from leaves, flowers, and leaves infected with Papaya Ringspot Virus (PRSV). We also used the deep sequencing data to annotate the micro RNA (miRNA) in papaya. We identified 60 miRNAs, 24 of which were conserved in other species, and 36 of which were novel miRNAs specific to papaya. In contrast to the Chargaff’s purine-pyrimidine equilibrium, cellular sRNA was significantly biased towards a purine rich population. Of the two purine bases, higher frequency of adenine was present in 23nt or longer sRNAs, while 22nt or shorter sRNAs were over represented by guanine bases. However, this bias was not observed in the annotated miRNAs in plants. The 21nt species were expressed from fewer loci but expressed at higher levels relative to the 24nt species. The highly expressed 21nt species were clustered in a few isolated locations of the genome. The PRSV infected leaves showed higher accumulation of 21 and 22nt sRNA compared to uninfected leaves. We observed higher accumulation of miRNA* of seven annotated miRNAs in virus-infected tissue, indicating the potential function of miRNA* under stressed conditions.

Conclusions

We have identified 60 miRNAs in papaya. Our study revealed the asymmetric purine-pyrimidine distribution in cellular sRNA population. The 21nt species of sRNAs have higher expression levels than 24nt sRNA. The miRNA* of some miRNAs shows higher accumulation in PRSV infected tissues, suggesting that these strands are not totally functionally redundant. The findings open a new avenue for further investigation of the sRNA silencing pathway in plants.

Genome-wide identification of miRNAs responsive to drought in peach (Prunus persica) by high-throughput deep sequencing

Peach (Prunus persica L.) is one of the most important worldwide fresh fruits. Since fruit growth largely depends on adequate water supply, drought stress is considered as the most important abiotic stress limiting fleshy fruit production and quality in peach. Plant responses to drought stress are regulated both at transcriptional and post-transcriptional level. As post-transcriptional gene regulators, miRNAs (miRNAs) are small (19–25 nucleotides in length), endogenous, non-coding RNAs. Recent studies indicate that miRNAs are involved in plant responses to drought. Therefore, Illumina deep sequencing technology was used for genome-wide identification of miRNAs and their expression profile in response to drought in peach. In this study, four sRNA libraries were constructed from leaf control (LC), leaf stress (LS), root control (RC) and root stress (RS) samples. We identified a total of 531, 471, 535 and 487 known mature miRNAs in LC, LS, RC and RS libraries, respectively. The expression level of 262 (104 up-regulated, 158 down-regulated) of the 453 miRNAs changed significantly in leaf tissue, whereas 368 (221 up-regulated, 147 down-regulated) of the 465 miRNAs had expression levels that changed significantly in root tissue upon drought stress. Additionally, a total of 197, 221, 238 and 265 novel miRNA precursor candidates were identified from LC, LS, RC and RS libraries, respectively. Target transcripts (137 for LC, 133 for LS, 148 for RC and 153 for RS) generated significant Gene Ontology (GO) terms related to DNA binding and catalytic activites. Genome-wide miRNA expression analysis of peach by deep sequencing approach helped to expand our understanding of miRNA function in response to drought stress in peach and Rosaceae. A set of differentially expressed miRNAs could pave the way for developing new strategies to alleviate the adverse effects of drought stress on plant growth and development.

miR156 and miR390 regulate tasiRNA accumulation and developmental timing in Physcomitrella patens

microRNA156 (miR156) affects developmental timing in flowering plants. miR156 and its target relationships with members of the SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) gene family appear universally conserved in land plants, but the specific functions of miR156 outside of flowering plants are unknown. We find that miR156 promotes a developmental change from young filamentous protonemata to leafy gametophores in the moss Physcomitrella patens, opposite to its role as an inhibitor of development in flowering plants. P. patens miR156 also influences accumulation of trans-acting small interfering RNAs (tasiRNAs) dependent upon a second ancient microRNA, miR390. Both miR156 and miR390 directly target a single major tasiRNA primary transcript. Inhibition of miR156 function causes increased miR390-triggered tasiRNA accumulation and decreased accumulation of tasiRNA targets. Overexpression of miR390 also caused a slower formation of gametophores, elevated miR390-triggered tasiRNA accumulation, and reduced level of tasiRNA targets. We conclude that a gene regulatory network controlled by miR156, miR390, and their targets controls developmental change in P. patens. The broad outlines and regulatory logic of this network are conserved in flowering plants, albeit with some modifications. Partially conserved small RNA networks thus influence developmental timing in plants with radically different body plans.

Identification and profiling of novel microRNAs in the Brassica rapa genome based on small RNA deep sequencing

BACKGROUND

MicroRNAs (miRNAs) are one of the functional non-coding small RNAs involved in the epigenetic control of the plant genome. Although plants contain both evolutionary conserved miRNAs and species-specific miRNAs within their genomes, computational methods often only identify evolutionary conserved miRNAs. The recent sequencing of the Brassica rapa genome enables us to identify miRNAs and their putative target genes. In this study, we sought to provide a more comprehensive prediction of B. rapa miRNAs based on high throughput small RNA deep sequencing.

RESULTS

We sequenced small RNAs from five types of tissue: seedlings, roots, petioles, leaves, and flowers. By analyzing 2.75 million unique reads that mapped to the B. rapa genome, we identified 216 novel and 196 conserved miRNAs that were predicted to target approximately 20% of the genome's protein coding genes. Quantitative analysis of miRNAs from the five types of tissue revealed that novel miRNAs were expressed in diverse tissues but their expression levels were lower than those of the conserved miRNAs. Comparative analysis of the miRNAs between the B. rapa and Arabidopsis thaliana genomes demonstrated that redundant copies of conserved miRNAs in the B. rapa genome may have been deleted after whole genome triplication. Novel miRNA members seemed to have spontaneously arisen from the B. rapa and A. thaliana genomes, suggesting the species-specific expansion of miRNAs. We have made this data publicly available in a miRNA database of B. rapa called BraMRs. The database allows the user to retrieve miRNA sequences, their expression profiles, and a description of their target genes from the five tissue types investigated here.

CONCLUSIONS

This is the first report to identify novel miRNAs from Brassica crops using genome-wide high throughput techniques. The combination of computational methods and small RNA deep sequencing provides robust predictions of miRNAs in the genome. The finding of numerous novel miRNAs, many with few target genes and low expression levels, suggests the rapid evolution of miRNA genes. The development of a miRNA database, BraMRs, enables us to integrate miRNA identification, target prediction, and functional annotation of target genes. BraMRs will represent a valuable public resource with which to study the epigenetic control of B. rapa and other closely related Brassica species. The database is available at the following link: http://bramrs.rna.kr [1].

Arabidopsis ribosomal proteins control developmental programs through translational regulation of auxin response factors

Upstream ORFs are elements found in the 5'-leader sequences of specific mRNAs that modulate the translation of downstream ORFs encoding major gene products. In Arabidopsis, the translational control of auxin response factors (ARFs) by upstream ORFs has been proposed as a regulatory mechanism required to respond properly to complex auxin-signaling inputs. In this study, we identify and characterize the aberrant auxin responses in specific ribosomal protein mutants in which multiple ARF transcription factors are simultaneously repressed at the translational level. This characteristic lends itself to the use of these mutants as genetic tools to bypass the genetic redundancy among members of the ARF family in Arabidopsis. Using this approach, we were able to assign unique functions for ARF2, ARF3, and ARF6 in plant development.

Global alteration of microRNAs and transposon-derived small RNAs in cotton (Gossypium hirsutum) during Cotton leafroll dwarf polerovirus (CLRDV) infection

Small RNAs (sRNAs) are a class of non-coding RNAs ranging from 20- to 40-nucleotides (nts) that are present in most eukaryotic organisms. In plants, sRNAs are involved in the regulation of development, the maintenance of genome stability and the antiviral response. Viruses, however, can interfere with and exploit the silencing-based regulatory networks, causing the deregulation of sRNAs, including small interfering RNAs (siRNAs) and microRNAs (miRNAs). To understand the impact of viral infection on the plant sRNA pathway, we deep sequenced the sRNAs in cotton leaves infected with Cotton leafroll dwarf virus (CLRDV), which is a member of the economically important virus family Luteoviridae. A total of 60 putative conserved cotton miRNAs were identified, including 19 new miRNA families that had not been previously described in cotton. Some of these miRNAs were clearly misregulated during viral infection, and their possible role in symptom development and disease progression is discussed. Furthermore, we found that the 24-nt heterochromatin-associated siRNAs were quantitatively and qualitatively altered in the infected plant, leading to the reactivation of at least one cotton transposable element. This is the first study to explore the global alterations of sRNAs in virus-infected cotton plants. Our results indicate that some CLRDV-induced symptoms may be correlated with the deregulation of miRNA and/or epigenetic networks.

Overexpression of gma-MIR394a confers tolerance to drought in transgenic Arabidopsis thaliana

MicroRNAs, key posttranscriptional regulators of eukaryotic gene expression, play important roles in plant development and response to stress. In this study, a soybean gma-MIR394a gene was functionally characterized, especially with regard to its role in drought stress resistance. Expression analysis revealed that gma-MIR394a was expressed differentially in various soybean tissues and was induced by drought, high salinity, low temperature stress, and abscisic acid treatment in leaves. One target gene of gma-miR394a, Glyma08g11030, was predicted and verified using a modified 5' RLM-RACE (RNA ligase-mediated rapid amplification of 5' cDNA ends) assay. Overexpression of gma-MIR394a resulted in plants with lowered leaf water loss and enhanced drought tolerance. Furthermore, overexpression of gma-MIR394a in Arabidopsis reduced the transcript of an F-box gene (At1g27340) containing a miR394 complementary target site. These results suggest that the gma-MIR394a gene functions in positive modulation of drought stress tolerance and has potential applications in molecular breeding to enhance drought tolerance in crops.

Identification of conserved and novel microRNAs from Liriodendron chinense floral tissues

Background

Liriodendron chinense (L. chinense) is an endangered basal angiosperm plant in China because of its low reproductive efficiency. Recently, miRNAs have obtained great attention because they can play important roles. Through high throughput sequencing technique, large amount of miRNAs were identified from different plant species. But there were few studies about the miRNAs in the basal angiosperms especially in the sexual reproduction process.

Results

Deep sequencing technology was applied to discover miRNAs in L. chinense flowers at different stages. After bioinformatic analysis, 496 putative conserved miRNAs representing 97 families and 2 novel miRNAs were found. Among them, one is previously regarded as gymnosperm specific. Their expressions were further validated by Real-time PCR for 13 selected miRNAs. Putative targeting genes were predicted and categorized with gene ontology (GO) analysis. About ten percents of the targets are involved in the reproduction process. Further expressional analysis showed that many of these miRNAs were highly related to the reproductive growth.

Conclusions

This is the first comprehensive identification of conserved and novel miRNAs in L. chinense. The data presented here might not only help to fill the gap of miRNA registered about basal angiosperm plants but also contribute to understanding the evolution of miRNAs. The differential expression of some of the miRNAs and the prediction of their target genes are also helpful in understanding the regulation of L. chinense sexual reproduction.

Evolution of the ARF gene family in land plants: old domains, new tricks

Auxin response factors (ARF) are key players in plant development. They mediate the cellular response to the plant hormone auxin by activating or repressing the expression of downstream developmental genes. The pivotal activation function of ARF proteins is enabled by their four-domain architecture, which includes both DNA-binding and protein dimerization motifs. To determine the evolutionary origin of this characteristic architecture, we built a comprehensive data set of 224 ARF-related protein sequences that represents all major living divisions of land plants, except hornworts. We found that ARFs are split into three subfamilies that could be traced back to the origin of the land plants. We also show that repeated events of extensive gene duplication contributed to the expansion of those three original subfamilies. Further examination of our data set uncovered a broad diversity in the structure of ARF transcripts and allowed us to identify an additional conserved motif in ARF proteins. We found that additional structural diversity in ARF proteins is mainly generated by two mechanisms: genomic truncation and alternative splicing. We propose that the loss of domains from the canonical, four-domain ARF structure has promoted functional shifts within the ARF family by disrupting either dimerization or DNA-binding capabilities. For instance, the loss of dimerization domains in some ARFs from moss and spikemoss genomes leads to proteins that are reminiscent of Aux/IAA proteins, possibly providing a clue on the evolution of these modulators of ARF function. We also assessed the functional impact of alternative splicing in the case of ARF4, for which we have identified a novel isoform in Arabidopsis thaliana. Genetic analysis showed that these two transcripts exhibit markedly different developmental roles in A. thaliana. Gene duplications, domain rearrangement, and post-transcriptional regulation have thus enabled a subtle control of auxin signaling through ARF proteins that may have contributed to the critical importance of these regulators in plant development and evolution.