Computational identification of microRNAs in apple expressed sequence tags and validation of their precise sequences by miR-RACE

Small RNA and Degradome Sequencing Reveal Roles of miRNAs in the Petal Color Fading of Malus Crabapple

The Malus crabapple is an important woody ornamental plant. The fading of petals during its development significantly affects their ornamental value. Petal color is related to anthocyanin content and miRNAs play an important role in the post-transcriptional regulation of anthocyanin synthesis. However, the mechanisms underlying miRNA regulation of petal fading have rarely been studied. Transcriptome and small RNA sequencing of petals from the blooming phases of Malus. 'Indian Summer' varieties S1 (small bud), S2 (initial-flowering), and S3 (late-flowering) allowed us to identify 230 known miRNAs and 17 novel miRNAs, including 52 differentially expressed miRNAs which targeted 494 genes and formed 823 miRNA-target pairs. Based on the target gene annotation results, miRNA-target pairs were screened that may be involved in the fading process of Malus crabapple petals through three different pathways: anthocyanin synthesis, transport, and degradation, involving mcr-miR858-MYB1\MYB5 and mcr-miR396-McCHI inhibiting anthocyanin synthesis; mcr-miR167, mcr-miR390, mcr-miR535, and mcr-miR858 inhibiting anthocyanin transport from the cytoplasm to the vacuole by targeting ABC transporter genes (ABCB, ABCC, ABCD, and ABCG); and mcr-miR398 targeting the superoxide dismutase genes (CZSOD2 and CCS) to accelerate anthocyanin degradation. These findings offer a novel approach to understanding the mechanism of petal fading and serve as a reference for other plants with floral fading.

Integrated miRNA, target mRNA, and metabolome profiling of Tinospora cordifolia with reference to berberine biosynthesis

MicroRNAs play a central role in gene regulation and emerge as novel targets for secondary metabolites improvement in plants. The crops thus can be improved through knowledge obtained by the study of miRNAs because of their conserved nature in gene regulation. The present study has been carried out on Tinospora cordifolia (T. cordifolia), because of its illimitable application for the treatment of various diseases. This plant has tremendous medicinal properties, yet unexplored at the molecular level, and has not received much recognition in the scientific field. Thus, here computational analysis was performed to identify T. cordifolia miRNAs using EST database. Using these miRNAs, we predicted their targets which were found to be associated with the regulation of diverse gene networks including 433 berberine biosynthesis genes in T. cordifolia. Further, selected miRNAs were validated and their expression was detected in different T. cordifolia tissues followed by expression analysis of their target mRNAs. These data were then compared with the metabolic profile of T. cordifolia with an emphasis on therapeutically important compound berberine. In this study, we did simultaneous miRNA/target gene expression and metabolome analysis which opens a new way for initiating new proposition and prioritization of miRNAs/genes/metabolites for targeted follow‑up metabolic engineering experimentations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03342-9.

Malus hupehensis miR168 Targets to ARGONAUTE1 and Contributes to the Resistance against Botryosphaeria dothidea Infection by Altering Defense Responses

MicroRNA (miRNA)-mediated post-transcriptional regulation plays a fundamental role in various plant physiological processes, including responses to pathogens. MicroRNA168 has been implicated as an essential factor of miRNA pathways by targeting ARGONAUTE1 (AGO1), the core component of the RNA-induced silencing complex (RISC). A fluctuation in AGO1 expression influences various plant-pathogen interactions, and the homeostasis of AGO1 and miR168 accumulation is maintained by a complicated feedback regulatory loop. In this study, the connection between miR168 and the resistance of Malus hupehensis to Botryosphaeria dothidea is revealed. The induction of both the mature miR168 and its precursor in plants subjected to B. dothidea infection indicate the transcriptional activation of MIR168a. MIR168a promoter analysis demonstrates that the promoter can be activated by B. dothidea and salicylic acid (SA). However, the direct target of miR168, M. hupehensis ARGONAUTE1 (MhAGO1), is shown to be induced under the infection. Expression and transcription activity analysis demonstrate the transcriptional activation and the post-transcriptional suppression of MhAGO1 in response to B. dothidea infection. By inhibiting reactive oxygen species (ROS) production and enhancing SA-mediated defense responses, miR168a delays the symptom development of leaves inoculated with B. dothidea and impedes the pathogen growth, while MhAGO1 is found to have the opposite effects. Collectively, these findings suggest that the expression of miR168 and MhAGO1 in M. hupehensis in response to B. dothidea infection is regulated by a complicated mechanism. Targeting to MhAGO1, a negative regulator, miR168 plays a positive role in the resistance by alterations in diverse defense responses.

Dme-miR-314-3p modulation in Cr(VI) exposed Drosophila affects DNA damage repair by targeting mus309

microRNAs (miRNAs) as one of the major epigenetic modulators negatively regulate mRNAs at post transcriptional level. It was therefore hypothesized that modulation of miRNAs by hexavalent Chromium [Cr(VI)], a priority environmental chemical, can affect DNA damage. In a genetically tractable model, Drosophila melanogaster, role of maximally up-regulated miRNA, dme-miR-314-3p, on DNA damage was examined by exposing the third instar larvae to 5.0-20.0 μg/ml Cr(VI) for 24 and 48 h. mus309, a Drosophila homologue of human Bloom's syndrome and predicted as one of the potential targets of this miRNA, was confirmed as its target by 5'RLM-RACE assay. A significant down-regulation of mus309 was observed in dme-miR-314-3p overexpression strain (myo-gal4>UAS-miR-314-3p) as compared with that in parental strains (myo-gal4 and UAS-miR-314-3p) and in w(1118). A significant increase in DNA damage including double strand breaks generation was observed in exposed myo-gal4>UAS-miR-314 and mus309 mutants as compared with that in parental strain and in unexposed control. A significant down-regulation of cell cycle regulation genes (CycA, CycB and cdc2) was observed in these exposed genotypes. Collectively, the study demonstrates that dme-miR-314-3p can mediate the downregulation of repair deficient gene mus309 leading to increased DNA damage and cell cycle arrest in exposed organism which may affect Cr(VI) mediated carcinogenesis.

Selection of Reference Genes for Normalization of MicroRNA Expression by RT-qPCR in Sugarcane Buds under Cold Stress

Sugarcane, accounting for 80% of world's sugar, originates in the tropics but is cultivated mainly in the subtropics. Therefore, chilling injury frequently occurs and results in serious losses. Recent studies in various plant species have established microRNAs as key elements in the post-transcriptional regulation of response to biotic and abiotic stresses including cold stress. Though, its accuracy is largely influenced by the use of reference gene for normalization, quantitative PCR is undoubtedly a popular method used for identification of microRNAs. For identifying the most suitable reference genes for normalizing miRNAs expression in sugarcane under cold stress, 13 candidates among 17 were investigated using four algorithms: geNorm, NormFinder, deltaCt, and Bestkeeper, and four candidates were excluded because of unsatisfactory efficiency and specificity. Verification was carried out using cold-related genes miR319 and miR393 in cold-tolerant and sensitive cultivars. The results suggested that miR171/18S rRNA and miR171/miR5059 were the best reference gene sets for normalization for miRNA RT-qPCR, followed by the single miR171 and 18S rRNA. These results can aid research on miRNA responses during sugarcane stress, and the development of sugarcane tolerant to cold stress. This study is the first report concerning the reference gene selection of miRNA RT-qPCR in sugarcane.

Determination of the precise sequences of computationally predicted miRNAs in Citrus reticulata by miR-RACE and characterization of the related target genes using RLM-RACE

MicroRNAs play vital roles in various biological and metabolic processes by regulating the expression of their target genes in model plants. Since there are limited reports on miRNAs in Citrus reticulata (Crt-miRNAs), the determination of precise sequences of miRNAs is essential to further analyze the functions of miRNAs in Citrus reticulata. Here, miR-RACE, a recently developed technique for determination of the potential miRNAs computationally, was employed to identify the precise sequences of Crt-miRNAs. Tissue- and development-specific expression of nine miRNAs were identified by quantitative RT-PCR in the leaves, stems, flowers and fruits Subsequently, 10 potential target genes were predicated for the eight Crt-miRNAs, most of which were transcription factors and disease resistance proteins. Four target genes were experimentally validated by Poly (A) polymerase-mediated 3′ rapid amplification of cDNA ends and RNA ligase-mediated 5′ rapid amplification of cDNA ends (PPM-RACE and RLM-RACE). Our findings showed that regulatory miRNAs in C. reticulata may play a key role in regulating growth, development, and response to disease. Future work is required to study the functions of miRNAs and their targets of C. reticulata.

A method for validating microRNAs in plants by miR-RACE

miRNA prediction algorithms often fail to predict the accurate location of the mature miRNA in a precursor sequence with nucleotide-level precision. miRNAs-rapid amplification of cDNA ends (miR-RACE) is an efficient method to determine the precise sequences of computationally predicted microRNAs (miRNAs). miR-RACE includes the following steps: miRNA-enriched library preparation, two specific 5'- and 3'-miRNA RACE (miR-RACE) PCR reactions, and sequence-directed cloning. The most challenging step is the two gene-specific primers designed for the two RACE reactions. The miR-RACE protocol is rapid and can be executed and completed in 2-3 days.

Advances in identification and validation of plant microRNAs and their target genes

Developments in the field of molecular biology and genetics, such as microarray, gene transfer and discovery of small regulatory RNAs, have led to significant advances in plant biotechnology. Among the small RNAs, microRNAs (miRNAs) have elicited much interest as key post-transcriptional regulators in eukaryotic gene expression. Advances in genome and transcriptome sequencing of plants have facilitated the generation of a huge wealth of sequence information that can find much use in the discovery of novel miRNAs and their target genes. In this review, we present an overview of the developments in the strategies and methods used to identify and study miRNAs, their target genes and the mechanisms by which these miRNAs interact with their target genes since the discovery of the first miRNA. The approaches discussed include both reverse and forward genetics. We observed that despite the availability of advanced methods, certain limitations ranging from the cost of materials, equipment and personnel to the availability of genome sequences for many plant species present a number of challenges for the development and utilization of modern scientific methods for the elucidation and development of miRNAs in many important plant species.

MicroRNAs in fruit trees: discovery, diversity and future research directions

Since the first description of microRNAs (miRNAs) 20 years ago, the number of miRNAs identified in different eukaryotic organisms has exploded, largely due to the recent advances in DNA sequencing technologies. Functional studies, mostly from model species, have revealed that miRNAs are major post-transcriptional regulators of gene expression in eukaryotes. In plants, they are implicated in fundamental biological processes, from plant development and morphogenesis, to regulation of plant pathogen and abiotic stress responses. Although a substantial number of miRNAs have been identified in fruit trees to date, their functions remain largely uncharacterised. The present review aims to summarise the progress made in miRNA research in fruit trees, focusing specifically on the economically important species Prunus persica, Malus domestica, Citrus spp, and Vitis vinifera. We also discuss future miRNA research prospects in these plants and highlight potential applications of miRNAs in the on-going improvement of fruit trees.

Identification of miRNAs and their targets in tea (Camellia sinensis)

MicroRNAs (miRNAs) are endogenous small RNAs playing a crucial role in plant growth and development, as well as stress responses. Among them, some are highly evolutionally conserved in the plant kingdom, this provide a powerful strategy for identifying miRNAs in a new species. Tea (Camellia sinensis) is one of the most important commercial beverage crops in the world, but only a limited number of miRNAs have been identified. In the present study, a total of 14 new C. sinensis miRNAs were identified by expressed sequence tag (EST) analysis from 47452 available C. sinensis ESTs. These miRNAs potentially target 51 mRNAs, which can act as transcription factors, and participate in stress response, transmembrane transport, and signal transduction. Analysis of gene ontology (GO), based on these targets, suggested that 37 biological processes were involved, such as oxidation-reduction process, stress response, and transport. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis inferred that the identified miRNAs took part in 13 metabolic networks. Our study will help further understanding of the essential roles of miRNAs in C. sinensis growth and development, and stress response.

Genomic resources in fruit plants: an assessment of current status

The availability of many genomic resources such as genome sequences, functional genomics resources including microarrays and RNA-seq, sufficient numbers of molecular markers, express sequence tags (ESTs) and high-density genetic maps is causing a rapid acceleration of genetics and genomic research of many fruit plants. This is leading to an increase in our knowledge of the genes that are linked to many horticultural and agronomically important traits. Recently, some progress has also been made on the identification and functional analysis of miRNAs in some fruit plants. This is one of the most active research fields in plant sciences. The last decade has witnessed development of genomic resources in many fruit plants such as apple, banana, citrus, grapes, papaya, pears, strawberry etc.; however, many of them are still not being exploited. Furthermore, owing to lack of resources, infrastructure and research facilities in many lesser-developed countries, development of genomic resources in many underutilized or less-studied fruit crops, which grow in these countries, is limited. Thus, research emphasis should be given to those fruit crops for which genomic resources are relatively scarce. The development of genomic databases of these less-studied fruit crops will enable biotechnologists to identify target genes that underlie key horticultural and agronomical traits. This review presents an overview of the current status of the development of genomic resources in fruit plants with the main emphasis being on genome sequencing, EST resources, functional genomics resources including microarray and RNA-seq, identification of quantitative trait loci and construction of genetic maps as well as efforts made on the identification and functional analysis of miRNAs in fruit plants.

Computational identification and characterization of conserved miRNAs and their target genes in garlic (Allium sativum L.) expressed sequence tags

The endogenous small non-coding functional microRNAs (miRNAs) are short in size, range from ~21 to 24 nucleotides in length, play a pivotal role in gene expression in plants and animals by silencing genes either by destructing or blocking of translation of homologous mRNA. Although various high-throughput, time consuming and expensive techniques like forward genetics and direct cloning are employed to detect miRNAs in plants but comparative genomics complemented with novel bioinformatic tools pave the way for efficient and cost-effective identification of miRNAs through homologous sequence search with previously known miRNAs. In this study, an attempt was made to identify and characterize conserved miRNAs in garlic expressed sequence tags (ESTs) through computational means. For identification of novel miRNAs in garlic, a total 3227 known mature miRNAs of plant kingdom Viridiplantae were searched for homology against 21,637 EST sequences resulting in identification of 6 potential miRNA candidates belonging to 6 different miRNA families. The psRNATarget server predicted 33 potential target genes and their probable functions for the six identified miRNA families in garlic. Most of the garlic miRNA target genes seem to encode transcription factors as well as genes involved in stress response, metabolism, plant growth and development. The results from the present study will shed more light on the understanding of molecular mechanisms of miRNA in garlic which may aid in the development of novel and precise techniques to understand some post-transcriptional gene silencing mechanism in response to stress tolerance.

Computational identification of conserved microRNAs and their targets from expression sequence tags of blueberry (Vaccinium corybosum)

MicroRNAs (miRNAs) are a class of endogenous, approximately 21nt in length, non-coding RNA, which mediate the expression of target genes primarily at post-transcriptional levels. miRNAs play critical roles in almost all plant cellular and metabolic processes. Although numerous miRNAs have been identified in the plant kingdom, the miRNAs in blueberry, which is an economically important small fruit crop, still remain totally unknown. In this study, we reported a computational identification of miRNAs and their targets in blueberry. By conducting an EST-based comparative genomics approach, 9 potential vco-miRNAs were discovered from 22,402 blueberry ESTs according to a series of filtering criteria, designated as vco-miR156-5p, vco-miR156-3p, vco-miR1436, vco-miR1522, vco-miR4495, vco-miR5120, vco-miR5658, vco-miR5783, and vco-miR5986. Based on sequence complementarity between miRNA and its target transcript, 34 target ESTs from blueberry and 70 targets from other species were identified for the vco-miRNAs. The targets were found to be involved in transcription, RNA splicing and binding, DNA duplication, signal transduction, transport and trafficking, stress response, as well as synthesis and metabolic process. These findings will greatly contribute to future research in regard to functions and regulatory mechanisms of blueberry miRNAs.

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.

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

Characterization and expression profiling of selected microRNAs in tomato (Solanum lycopersicon) 'Jiangshu14'

Presence of selected tomato (Solanum lycopersicon) microRNAs (sly-miRNAs) was validated and their expression profiles established in roots, stems, leaves, flowers and fruits of tomato variety Jiangshu14 by quantitative RT-PCR (qRT-PCR). In addition conservation characteristics these sly-miRNAs were analyzed and target genes predicted bioinformatically. Results indicate that some of these miRNAs are specific to tomato while most are conserved in other plant species. Predicted sly-miRNA targets genes were shown to be targeted by either by a single or more miRNAs and are involved in diverse processes in tomato plant growth and development. All the 36 miRNAs were present in the cDNA of mixed tissues and qRT-PCR revealed that some of these sly-miRNAs are ubiquitous in tomato while others have tissue-specific expression. The experimental validation and expression profiling as well target gene prediction of these miRNAs in tomato as done in this study can add to the knowledge on the important roles played by these sly-miRNAs in the growth and development, environmental stress tolerance as well as pest and disease resistance in tomatoes and related species. In addition these findings broaden the knowledge of small RNA-mediated regulation in S. lycopersicon. It is recommended that experimental validation of the target genes be done so as to give a much more comprehensive information package on these miRNAs in tomato and specifically in the selected variety.

In silico identification and characterization of conserved miRNAs and their target genes in sweet potato (Ipomoea batatas L.) expressed sequence tags (ESTs)

The endogenous small non-coding micro RNAs (miRNAs), which are typically ~21-24 nt nucleotides, play a crucial role in regulating the intrinsic normal growth of cells and development of the plants as well as in maintaining the integrity of genomes. These small non-coding RNAs function as the universal specificity factors in post-transcriptional gene silencing. Discovering miRNAs, identifying their targets, and further inferring miRNA functions is a routine process to understand normal biological processes of miRNAs and their roles in the development of plants. Comparative genomics based approach using expressed sequence tags (EST) and genome survey sequences (GSS) offer a cost-effective platform for identification and characterization of miRNAs and their target genes in plants. Despite the fact that sweet potato (Ipomoea batatas L.) is an important staple food source for poor small farmers throughout the world, the role of miRNA in various developmental processes remains largely unknown. In this paper, we report the computational identification of miRNAs and their target genes in sweet potato from their ESTs. Using comparative genomics-based approach, 8 potential miRNA candidates belonging to miR168, miR2911, and miR156 families were identified from 23 406 ESTs in sweet potato. A total of 42 target genes were predicted and their probable functions were illustrated. Most of the newly identified miRNAs target transcription factors as well as genes involved in plant growth and development, signal transduction, metabolism, defense, and stress response. The identification of miRNAs and their targets is expected to accelerate the pace of miRNA discovery, leading to an improved understanding of the role of miRNA in development and physiology of sweet potato, as well as stress response.

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.

Genome wide identification of chilling responsive microRNAs in Prunus persica

BACKGROUND

MicroRNAs (miRNAs) are small RNAs (sRNAs) approximately 21 nucleotides in length that negatively control gene expression by cleaving or inhibiting the translation of target gene transcripts. Within this context, miRNAs and siRNAs are coming to the forefront as molecular mediators of gene regulation in plant responses to annual temperature cycling and cold stress. For this reason, we chose to identify and characterize the conserved and non-conserved miRNA component of peach (Prunus persica (L.) Batsch) focusing our efforts on both the recently released whole genome sequence of peach and sRNA transcriptome sequences from two tissues representing non-dormant leaves and dormant leaf buds. Conserved and non-conserved miRNAs, and their targets were identified. These sRNA resources were used to identify cold-responsive miRNAs whose gene targets co-localize with previously described QTLs for chilling requirement (CR).

RESULTS

Analysis of 21 million peach sRNA reads allowed us to identify 157 and 230 conserved and non-conserved miRNA sequences. Among the non-conserved miRNAs, we identified 205 that seem to be specific to peach. Comparative genome analysis between peach and Arabidopsis showed that conserved miRNA families, with the exception of miR5021, are similar in size. Sixteen of these conserved miRNA families are deeply rooted in land plant phylogeny as they are present in mosses and/or lycophytes. Within the other conserved miRNA families, five families (miR1446, miR473, miR479, miR3629, and miR3627) were reported only in tree species (Populustrichocarpa, Citrus trifolia, and Prunus persica). Expression analysis identified several up-regulated or down-regulated miRNAs in winter buds versus young leaves. A search of the peach proteome allowed the prediction of target genes for most of the conserved miRNAs and a large fraction of non-conserved miRNAs. A fraction of predicted targets in peach have not been previously reported in other species. Several conserved and non-conserved miRNAs and miRNA-regulated genes co-localize with Quantitative Trait Loci (QTLs) for chilling requirement (CR-QTL) and bloom date (BD-QTL).

CONCLUSIONS

In this work, we identified a large set of conserved and non-conserved miRNAs and describe their evolutionary footprint in angiosperm lineages. Several of these miRNAs were induced in winter buds and co-localized with QTLs for chilling requirement and bloom date thus making their gene targets potential candidates for mediating plant responses to cold stress. Several peach homologs of genes participating in the regulation of vernalization in Arabidopsis were identified as differentially expressed miRNAs targets, potentially linking these gene activities to cold responses in peach dormant buds. The non-conserved miRNAs may regulate cellular, physiological or developmental processes specific to peach and/or other tree species.

Identification and validation of potential conserved microRNAs and their targets in peach (Prunus persica)

MicroRNAs are a class of small, endogenous, non-coding RNA molecules that negatively regulate gene expression at the transcriptional or the post-transcriptional level. Although a large number of miRNAs have been identified in many plant species, especially from model plants and crops, they remain largely unknown in peach. In this study, 110 potential miRNAs belonging to 37 families were identified using computational methods. A total of 43 potential targets were found for 21 families based on near-perfect or perfect complementarity between the plant miRNA and the target sequences. A majority of the targets were transcription factors which play important roles in peach development. qRT-PCR analysis of RNA samples prepared from different peach tissues for 25 miRNA families revealed that miRNAs were differentially expressed in different tissues. Furthermore, two target genes were experimentally verified by detection of the miRNA-mediated mRNA cleavage sites in peach using RNA ligase-mediated 5' rapid amplification of cDNA ends (RLM-RACE). Finally, we studied the expression pattern of the two target genes in three different tissues of peach to further understand the mechanism of the interaction between miRNAs and their target genes.

MicroRNAs and their diverse functions in plants

microRNAs (miRNAs) are an extensive class of newly identified small RNAs, which regulate gene expression at the post-transcriptional level by mRNA cleavage or translation inhibition. Currently, there are 3,070 miRNAs deposited in the public available miRNA database; these miRNAs were obtained from 43 plant species using both computational (comparative genomics) and experimental (direct cloning and deep sequencing) approaches. Like other signaling molecules, plant miRNAs can also be moved from one tissue to another through the vascular system. These mobile miRNAs may play an important role in plant nutrient homeostasis and response to environmental biotic and abiotic stresses. In addition, miRNAs also control a wide range of biological and metabolic processes, including developmental timing, tissue-specific development, and stem cell maintenance and differentiation. Currently, a majority of plant miRNA-related researches are purely descriptive, and provide no further detailed mechanistic insight into miRNA-mediated gene regulation and other functions. To better understand the function and regulatory mechanisms of plant miRNAs, more strategies need to be employed to investigate the functions of miRNAs and their associated signaling pathways and gene networks. Elucidating the evolutionary mechanism of miRNAs is also important. It is possible to develop a novel miRNA-based biotechnology for improving plant yield, quality and tolerance to environmental biotic and abiotic stresses besides focusing on basic genetic studies.

MicroRNAs in trees

MicroRNAs (miRNAs) are 20-24 nucleotide long molecules processed from a specific class of RNA polymerase II transcripts that mainly regulate the stability of mRNAs containing a complementary sequence by targeted degradation in plants. Many features of tree biology are regulated by miRNAs affecting development, metabolism, adaptation and evolution. MiRNAs may be modified and harnessed for controlled suppression of specific genes to learn about gene function, or for practical applications through genetic engineering. Modified (artificial) miRNAs act as dominant suppressors and are particularly useful in tree genetics because they bypass the generations of inbreeding needed for fixation of recessive mutations. The purpose of this review is to summarize the current status of information on miRNAs in trees and to guide future studies on the role of miRNAs in the biology of woody perennials and to illustrate their utility in directed genetic modification of trees.

Identification and comparative profiling of miRNAs in an early flowering mutant of trifoliate orange and its wild type by genome-wide deep sequencing

MicroRNAs (miRNAs) are a new class of small, endogenous RNAs that play a regulatory role in various biological and metabolic processes by negatively affecting gene expression at the post-transcriptional level. While the number of known Arabidopsis and rice miRNAs is continuously increasing, information regarding miRNAs from woody plants such as citrus remains limited. Solexa sequencing was performed at different developmental stages on both an early flowering mutant of trifoliate orange (precocious trifoliate orange, Poncirus trifoliata L. Raf.) and its wild-type in this study, resulting in the obtainment of 141 known miRNAs belonging to 99 families and 75 novel miRNAs in four libraries. A total of 317 potential target genes were predicted based on the 51 novel miRNAs families, GO and KEGG annotation revealed that high ranked miRNA-target genes are those implicated in diverse cellular processes in plants, including development, transcription, protein degradation and cross adaptation. To characterize those miRNAs expressed at the juvenile and adult development stages of the mutant and its wild-type, further analysis on the expression profiles of several miRNAs through real-time PCR was performed. The results revealed that most miRNAs were down-regulated at adult stage compared with juvenile stage for both the mutant and its wild-type. These results indicate that both conserved and novel miRNAs may play important roles in citrus growth and development, stress responses and other physiological processes.

Unique expression, processing regulation, and regulatory network of peach (Prunus persica) miRNAs

BACKGROUND

MicroRNAs (miRNAs) have recently emerged as important gene regulators in plants. MiRNAs and their targets have been extensively studied in Arabidopsis and rice. However, relatively little is known about the characterization of miRNAs and their target genes in peach (Prunus persica), which is a complex crop with unique developmental programs.

RESULTS

We performed small RNA deep sequencing and identified 47 peach-specific and 47 known miRNAs or families with distinct expression patterns. Together, the identified miRNAs targeted 80 genes, many of which have not been reported previously. Like the model plant systems, peach has two of the three conserved trans-acting siRNA biogenesis pathways with similar mechanistic features and target specificity. Unique to peach, three of the miRNAs collectively target 49 MYBs, 19 of which are known to regulate phenylpropanoid metabolism, a key pathway associated with stone hardening and fruit color development, highlighting a critical role of miRNAs in the regulation of peach fruit development and ripening. We also found that the majority of the miRNAs were differentially regulated in different tissues, in part due to differential processing of miRNA precursors. Up to 16% of the peach-specific miRNAs were differentially processed from their precursors in a tissue specific fashion, which has been rarely observed in plant cells. The miRNA precursor processing activity appeared not to be coupled with its transcriptional activity but rather acted independently in peach.

CONCLUSIONS

Collectively, the data characterizes the unique expression pattern and processing regulation of peach miRNAs and demonstrates the presence of a complex, multi-level miRNA regulatory network capable of targeting a wide variety of biological functions, including phenylpropanoid pathways which play a multifaceted spatial-temporal role in peach fruit development.

High-throughput sequencing of small RNAs and analysis of differentially expressed microRNAs associated with pistil development in Japanese apricot

BACKGROUND

MicroRNAs (miRNAs) are a class of endogenous, small, non-coding RNAs that regulate gene expression by mediating gene silencing at transcriptional and post-transcriptional levels in high plants. However, the diversity of miRNAs and their roles in floral development in Japanese apricot (Prunus mume Sieb. et Zucc) remains largely unexplored. Imperfect flowers with pistil abortion seriously decrease production yields. To understand the role of miRNAs in pistil development, pistil development-related miRNAs were identified by Solexa sequencing in Japanese apricot.

RESULTS

Solexa sequencing was used to identify and quantitatively profile small RNAs from perfect and imperfect flower buds of Japanese apricot. A total of 22,561,972 and 24,952,690 reads were sequenced from two small RNA libraries constructed from perfect and imperfect flower buds, respectively. Sixty-one known miRNAs, belonging to 24 families, were identified. Comparative profiling revealed that seven known miRNAs exhibited significant differential expression between perfect and imperfect flower buds. A total of 61 potentially novel miRNAs/new members of known miRNA families were also identified by the presence of mature miRNAs and corresponding miRNA*s in the sRNA libraries. Comparative analysis showed that six potentially novel miRNAs were differentially expressed between perfect and imperfect flower buds. Target predictions of the 13 differentially expressed miRNAs resulted in 212 target genes. Gene ontology (GO) annotation revealed that high-ranking miRNA target genes are those implicated in the developmental process, the regulation of transcription and response to stress.

CONCLUSIONS

This study represents the first comparative identification of miRNAomes between perfect and imperfect Japanese apricot flowers. Seven known miRNAs and six potentially novel miRNAs associated with pistil development were identified, using high-throughput sequencing of small RNAs. The findings, both computationally and experimentally, provide valuable information for further functional characterisation of miRNAs associated with pistil development in plants.

Apple miRNAs and tasiRNAs with novel regulatory networks

BACKGROUND

MicroRNAs (miRNAs) and their regulatory functions have been extensively characterized in model species but whether apple has evolved similar or unique regulatory features remains unknown.

RESULTS

We performed deep small RNA-seq and identified 23 conserved, 10 less-conserved and 42 apple-specific miRNAs or families with distinct expression patterns. The identified miRNAs target 118 genes representing a wide range of enzymatic and regulatory activities. Apple also conserves two TAS gene families with similar but unique trans-acting small interfering RNA (tasiRNA) biogenesis profiles and target specificities. Importantly, we found that miR159, miR828 and miR858 can collectively target up to 81 MYB genes potentially involved in diverse aspects of plant growth and development. These miRNA target sites are differentially conserved among MYBs, which is largely influenced by the location and conservation of the encoded amino acid residues in MYB factors. Finally, we found that 10 of the 19 miR828-targeted MYBs undergo small interfering RNA (siRNA) biogenesis at the 3' cleaved, highly divergent transcript regions, generating over 100 sequence-distinct siRNAs that potentially target over 70 diverse genes as confirmed by degradome analysis.

CONCLUSIONS

Our work identified and characterized apple miRNAs, their expression patterns, targets and regulatory functions. We also discovered that three miRNAs and the ensuing siRNAs exploit both conserved and divergent sequence features of MYB genes to initiate distinct regulatory networks targeting a multitude of genes inside and outside the MYB family.

Identification and characterization of microRNAs and their target genes in Brassica oleracea

The microRNAs are a new class of small non-coding endogenous RNAs with lengths of approximately ~21 nt. MicroRNAs perform their biological function via the degradation of the target mRNAs or by inhibiting protein translation. Until recently, only limited numbers of miRNAs were identified in Brassica oleracea, a vegetable widely cultivated around the world. In present study, 193 potential miRNA candidates were identified from 17 expressed sequence tag (ESTs) and 152 genome survey sequences (GSSs) in B. oleracea. These miRNA candidates were classified into 70 families using a well-defined comparative genome-based computational analysis. Most miRNAs belong to the miRNA169, miR5021, miR156 and miR158 families. Of these, 36 miRNA families are firstly found in Brassica species. Around 1393 B. oleracea genes were predicted as candidate targets of 175 miRNAs. The mutual relationship between miRNAs and the candidate target genes was verified by checking differentially expression levels using quantitative real-time polymerase chain reaction (qRT-PCR) and 5' RLM-RACE analyses. These target genes participate in multiple biological and metabolic processes, including signal transduction, stress response, and plant development. Gene Ontology analysis shows that the 818, 514, and 265 target genes are involved in molecular functions, biological processes, and cellular component respectively. The Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway enrichment analysis suggests that these miRNAs might regulate 186 metabolic pathways, including those of lipid, energy, starch and sucrose, fatty acid and nitrogen.

Identification of microRNAs from Amur grape (Vitis amurensis Rupr.) by deep sequencing and analysis of microRNA variations with bioinformatics

BACKGROUND

MicroRNA (miRNA) is a class of functional non-coding small RNA with 19-25 nucleotides in length while Amur grape (Vitis amurensis Rupr.) is an important wild fruit crop with the strongest cold resistance among the Vitis species, is used as an excellent breeding parent for grapevine, and has elicited growing interest in wine production. To date, there is a relatively large number of grapevine miRNAs (vv-miRNAs) from cultivated grapevine varieties such as Vitis vinifera L. and hybrids of V. vinifera and V. labrusca, but there is no report on miRNAs from Vitis amurensis Rupr, a wild grapevine species.

RESULTS

A small RNA library from Amur grape was constructed and Solexa technology used to perform deep sequencing of the library followed by subsequent bioinformatics analysis to identify new miRNAs. In total, 126 conserved miRNAs belonging to 27 miRNA families were identified, and 34 known but non-conserved miRNAs were also found. Significantly, 72 new potential Amur grape-specific miRNAs were discovered. The sequences of these new potential va-miRNAs were further validated through miR-RACE, and accumulation of 18 new va-miRNAs in seven tissues of grapevines confirmed by real time RT-PCR (qRT-PCR) analysis. The expression levels of va-miRNAs in flowers and berries were found to be basically consistent in identity to those from deep sequenced sRNAs libraries of combined corresponding tissues. We also describe the conservation and variation of va-miRNAs using miR-SNPs and miR-LDs during plant evolution based on comparison of orthologous sequences, and further reveal that the number and sites of miR-SNP in diverse miRNA families exhibit distinct divergence. Finally, 346 target genes for the new miRNAs were predicted and they include a number of Amur grape stress tolerance genes and many genes regulating anthocyanin synthesis and sugar metabolism.

CONCLUSIONS

Deep sequencing of short RNAs from Amur grape flowers and berries identified 72 new potential miRNAs and 34 known but non-conserved miRNAs, indicating that specific miRNAs exist in Amur grape. These results show that a number of regulatory miRNAs exist in Amur grape and play an important role in Amur grape growth, development, and response to abiotic or biotic stress.

Identification of microRNAs from Amur grape (Vitis amurensis Rupr.) by deep sequencing and analysis of microRNA variations with bioinformatics

Background

MicroRNA (miRNA) is a class of functional non-coding small RNA with 19-25 nucleotides in length while Amur grape (Vitis amurensis Rupr.) is an important wild fruit crop with the strongest cold resistance among the Vitis species, is used as an excellent breeding parent for grapevine, and has elicited growing interest in wine production. To date, there is a relatively large number of grapevine miRNAs (vv-miRNAs) from cultivated grapevine varieties such as Vitis vinifera L. and hybrids of V. vinifera and V. labrusca, but there is no report on miRNAs from Vitis amurensis Rupr, a wild grapevine species.

Results

A small RNA library from Amur grape was constructed and Solexa technology used to perform deep sequencing of the library followed by subsequent bioinformatics analysis to identify new miRNAs. In total, 126 conserved miRNAs belonging to 27 miRNA families were identified, and 34 known but non-conserved miRNAs were also found. Significantly, 72 new potential Amur grape-specific miRNAs were discovered. The sequences of these new potential va-miRNAs were further validated through miR-RACE, and accumulation of 18 new va-miRNAs in seven tissues of grapevines confirmed by real time RT-PCR (qRT-PCR) analysis. The expression levels of va-miRNAs in flowers and berries were found to be basically consistent in identity to those from deep sequenced sRNAs libraries of combined corresponding tissues. We also describe the conservation and variation of va-miRNAs using miR-SNPs and miR-LDs during plant evolution based on comparison of orthologous sequences, and further reveal that the number and sites of miR-SNP in diverse miRNA families exhibit distinct divergence. Finally, 346 target genes for the new miRNAs were predicted and they include a number of Amur grape stress tolerance genes and many genes regulating anthocyanin synthesis and sugar metabolism.

Conclusions

Deep sequencing of short RNAs from Amur grape flowers and berries identified 72 new potential miRNAs and 34 known but non-conserved miRNAs, indicating that specific miRNAs exist in Amur grape. These results show that a number of regulatory miRNAs exist in Amur grape and play an important role in Amur grape growth, development, and response to abiotic or biotic stress.

Computational identification of MicroRNAs in strawberry expressed sequence tags and validation of their precise sequences by miR-RACE

MicroRNAs (miRNAs) are small, endogenously expressed, nonprotein-coding RNAs that regulate gene expression at the post-transcriptional level in both animals and plants through repressing translation or inducing mRNA degradation. A comprehensive strategy to identify new miRNA homologs by mining the repository of available strawberry expressed sequence tags (ESTs) was developed. By adopting a range of filtering criteria, we identified 11 potential miRNAs belonging to 5 miRNA families from 47 890 Fragaria vesca EST sequences. Using 2 specific 5' and 3' miRNA RACE PCR reactions and a sequence-directed cloning method, we accurately determined both end sequences of 5 candidate miRNAs. Meanwhile, qRT-PCR was used to detect the expression of these 5 miRNAs in different strawberry organs and tissues at several growing stages. These newly identified F. vesca miRNAs (fve-miRNAs) and their expression information can improve our understanding of possible roles of fve-miRNAs in regulating the growth and development of F. vesca.

Identification of conserved microRNAs and their targets in Chinese cabbage (Brassica rapa subsp. pekinensis)

The microRNAs (miRNAs) are a new class of small nonprotein-coding RNAs that have been identified to regulate gene expression at the post-transcriptional level by targeting mRNAs for degradation or by inhibiting protein translation. Until now, thousands of miRNAs have been identified in many plants species. However, only 23 miRNAs have been reported from the microRNA database in Chinese cabbage ( Brassica rapa subsp. pekinensis ), one of the most widely cultivated vegetables in China and East Asia. In the present study, 168 potential miRNAs, derived from 22 EST and 119 GSS sequences in Chinese cabbage were identified and classified into 38 miRNA families by well-defined computational analysis, in which most belonged to the miRNA1533, miRNA156, and miRNA2911 families. Totally, there are 129 identified miRNAs potentially targeting 1386 Chinese cabbage EST genes, which play roles in multiple biological and metabolic processes including metabolism, cell growth, signal transduction, stress response, and plant development. Gene ontology analysis, based on these target proteins, showed that 688, 532, and 287 genes were involved in molecular functions, biological processes, and cellular components, respectively. KEGG pathway analysis demonstrated that these miRNAs participated in 214 metabolism pathways, including, amongst others, plant–pathogen interaction, fatty acid metabolism, amino acid metabolism, nitrogen metabolism, plant hormone signal transduction.

Characterization of microRNAs identified in a table grapevine cultivar with validation of computationally predicted grapevine miRNAs by miR-RACE

Background

Alignment analysis of the Vv-miRNAs identified from various grapevine cultivars indicates that over 30% orthologous Vv-miRNAs exhibit a 1–3 nucleotide discrepancy only at their ends, suggesting that this sequence discrepancy is not a random event, but might mainly derive from divergence of cultivars. With advantages of miR-RACE technology in determining precise sequences of potential miRNAs from bioinformatics prediction, the precise sequences of vv-miRNAs predicted computationally can be verified with miR-RACE in a different grapevine cultivar. This presents itself as a new approach for large scale discovery of precise miRNAs in different grapevine varieties.

Methodology/Principal Findings

Among 88 unique sequences of Vv-miRNAs from bioinformatics prediction, 83 (96.3%) were successfully validated with MiR-RACE in grapevine cv. ‘Summer Black’. All the validated sequences were identical to their corresponding ones obtained from deep sequencing of the small RNA library of ‘Summer Black’. Quantitative RT-PCR analysis of the expressions levels of 10 Vv-miRNA/target gene pairs in grapevine tissues showed some negative correlation trends. Finally, comparison of Vv-miRNA sequences with their orthologs in Arabidopsis and study on the influence of divergent bases of the orthologous miRNAs on their targeting patterns in grapevine were also done.

Conclusion

The validation of precise sequences of potential Vv-miRNAs from computational prediction in a different grapevine cultivar can be a new way to identify the orthologous Vv-miRNAs. Nucleotide discrepancy of orthologous Vv-miRNAs from different grapevine cultivars normally does not change their target genes. However, sequence variations of some orthologous miRNAs in grapevine and Arabidopsis can change their targeting patterns. These precise Vv-miRNAs sequences validated in our study could benefit some further study on grapevine functional genomics.

Computational identification of microRNAs in peach expressed sequence tags and validation of their precise sequences by miR-RACE

Twenty-two potential miRNAs from seven miRNA families were first predicted from more than 80,857 EST sequences of peach (Prunus persica). Using two specific 5' and 3' miRNA RACE (miR-RACE) PCR reactions and sequence-directed cloning, we accurately determined the precise sequences, especially both ends, of eight candidate miRNAs. The sequencing results demonstrated that the ppe-miRNAs were conserved to those that were predicted computationally except ppe-miR171b. We validated the existence of two members (ppe-miR171a and miR171b) of the miR171 family in peach that belonged to different precursors. qRT-PCR was further employed in analyzing expression of the eight miRNAs in peach leaves, flowers, and fruits at different developing stages, where some of the miRNAs showed tissue-specific expression.

Computational identification of microRNAs in peach expressed sequence tags and validation of their precise sequences by miR-RACE

Twenty-two potential miRNAs from seven miRNA families were first predicted from more than 80,857 EST sequences of peach (Prunus persica). Using two specific 5' and 3' miRNA RACE (miR-RACE) PCR reactions and sequence-directed cloning, we accurately determined the precise sequences, especially both ends, of eight candidate miRNAs. The sequencing results demonstrated that the ppe-miRNAs were conserved to those that were predicted computationally except ppe-miR171b. We validated the existence of two members (ppe-miR171a and miR171b) of the miR171 family in peach that belonged to different precursors. qRT-PCR was further employed in analyzing expression of the eight miRNAs in peach leaves, flowers, and fruits at different developing stages, where some of the miRNAs showed tissue-specific expression.