Deep sequencing of tomato short RNAs identifies microRNAs targeting genes involved in fruit ripening

Interaction of small RNA-8105 and the intron of IbMYB1 RNA regulates IbMYB1 family genes through secondary siRNAs and DNA methylation after wounding

Small RNAs (sRNAs) play important roles in plants under stress conditions. However, limited research has been performed on the sRNAs involved in plant wound responses. In the present study, a novel wounding-induced sRNA, sRNA8105, was identified in sweet potato (Ipomoea batatas cv. Tainung 57) using microarray analysis. It was found that expression of sRNA8105 increased after mechanical wounding. Furthermore, Dicer-like 1 (DCL1) is required for the sRNA8105 precursor (pre-sRNA8105) to generate 22 and 24 nt mature sRNA8105. sRNA8105 targeted the first intron of IbMYB1 (MYB domain protein 1) before RNA splicing, and mediated RNA cleavage and DNA methylation of IbMYB1. The interaction between sRNA8105 and IbMYB1 was confirmed by cleavage site mapping, agro-infiltration analyses, and use of a transgenic sweet potato over-expressing pre-sRNA8105 gene. Induction of IbMYB1-siRNA was observed in the wild-type upon wounding and in transgenic sweet potato over-expressing pre-sRNA8105 gene without wounding, resulting in decreased expression of the whole IbMYB1 gene family, i.e. IbMYB1 and the IbMYB2 genes, and thus directing metabolic flux toward biosynthesis of lignin in the phenylpropanoid pathway. In conclusion, sRNA8105 induced by wounding binds to the first intron of IbMYB1 RNA to methylate IbMYB1, cleave IbMYB1 RNA, and trigger production of secondary siRNAs, further repressing the expression of the IbMYB1 family genes and regulating the phenylpropanoid pathway.

Analysis of the global transcriptome of longan (Dimocarpus longan Lour.) embryogenic callus using Illumina paired-end sequencing

Background

Longan is a tropical/subtropical fruit tree of great economic importance in Southeast Asia. Progress in understanding molecular mechanisms of longan embryogenesis, which is the primary influence on fruit quality and yield, is slowed by lack of transcriptomic and genomic information. Illumina second generation sequencing, which is suitable for generating enormous numbers of transcript sequences that can be used for functional genomic analysis of longan.

Results

In this study, a longan embryogenic callus (EC) cDNA library was sequenced using an Illumina HiSeq 2000 system. A total of 64,876,258 clean reads comprising 5.84 Gb of nucleotides were assembled into 68,925 unigenes of 448-bp mean length, with unigenes ≥1000 bp accounting for 8.26% of the total. Using BLASTx, 40,634 unigenes were found to have significant similarity with accessions in Nr and Swiss- Prot databases. Of these, 38,845 unigenes were assigned to 43 GO sub-categories and 17,118 unigenes were classified into 25 COG sub-groups. In addition, 17,306 unigenes mapped to 199 KEGG pathways, with the categories of Metabolic pathways, Plant-pathogen interaction, Biosynthesis of secondary metabolites, and Genetic information processing being well represented. Analyses of unigenes ≥1000 bp revealed 328 embryogenesis-related unigenes as well as numerous unigenes expressed in EC associated with functions of reproductive growth, such as flowering, gametophytogenesis, and fertility, and vegetative growth, such as root and shoot growth. Furthermore, 23 unigenes related to embryogenesis and reproductive and vegetative growth were validated by quantitative real time PCR (qPCR) in samples from different stages of longan somatic embryogenesis (SE); their differentially expressions in the various embryogenic cultures indicated their possible roles in longan SE.

Conclusions

The quantity and variety of expressed EC genes identified in this study is sufficient to serve as a global transcriptome dataset for longan EC and to provide more molecular resources for longan functional genomics.

Molecular and genetic regulation of fruit ripening

Fleshy fruit undergo a novel developmental program that ends in the irreversible process of ripening and eventual tissue senescence. During this maturation process, fruit undergo numerous physiological, biochemical and structural alterations, making them more attractive to seed dispersal organisms. In addition, advanced or over-ripening and senescence, especially through tissue softening and eventual decay, render fruit susceptible to invasion by opportunistic pathogens. While ripening and senescence are often used interchangeably, the specific metabolic activities of each would suggest that ripening is a distinct process of fleshy fruits that precedes and may predispose the fruit to subsequent senescence.

Molecular and genetic regulation of fruit ripening

Fleshy fruit undergo a novel developmental program that ends in the irreversible process of ripening and eventual tissue senescence. During this maturation process, fruit undergo numerous physiological, biochemical and structural alterations, making them more attractive to seed dispersal organisms. In addition, advanced or over-ripening and senescence, especially through tissue softening and eventual decay, render fruit susceptible to invasion by opportunistic pathogens. While ripening and senescence are often used interchangeably, the specific metabolic activities of each would suggest that ripening is a distinct process of fleshy fruits that precedes and may predispose the fruit to subsequent senescence.

Characterization of target mRNAs for grapevine microRNAs with an integrated strategy of modified RLM-RACE, newly developed PPM-RACE and qPCRs

MicroRNAs (miRNAs) regulate target gene expression by mediating target gene cleavage or inhibition of translation at transcriptional and post-transcriptional levels in higher plants. Until now, many grapevine microRNAs (Vv-miRNAs) have been identified and quite a number of miRNA target genes were also verified by various analysis. However, global interaction of miRNAs with their target genes still remained to perform more research. We reported experimental validation of a number of miRNA target genes in table grapevine that had been previously identified by bioinformatics in our earlier studies. To verify more predicted target genes of Vv-miRNAs and elucidate the modes by which these Vv-miRNAs work on their target genes, 31 unverified potential target genes for 18 Vv-miRNAs were experimentally verified by a new integrated strategy employing a modified 5'-RLM-RACE (RNA ligase-mediated 5' rapid amplification of cDNA ends), 3'-PPM-RACE (poly(A) polymerase-mediated 3' rapid amplification of cDNA ends) and qRT-PCRs of cleavage products. The results showed that these Vv-miRNAs negatively regulated expression of their target messenger RNAs (mRNAs) through guiding corresponding target mRNA cleavage, of which about 94.4% Vv-miRNAs cleaved their target mRNAs mainly at the tenth nucleotide of 5'-end of miRNAs. Expression levels of both miRNAs and their target mRNAs in eight tissues exhibited inverse relationships, and expressions both of cleaved targets and miRNAs indicated a cleavage mode of Vv-miRNAs on their target genes. Our results confirm the importance of Vv-miRNAs in grapevine growth and development, and suggest more study on Vv-miRNAs and targets can enrich the knowledge of miRNA mediated-regulation in grapevine.

Insilico profiling of microRNAs in Korean ginseng (Panax ginseng Meyer)

MicroRNAs (miRNAs) are a class of recently discovered non-coding small RNA molecules, on average approximately 21 nucleotides in length, which underlie numerous important biological roles in gene regulation in various organisms. The miRNA database (release 18) has 18,226 miRNAs, which have been deposited from different species. Although miRNAs have been identified and validated in many plant species, no studies have been reported on discovering miRNAs in Panax ginseng Meyer, which is a traditionally known medicinal plant in oriental medicine, also known as Korean ginseng. It has triterpene ginseng saponins called ginsenosides, which are responsible for its various pharmacological activities. Predicting conserved miRNAs by homology-based analysis with available expressed sequence tag (EST) sequences can be powerful, if the species lacks whole genome sequence information. In this study by using the EST based computational approach, 69 conserved miRNAs belonging to 44 miRNA families were identified in Korean ginseng. The digital gene expression patterns of predicted conserved miRNAs were analyzed by deep sequencing using small RNA sequences of flower buds, leaves, and lateral roots. We have found that many of the identified miRNAs showed tissue specific expressions. Using the insilico method, 346 potential targets were identified for the predicted 69 conserved miRNAs by searching the ginseng EST database, and the predicted targets were mainly involved in secondary metabolic processes, responses to biotic and abiotic stress, and transcription regulator activities, as well as a variety of other metabolic processes.

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.

Ineffective delivery of diet-derived microRNAs to recipient animal organisms

Cross-kingdom delivery of specific microRNAs to recipient organisms via food ingestion has been reported recently. However, it is unclear if such delivery of microRNAs occurs frequently in animal organisms after typical dietary intake. We found substantial levels of specific microRNAs in diets commonly consumed orally by humans, mice, and honey bees. Yet, after ingestion of fruit replete with plant microRNAs (MIR156a, MIR159a, and MIR169a), a cohort of healthy athletes did not carry detectable plasma levels of those molecules. Similarly, despite consumption of a diet with animal fat replete in endogenous miR-21, negligible expression of miR-21 in plasma or organ tissue was observed in miR-21 -/- recipient mice. Correspondingly, when fed vegetarian diets containing the above plant microRNAs, wild-type recipient mice expressed insignificant levels of these microRNAs. Finally, despite oral uptake of pollen containing these plant microRNAs, negligible delivery of these molecules was observed in recipient honeybees. Therefore, we conclude that horizontal delivery of microRNAs via typical dietary ingestion is neither a robust nor a frequent mechanism to maintain steady-state microRNA levels in a variety of model animal organisms, thus defining the biological limits of these molecules in vivo.

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.

Identification of miRNAs involved in long-term simulated microgravity response in Solanum lycopersicum

To identify the miRNAs associated with the simulated microgravity response in plants and to ascertain the regulation network mediated by miRNAs under simulated microgravity conditions, we constructed a miRNA library by direct cloning method and analyzed the library. Seven miRNAs that are conserved in other plants were cloned for the first time in Solanum lycopersicum under simulated microgravity condition. The expressions of six of the seven miRNAs were up-regulated, especially by long-term simulated microgravity. Gene ontology analysis showed that most of the predicted targeted genes were involved in transcription regulation, signal transduction and stress response, implying a complicated relationship among the external signal, internal transduction and final phenotype. Six of the predicted targets were validated by 5' RACE and reverse transcription real-time quantitative PCR. The results showed that with increasing miRNA expression levels, the corresponding target genes were down-regulated. The target gene of one of miRNAs, miR159e*, was thought to be associated with an increasing of starch amount under microgravity condition. A multi-stresses response network mediated by miRNAs under simulated microgravity condition was proposed. Cis-elements located in the upstream sequences of each miRNA were identified and their roles in gene regulation were investigated. In addition to the seven miRNAs that had homologs in other plants, six conserved S. lycopersicum miRNAs were identified. In the study, miRNAs were identified in S. lycopersicum for the first time under long-term simulated microgravity condition, which will help reveal the regulation mechanism mediated by miRNAs under simulated microgravity condition and adaptation to Earth's gravity.

Comparative analysis reveals dynamic changes in miRNAs and their targets and expression during somatic embryogenesis in longan (Dimocarpus longan Lour.)

Somatic embryogenesis (SE), which resembles zygotic embryogenesis, is an essential component of the process of plant cell differentiation and embryo development. Although microRNAs (miRNAs) are important regulators of many plant develop- mental processes, their roles in SE have not been thoroughly investigated. In this study, we used deep-sequencing, computational, and qPCR methods to identify, profile, and describe conserved and novel miRNAs involved in longan (Dimocarpus longan) SE. A total of 643 conserved and 29 novel miRNAs (including star strands) from more than 169 miRNA families were identified in longan embryogenic tissue using Solexa sequencing. By combining computational and degradome sequencing approaches, we were able to predict 2063 targets of 272 miRNAs and verify 862 targets of 181 miRNAs. Target annotation revealed that the putative targets were involved in a broad variety of biological processes, including plant metabolism, signal transduction, and stimulus response. Analysis of stage- and tissue-specific expressions of 20 conserved and 4 novel miRNAs indicated their possible roles in longan SE. These miRNAs were dlo-miR156 family members and dlo-miR166c* associated with early embryonic culture developmental stages; dlo-miR26, dlo-miR160a, and families dlo-miR159, dlo-miR390, and dlo-miR398b related to heart-shaped and torpedo- shaped embryo formation; dlo-miR4a, dlo-miR24, dlo-miR167a, dlo-miR168a*, dlo-miR397a, dlo-miR398b.1, dlo-miR398b.2, dlo-miR808 and dlo-miR5077 involved in cotyledonary embryonic development; and dlo-miR17 and dlo-miR2089*-1 that have regulatory roles during longan SE. In addition, dlo-miR167a, dlo-miR808, and dlo-miR5077 may be required for mature embryo formation. This study is the first reported investigation of longan SE involving large-scale cloning, characterization, and expression profiling of miRNAs and their targets. The reported results contribute to our knowledge of somatic embryo miRNAs and provide insights into miRNA biogenesis and expression in plant somatic embryo development.

Identification of microRNA targets in tomato fruit development using high-throughput sequencing and degradome analysis

MicroRNAs (miRNAs) play important roles in plant development through regulation of gene expression by mRNA degradation or translational inhibition. Despite the fact that tomato (Solanum lycopersicum) is the model system for studying fleshy fruit development and ripening, only a few experimentally proven miRNA targets are known, and the role of miRNA action in these processes remains largely unknown. Here, by using parallel analysis of RNA ends (PARE) for global identification of miRNA targets and comparing four different stages of tomato fruit development, a total of 119 target genes of miRNAs were identified. Of these, 106 appeared to be new targets. A large part of the identified targets (56) coded for transcription factors. Auxin response factors, as well as two known ripening regulators, colorless non-ripening (CNR) and APETALA2a (SlAP2a), with developmentally regulated degradation patterns were identified. The levels of the intact messenger of both CNR and AP2a are actively modulated during ripening, by miR156/157 and miR172, respectively. Additionally, two TAS3-mRNA loci were identified as targets of miR390. Other targets such as Argonaute 1 (AGO1), shown to be involved in miRNA biogenesis in other plant species, were identified, which suggests a feedback loop regulation of this process. In this study, it is shown that miRNA-guided cleavage of mRNAs is likely to play an important role in tomato fruit development and ripening.

Catalyzing plant science research with RNA-seq

Next generation DNA sequencing technologies are driving increasingly rapid, affordable and high resolution analyses of plant transcriptomes through sequencing of their associated cDNA (complementary DNA) populations; an analytical platform commonly referred to as RNA-sequencing (RNA-seq). Since entering the arena of whole genome profiling technologies only a few years ago, RNA-seq has proven itself to be a powerful tool with a remarkably diverse range of applications, from detailed studies of biological processes at the cell type-specific level, to providing insights into fundamental questions in plant biology on an evolutionary time scale. Applications include generating genomic data for heretofore unsequenced species, thus expanding the boundaries of what had been considered "model organisms," elucidating structural and regulatory gene networks, revealing how plants respond to developmental cues and their environment, allowing a better understanding of the relationships between genes and their products, and uniting the "omics" fields of transcriptomics, proteomics, and metabolomics into a now common systems biology paradigm. We provide an overview of the breadth of such studies and summarize the range of RNA-seq protocols that have been developed to address questions spanning cell type-specific-based transcriptomics, transcript secondary structure and gene mapping.

Regulation of ripening and opportunities for control in tomato and other fruits

Fruits are an important part of a healthy diet. They provide essential vitamins and minerals, and their consumption is associated with a reduced risk of heart disease and certain cancers. These important plant products can, however, be expensive to purchase, may be of disappointing quality and often have a short shelf life. A major challenge for crop improvement in fleshy fruit species is the enhancement of their health-promoting attributes while improving quality and reducing postharvest waste. To achieve these aims, a sound mechanistic understanding of the processes involved in fruit development and ripening is needed. In recent years, substantial insights have been made into the mechanistic basis of ethylene biosynthesis, perception and signalling and the identity of master regulators of ripening that operate upstream of, or in concert with a regulatory pathway mediated by this plant hormone. The role of other plant hormones in the ripening process has, however, remained elusive, and the links between regulators and downstream processes are still poorly understood. In this review, we focus on tomato as a model for fleshy fruit and provide an overview of the molecular circuits known to be involved in ripening, especially those controlling pigment accumulation and texture changes. We then discuss how this information can be used to understand ripening in other fleshy fruit-bearing species. Recent developments in comparative genomics and systems biology approaches are discussed. The potential role of epigenetic changes in generating useful variation is highlighted along with opportunities for enhancing the level of metabolites that have a beneficial effect on human health.

Functional genomics of tomato in a post-genome-sequencing phase

Completion of tomato genome sequencing project has broad impacts on genetic and genomic studies of tomato and Solanaceae plants. The reference genome sequence derived from Solanum lycopersicum cv 'Heinz 1706' serves as the firm basis for sequencing-based approaches to tomato genomics. In this article, we first present a brief summary of the genome sequencing project and a summary of the reference genome sequence. We then focus on recent progress in transcriptome sequencing and small RNA sequencing and show how the reference genome sequence makes these analyses more comprehensive than before. We discuss the potential of in-depth analysis that is based on DNA methylome sequencing and transcription start-site detection. Finally, we describe the current status of efforts to resequence S. lycopersicum cultivars to demonstrate how resequencing can allow the use of intraspecific genomic diversity for detailed phenotyping and breeding.

Identification and characterization of novel and conserved microRNAs in radish (Raphanus sativus L.) using high-throughput sequencing

MicroRNAs (miRNAs) are endogenous, non-coding, small RNAs that play significant regulatory roles in plant growth, development, and biotic and abiotic stress responses. To date, a great number of conserved and species-specific miRNAs have been identified in many important plant species such as Arabidopsis, rice and poplar. However, little is known about identification of miRNAs and their target genes in radish (Raphanus sativus L.). In the present study, a small RNA library from radish root was constructed and sequenced using the high-throughput Solexa sequencing. Through sequence alignment and secondary structure prediction, a total of 545 conserved miRNA families as well as 15 novel (with their miRNA* strand) and 64 potentially novel miRNAs were identified. Quantitative real-time PCR (qRT-PCR) analysis confirmed that both conserved and novel miRNAs were expressed in radish, and some of them were preferentially expressed in certain tissues. A total of 196 potential target genes were predicted for 42 novel radish miRNAs. Gene ontology (GO) analysis showed that most of the targets were involved in plant growth, development, metabolism and stress responses. This study represents a first large-scale identification and characterization of radish miRNAs and their potential target genes. These results could lead to the further identification of radish miRNAs and enhance our understanding of radish miRNA regulatory mechanisms in diverse biological and metabolic processes.

A deep sequencing approach to uncover the miRNOME in the human heart

MicroRNAs (miRNAs) are a class of non-coding RNAs of ∼22 nucleotides in length, and constitute a novel class of gene regulators by imperfect base-pairing to the 3′UTR of protein encoding messenger RNAs. Growing evidence indicates that miRNAs are implicated in several pathological processes in myocardial disease. The past years, we have witnessed several profiling attempts using high-density oligonucleotide array-based approaches to identify the complete miRNA content (miRNOME) in the healthy and diseased mammalian heart. These efforts have demonstrated that the failing heart displays differential expression of several dozens of miRNAs. While the total number of experimentally validated human miRNAs is roughly two thousand, the number of expressed miRNAs in the human myocardium remains elusive. Our objective was to perform an unbiased assay to identify the miRNOME of the human heart, both under physiological and pathophysiological conditions. We used deep sequencing and bioinformatics to annotate and quantify microRNA expression in healthy and diseased human heart (heart failure secondary to hypertrophic or dilated cardiomyopathy). Our results indicate that the human heart expresses >800 miRNAs, the majority of which not being annotated nor described so far and some of which being unique to primate species. Furthermore, >250 miRNAs show differential and etiology-dependent expression in human dilated cardiomyopathy (DCM) or hypertrophic cardiomyopathy (HCM). The human cardiac miRNOME still possesses a large number of miRNAs that remain virtually unexplored. The current study provides a starting point for a more comprehensive understanding of the role of miRNAs in regulating human heart disease.

Identification and characterization of miRNA transcriptome in potato by high-throughput sequencing

Micro RNAs (miRNAs) represent a class of short, non-coding, endogenous RNAs which play important roles in post-transcriptional regulation of gene expression. While the diverse functions of miRNAs in model plants have been well studied, the impact of miRNAs in crop plant biology is poorly understood. Here we used high-throughput sequencing and bioinformatics analysis to analyze miRNAs in the tuber bearing crop potato (Solanum tuberosum). Small RNAs were analysed from leaf and stolon tissues. 28 conserved miRNA families were found and potato-specific miRNAs were identified and validated by RNA gel blot hybridization. The size, origin and predicted targets of conserved and potato specific miRNAs are described. The large number of miRNAs and complex population of small RNAs in potato suggest important roles for these non-coding RNAs in diverse physiological and metabolic pathways.

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.

Deep sequencing of maize small RNAs reveals a diverse set of microRNA in dry and imbibed seeds

Seed germination plays a pivotal role during the life cycle of plants. As dry seeds imbibe water, the resumption of energy metabolism and cellular repair occur and miRNA-mediated gene expression regulation is involved in the reactivation events. This research was aimed at understanding the role of miRNA in the molecular control during seed imbibition process. Small RNA libraries constructed from dry and imbibed maize seed embryos were sequenced using the Illumina platform. Twenty-four conserved miRNA families were identified in both libraries. Sixteen of them showed significant expression differences between dry and imbibed seeds. Twelve miRNA families, miR156, miR159, miR164, miR166, miR167, miR168, miR169, miR172, miR319, miR393, miR394 and miR397, were significantly down-regulated; while four families, miR398, miR408, miR528 and miR529, were significantly up-regulated in imbibed seeds compared to that in dry seeds. Furthermore, putative novel maize miRNAs and their target genes were predicted. Target gene GO analysis was performed for novel miRNAs that were sequenced more than 50 times in the normalized libraries. The result showed that carbohydrate catabolic related genes were specifically enriched in the dry seed, while in imbibed seed target gene enrichment covered a broad range of functional categories including genes in amino acid biosynthesis, isomerase activity, ligase activity and others. The sequencing results were partially validated by quantitative RT-PCR for both conserved and novel miRNAs and the predicted target genes. Our data suggested that diverse and complex miRNAs are involved in the seed imbibition process. That miRNA are involved in plant hormone regulation may play important roles during the dry-imbibed seed transition.

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.

Diverse correlation patterns between microRNAs and their targets during tomato fruit development indicates different modes of microRNA actions

MicroRNAs negatively regulate the accumulation of mRNAs therefore when they are expressed in the same cells their expression profiles show an inverse correlation. We previously described one positively correlated miRNA/target pair, but it is not known how widespread this phenomenon is. Here, we investigated the correlation between the expression profiles of differentially expressed miRNAs and their targets during tomato fruit development using deep sequencing, Northern blot and RT-qPCR. We found an equal number of positively and negatively correlated miRNA/target pairs indicating that positive correlation is more frequent than previously thought. We also found that the correlation between microRNA and target expression profiles can vary between mRNAs belonging to the same gene family and even for the same target mRNA at different developmental stages. Since microRNAs always negatively regulate their targets, the high number of positively correlated microRNA/target pairs suggests that mutual exclusion could be as widespread as temporal regulation. The change of correlation during development suggests that the type of regulatory circuit directed by a microRNA can change over time and can be different for individual gene family members. Our results also highlight potential problems for expression profiling-based microRNA target identification/validation.

Identifying conserved and novel microRNAs in developing seeds of Brassica napus using deep sequencing

MicroRNAs (miRNAs) are important post-transcriptional regulators of plant development and seed formation. In Brassica napus, an important edible oil crop, valuable lipids are synthesized and stored in specific seed tissues during embryogenesis. The miRNA transcriptome of B. napus is currently poorly characterized, especially at different seed developmental stages. This work aims to describe the miRNAome of developing seeds of B. napus by identifying plant-conserved and novel miRNAs and comparing miRNA abundance in mature versus developing seeds. Members of 59 miRNA families were detected through a computational analysis of a large number of reads obtained from deep sequencing two small RNA and two RNA-seq libraries of (i) pooled immature developing stages and (ii) mature B. napus seeds. Among these miRNA families, 17 families are currently known to exist in B. napus; additionally 29 families not reported in B. napus but conserved in other plant species were identified by alignment with known plant mature miRNAs. Assembled mRNA-seq contigs allowed for a search of putative new precursors and led to the identification of 13 novel miRNA families. Analysis of miRNA population between libraries reveals that several miRNAs and isomiRNAs have different abundance in developing stages compared to mature seeds. The predicted miRNA target genes encode a broad range of proteins related to seed development and energy storage. This work presents a comparative study of the miRNA transcriptome of mature and developing B. napus seeds and provides a basis for future research on individual miRNAs and their functions in embryogenesis, seed maturation and lipid accumulation in B. napus.

Small RNA and transcriptome deep sequencing proffers insight into floral gene regulation in Rosa cultivars

Background

Roses (Rosa sp.), which belong to the family Rosaceae, are the most economically important ornamental plants—making up 30% of the floriculture market. However, given high demand for roses, rose breeding programs are limited in molecular resources which can greatly enhance and speed breeding efforts. A better understanding of important genes that contribute to important floral development and desired phenotypes will lead to improved rose cultivars. For this study, we analyzed rose miRNAs and the rose flower transcriptome in order to generate a database to expound upon current knowledge regarding regulation of important floral characteristics. A rose genetic database will enable comprehensive analysis of gene expression and regulation via miRNA among different Rosa cultivars.

Results

We produced more than 0.5 million reads from expressed sequences, totalling more than 110 million bp. From these, we generated 35,657, 31,434, 34,725, and 39,722 flower unigenes from Rosa hybrid: ‘Vital’, ‘Maroussia’, and ‘Sympathy’ and Rosa rugosa Thunb. , respectively. The unigenes were assigned functional annotations, domains, metabolic pathways, Gene Ontology (GO) terms, Plant Ontology (PO) terms, and MIPS Functional Catalogue (FunCat) terms. Rose flower transcripts were compared with genes from whole genome sequences of Rosaceae members (apple, strawberry, and peach) and grape. We also produced approximately 40 million small RNA reads from flower tissue for Rosa, representing 267 unique miRNA tags. Among identified miRNAs, 25 of them were novel and 242 of them were conserved miRNAs. Statistical analyses of miRNA profiles revealed both shared and species-specific miRNAs, which presumably effect flower development and phenotypes.

Conclusions

In this study, we constructed a Rose miRNA and transcriptome database, and we analyzed the miRNAs and transcriptome generated from the flower tissues of four Rosa cultivars. The database provides a comprehensive genetic resource which can be used to better understand rose flower development and to identify candidate genes for important phenotypes.

Identification and characterization of salt-responsive microRNAs in Populus tomentosa by high-throughput sequencing

Salt is one of the main environmental factors limiting plant growth and a better understanding of mechanisms of salt stress would aid efforts to bolster plant salt tolerance. MicroRNAs are well known for their important regulatory roles in response to abiotic stress in plants. In this study, high-throughput sequencing was employed to identify miRNAs in Populus tomentosa plantlets treated or not with salt (200 mM for 10 h). We found 141 conserved miRNAs belonging to 31 families, 29 non-conserved but previously-known miRNAs belonging to 26 families, and 17 novel miRNAs. Under salt stress, 19 miRNAs belonging to seven conserved miRNA families were significantly downregulated, and two miRNAs belonging to two conserved miRNA families were upregulated. Of seven non-conserved miRNAs with significantly altered expression, five were downregulated and two were upregulated. Furthermore, eight miRNAs were validated by qRT-PCR and their dynamic differential expressions were analyzed. In addition, 269 target genes of identified miRNAs were predicted and categorized by function. These results provide new insights into salt-responsive miRNAs in Populus.

Genome-wide identification of reverse complementary microRNA genes in plants

MicroRNAs (miRNAs) are ∼21-nucleotide small RNAs (sRNAs) with essential regulatory roles in plants. They are generated from stem-loop-structured precursors through two sequential Dicer-like 1 (DCL1)-mediated cleavages. To date, hundreds of plant miRNAs have been uncovered. However, the question, whether the sequences reverse complementary (RC) to the miRNA precursors could form hairpin-like structures and produce sRNA duplexes similar to the miRNA/miRNA* pairs has not been solved yet. Here, we interrogated this possibility in 16 plant species based on sRNA high-throughput sequencing data and secondary structure prediction. A total of 59 RC sequences with great potential to form stem-loop structures and generate miRNA/miRNA*-like duplexes were identified in ten plants, which were named as RC-miRNA precursors. Unlike the canonical miRNAs, only a few cleavage targets of the RC-miRNAs were identified in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), and none in Soybean (Glycine max) based on degradome data. Surprisingly, the genomic regions surrounding some of the RC-miRNA target recognition sites were observed to be specifically methylated in both Arabidopsis and rice. Taken together, we reported a new class of miRNAs, called RC-miRNAs, which were generated from the antisense strands of the miRNA precursors. Based on the results, we speculated that the mature RC-miRNAs might have subtle regulatory activity through target cleavages, but might possess short interfering RNA-like activity by guiding sequence-specific DNA methylation.

MicroR828 regulates lignin and H2O2 accumulation in sweet potato on wounding

MicroRNAs (miRNAs) are small noncoding RNAs which post-transcriptionally regulate gene expression by directing mRNA cleavage or translational inhibition. miRNAs play multiple roles in the growth, development and stress responses in plants. However, little is known of the wounding-responsive miRNAs and their regulation. Here, we investigated the expression patterns of microR828 (miR828) on wounding in sweet potato (Ipomoea batatas cv Tainung 57). The expression of miR828 was only detected in leaves, and was induced by wounding rather than by ethylene, hydrogen peroxide (H2O2), methyl jasmonate or nitric oxide (NO). Moreover, cyclic guanosine monophosphate (cGMP) was necessary for miR828 accumulation in leaves on wounding. Two miR828 target candidates, named IbMYB and IbTLD, were obtained by cDNA cloning, and their mRNA cleavage caused by miR828 was confirmed by cleavage site mapping, agro-infiltration and transgenics studies. The reduction in IbMYB and IbTLD expression coincided with the induction of miR828, demonstrating that IbMYB and IbTLD might be miR828 targets. Furthermore, transgenic sweet potato overexpressing miR828 precursor affected lignin and H2O2 contents. These results showed that cGMP could regulate wounding-responsive miR828, which repressed the expression of IbMYB and IbTLD. Subsequently, lignin and H2O2 were accumulated to participate in defense mechanisms.

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

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

Deep sequencing and microarray hybridization identify conserved and species-specific microRNAs during somatic embryogenesis in hybrid yellow poplar

BACKGROUND

To date, several studies have indicated a major role for microRNAs (miRNAs) in regulating plant development, but miRNA-mediated regulation of the developing somatic embryo is poorly understood, especially during early stages of somatic embryogenesis in hardwood plants. In this study, Solexa sequencing and miRNA microfluidic chips were used to discover conserved and species-specific miRNAs during somatic embryogenesis of hybrid yellow poplar (Liriodendron tulipifera×L. chinense).

METHODOLOGY/PRINCIPAL FINDINGS

A total of 17,214,153 reads representing 7,421,623 distinct sequences were obtained from a short RNA library generated from small RNAs extracted from all stages of somatic embryos. Through a combination of deep sequencing and bioinformatic analyses, we discovered 83 sequences with perfect matches to known miRNAs from 33 conserved miRNA families and 273 species-specific candidate miRNAs. MicroRNA microarray results demonstrated that many conserved and species-specific miRNAs were expressed in hybrid yellow poplar embryos. In addition, the microarray also detected another 149 potential miRNAs, belonging to 29 conserved families, which were not discovered by deep sequencing analysis. The biological processes and molecular functions of the targets of these miRNAs were predicted by carrying out BLAST search against Arabidopsis thaliana GenBank sequences and then analyzing the results with Gene Ontology.

CONCLUSIONS

Solexa sequencing and microarray hybridization were used to discover 232 candidate conserved miRNAs from 61 miRNA families and 273 candidate species-specific miRNAs in hybrid yellow poplar. In these predicted miRNAs, 64 conserved miRNAs and 177 species-specific miRNAs were detected by both sequencing and microarray hybridization. Our results suggest that miRNAs have wide-ranging characteristics and important roles during all stages of somatic embryogenesis in this economically important species.

Identification of miRNAs and their targets from Brassica napus by high-throughput sequencing and degradome analysis

BACKGROUND

MicroRNAs (miRNAs) are endogenous regulators of a broad range of physiological processes and act by either degrading mRNA or blocking its translation. Oilseed rape (Brassica napus) is one of the most important crops in China, Europe and other Asian countries with publicly available expressed sequence tags (ESTs) and genomic survey sequence (GSS) databases, but little is known about its miRNAs and their targets. To date, only 46 miRNAs have been identified in B. napus.

RESULTS

Forty-one conserved and 62 brassica-specific candidate B. napus miRNAs, including 20 miRNA* sequences, were identified using Solexa sequencing technology. Furthermore, 33 non-redundant mRNA targets of conserved brassica miRNAs and 19 new non-redundant mRNA targets of novel brassica-specific miRNAs were identified by genome-scale sequencing of mRNA degradome.

CONCLUSIONS

This study describes large scale cloning and characterization of B. napus miRNAs and their potential targets, providing the foundation for further characterization of miRNA function in the regulation of diverse physiological processes in B. napus.

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.

Identification and characterization of small non-coding RNAs from Chinese fir by high throughput sequencing

BACKGROUND

Small non-coding RNAs (sRNAs) play key roles in plant development, growth and responses to biotic and abiotic stresses. At least four classes of sRNAs have been well characterized in plants, including repeat-associated siRNAs (rasiRNAs), microRNAs (miRNAs), trans-acting siRNAs (tasiRNAs) and natural antisense transcript-derived siRNAs. Chinese fir (Cunninghamia lanceolata) is one of the most important coniferous evergreen tree species in China. No sRNA from Chinese fir has been described to date.

RESULTS

To obtain sRNAs in Chinese fir, we sequenced a sRNA library generated from seeds, seedlings, leaves, stems and calli, using Illumina high throughput sequencing technology. A comprehensive set of sRNAs were acquired, including conserved and novel miRNAs, rasiRNAs and tasiRNAs. With BLASTN and MIREAP we identified a total of 115 conserved miRNAs comprising 40 miRNA families and one novel miRNA with precursor sequence. The expressions of 16 conserved and one novel miRNAs and one tasiRNA were detected by RT-PCR. Utilizing real time RT-PCR, we revealed that four conserved and one novel miRNAs displayed developmental stage-specific expression patterns in Chinese fir. In addition, 209 unigenes were predicted to be targets of 30 Chinese fir miRNA families, of which five target genes were experimentally verified by 5' RACE, including a squamosa promoter-binding protein gene, a pentatricopeptide (PPR) repeat-containing protein gene, a BolA-like family protein gene, AGO1 and a gene of unknown function. We also demonstrated that the DCL3-dependent rasiRNA biogenesis pathway, which had been considered absent in conifers, existed in Chinese fir. Furthermore, the miR390-TAS3-ARF regulatory pathway was elucidated.

CONCLUSIONS

We unveiled a complex population of sRNAs in Chinese fir through high throughput sequencing. This provides an insight into the composition and function of sRNAs in Chinese fir and sheds new light on land plant sRNA evolution.

Analyzing the microRNA Transcriptome in Plants Using Deep Sequencing Data

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

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.

Discovery and profiling of novel and conserved microRNAs during flower development in Carya cathayensis via deep sequencing

Hickory (Carya cathayensis Sarg.) is an economically important woody plant in China, but its long juvenile phase delays yield. MicroRNAs (miRNAs) are critical regulators of genes and important for normal plant development and physiology, including flower development. We used Solexa technology to sequence two small RNA libraries from two floral differentiation stages in hickory to identify miRNAs related to flower development. We identified 39 conserved miRNA sequences from 114 loci belonging to 23 families as well as two novel and ten potential novel miRNAs belonging to nine families. Moreover, 35 conserved miRNA*s and two novel miRNA*s were detected. Twenty miRNA sequences from 49 loci belonging to 11 families were differentially expressed; all were up-regulated at the later stage of flower development in hickory. Quantitative real-time PCR of 12 conserved miRNA sequences, five novel miRNA families, and two novel miRNA*s validated that all were expressed during hickory flower development, and the expression patterns were similar to those detected with Solexa sequencing. Finally, a total of 146 targets of the novel and conserved miRNAs were predicted. This study identified a diverse set of miRNAs that were closely related to hickory flower development and that could help in plant floral induction.

Characterization of the small RNA component of leaves and fruits from four different cucurbit species

Background

MicroRNAs (miRNAs) are a class of non-coding small RNAs involved in post-transcriptional regulation of gene expression critical for plant growth and development, stress responses and other diverse biological processes in plants. The Cucurbitaceae or cucurbit family represents some of economically important species, particularly those with edible and medicinal fruits. Genomic tools for the molecular analysis of members of this family are just emerging. Partial draft genome sequence became available recently for cucumber and watermelon facilitating investigation of the small RNA component of the transcriptomes in cucurbits.

Results

We generated four small RNA libraries from bottle gourd (Lagenaria siceraria), Cucurbita moschata, Cucurbita pepo, and, watermelon (Citrullus lanatus var. lanatus) in order to identify conserved and novel lineage specific miRNAs in these cucurbits. Deep sequencing of small RNA libraries from these species resulted in 1,597,263, 532,948, 601,388, and 493,384 unique sRNA reads from bottle gourd, moschata, pepo and watermelon, respectively. Sequence analysis of these four libraries resulted in identification of 21 miRNA families that are highly conserved and 8 miRNA families that are moderately conserved in diverse dicots. We also identified 4 putative novel miRNAs in these plant species. Furthermore, the tasiRNAs were identified and their biogenesis was determined in these cucurbits. Small RNA blot analysis or q-PCR analyses of leaf and fruit tissues of these cucurbits showed differential expression of several conserved miRNAs. Interestingly, the abundance of several miRNAs in leaves and fruits of closely related C. moschata and C. pepo was also distinctly different. Target genes for the most conserved miRNAs are also predicted.

Conclusion

High-throughput sequencing of small RNA libraries from four cucurbit species has provided a glimpse of small RNA component in their transcriptomes. The analysis also showed considerable variation within four cucurbit species with regards to expression of individual miRNAs.

Two strawberry miR159 family members display developmental-specific expression patterns in the fruit receptacle and cooperatively regulate Fa-GAMYB

• We have reported previously that the gibberellin (GA) content in strawberry receptacle is high, peaking at specific stages, pointing to a role of this hormone in fruit development. In Arabidopsis, miR159 levels are dependent on GA concentration. This prompted us to investigate the role of two members of the miR159 family and their putative strawberry target gene, GAMYB, in relation to changes in GA content during the course of fruit development. • The highest expression level of the two Fa-MIR159 genes was in the fruit's receptacle tissue, with dramatic changes observed throughout development. The lowest levels of total mature miR159 (a and b) were observed during the white stage of receptacle development, which was concurrent with the highest expression of Fa-GAMYB. A functional interaction between miR159 and Fa-GAMYB has been demonstrated in receptacle tissue. • The application of bioactive GA (i.e. GA(3) ) to strawberry plants caused the down-regulated expression of Fa-MIR159a, but the expression of Fa-MIR159b was not affected significantly. Clear discrepancies between Fa-MIR159b and mature Fa-miR159b levels were indicative of post-transcriptional regulation of Fa-MIR159b gene expression. • We propose that Fa-miR159a and Fa-miR159b interact with Fa-GAMYB during the course of strawberry receptacle development, and that they act in a cooperative fashion to respond, in part, to changes in GA endogenous levels.

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.

Uncovering small RNA-mediated responses to cold stress in a wheat thermosensitive genic male-sterile line by deep sequencing

The male sterility of thermosensitive genic male sterile (TGMS) lines of wheat (Triticum aestivum) is strictly controlled by temperature. The early phase of anther development is especially susceptible to cold stress. MicroRNAs (miRNAs) play an important role in plant development and in responses to environmental stress. In this study, deep sequencing of small RNA (smRNA) libraries obtained from spike tissues of the TGMS line under cold and control conditions identified a total of 78 unique miRNA sequences from 30 families and trans-acting small interfering RNAs (tasiRNAs) derived from two TAS3 genes. To identify smRNA targets in the wheat TGMS line, we applied the degradome sequencing method, which globally and directly identifies the remnants of smRNA-directed target cleavage. We identified 26 targets of 16 miRNA families and three targets of tasiRNAs. Comparing smRNA sequencing data sets and TaqMan quantitative polymerase chain reaction results, we identified six miRNAs and one tasiRNA (tasiRNA-ARF [for Auxin-Responsive Factor]) as cold stress-responsive smRNAs in spike tissues of the TGMS line. We also determined the expression profiles of target genes that encode transcription factors in response to cold stress. Interestingly, the expression of cold stress-responsive smRNAs integrated in the auxin-signaling pathway and their target genes was largely noncorrelated. We investigated the tissue-specific expression of smRNAs using a tissue microarray approach. Our data indicated that miR167 and tasiRNA-ARF play roles in regulating the auxin-signaling pathway and possibly in the developmental response to cold stress. These data provide evidence that smRNA regulatory pathways are linked with male sterility in the TGMS line during cold stress.

SoMART: a web server for plant miRNA, tasiRNA and target gene analysis

Plant microRNAs (miRNAs) and trans-acting small interfering RNAs (tasiRNAs) play important roles in a variety of biological processes. Bioinformatics prediction and small RNA (sRNA) cloning are the most important approaches for identification of miRNAs and tasiRNAs and their targets. However, these approaches are not readily accessible to every researcher. Here we present SoMART, a web server for miRNA/tasiRNA analysis resources and tools, which is designed for researchers who are interested in identifying miRNAs or tasiRNAs that potentially regulate genes of interest. The server includes four sets of tools: 'Slicer detector' for detecting sRNAs targeting input genes, 'dRNA mapper' for detecting degradome (d)RNA products derived from input genes, 'PreMIR detector' for identifying miRNA precursors (MIRs) or tasiRNA precursor (TASs) of input sRNAs, and 'sRNA mapper' for mapping sRNAs onto input genes. We also developed a dRNA-seq protocol to achieve longer dRNA reads for better characterization of miRNA precursors by dRNA mapper. To validate the server function and robustness, we installed sRNA, dRNA and collected genomic DNA or transcriptome databases from Arabidopsis and solanaceous plants, and characterized miR172-mediated regulation of the APETALA2 gene in potato (Solanum tuberosum) and demonstrated conservation of MIR390-triggered TAS3 in tomato (Solanum lycopersicum). More importantly, we predicted the existence of MIR482-triggered TAS5 in tomato. We further tested and confirmed the efficiency and accuracy of the server by analyses of 21 validated miRNA targets and 115 miRNA precursors in Arabidopsis thaliana. SoMART is available at http://somart.ist.berkeley.edu.

Genomic organization, phylogenetic comparison and differential expression of the SBP-box family of transcription factors in tomato

SBP-box genes represent transcription factors ubiquitously found in the plant kingdom and recognized as important regulators of many different aspects of plant development. In this study, 15 SBP-box gene family members were identified in tomato and analyzed with respect to their genomic organization and other structural features. Phylogenetic reconstruction based on the DNA-binding SBP-domain, allowed the classification of the SlySBP proteins into eight groups representing clear orthologous relationships to family members of other flowering plants and the moss Physcomitrella. In order to have a better understanding of their possible function in the development of a fleshy-fruit species like tomato, the mRNA expression levels of all SlySBP genes were quantified in vegetative and reproductive organs of plants, at different stages of growth. As transcripts of ten SlySBP genes were found to carry putative miR156- and miR157-response elements, the expression levels of the corresponding microRNAs were determined as well, revealing different patterns of expression. In addition, eight putative miR156 and four miR157 encoding loci could be identified in the tomato genome, four of them forming a polycistronic cluster. Whereas miR156 and miR157 levels were highest in seedlings, leaves and anthers of young flowers, most miR156-targeted SlySBP genes were found to be expressed in young inflorescences and during fruit development and ripening, suggesting a particularly important role during tomato reproductive growth. The data presented provide a basis for future clarification of the various functions that SBP-box gene family members play in tomato growth and development.

Genome-wide identification and expression analysis of heat-responsive and novel microRNAs in Populus tomentosa

Plant microRNAs have a vital role in various abiotic stress responses by regulating gene expression. Heat stress is one of the most severe abiotic stresses, and affects plant growth and development, even leading to death. To identify heat-responsive miRNAs at the genome-wide level in Populus, Solexa sequencing was employed to sequence two libraries from Populus tomentosa, treated and untreated by heat stress. Sequence analysis identified 134 conserved miRNAs belonging to 30 miRNA families, and 16 novel miRNAs belonging to 14 families. Among these miRNAs, 52 miRNAs from 15 families were responsive to heat stress and most of them were down-regulated. qRT-PCR analysis confirmed that the conserved and novel miRNAs were expressed in P. tomentosa, and revealed similar expression trends to the Solexa sequencing results obtained under heat stress. One hundred and nine targets of the novel miRNAs were predicted. This study opens up a new avenue for understanding the regulatory mechanisms of miRNAs involvement in the heat stress response of trees.

Genome-wide profiling of miRNAs and other small non-coding RNAs in the Verticillium dahliae-inoculated cotton roots

MicroRNAs (miRNAs) and small interfering RNAs (siRNAs) are short (19–25 nucleotides) non-coding RNA molecules that have large-scale regulatory effects on development and stress responses in plants. Verticillium wilt is a vascular disease in plants caused by the fungal pathogen Verticillium dahliae. The objective of this study is to investigate the transcriptional profile of miRNAs and other small non-coding RNAs in Verticillium–inoculated cotton roots. Four small RNA libraries were constructed from mocked and infected roots of two cotton cultured species which are with different Verticillium wilt tolerance (‘Hai-7124’, Gossypium barbadense L., a Verticillium-tolerant cultivar, and ‘Yi-11’, Gossypium hirsutum L. a Verticillium-sensitive cultivar). The length distribution of obtained small RNAs was significantly different between libraries. There were a total of 215 miRNA families identified in the two cotton species. Of them 14 were novel miRNAs. There were >65 families with different expression between libraries. We also identified two trans-acting siRNAs and thousands of endogenous siRNA candidates, and hundred of them exhibited altered expression after inoculation of Verticillium. Interesting, many siRNAs were found with a perfect match with retrotransposon sequences, suggested that retrotransposons maybe one of sources for the generation of plant endogenous siRNAs. The profiling of these miRNAs and other small non-coding RNAs lay the foundation for further understanding of small RNAs function in the regulation of Verticillium defence responses in cotton roots.

Transcriptome-wide identification and characterization of miRNAs from Pinus densata

Background

MicroRNAs (miRNAs) play key roles in diverse developmental processes, nutrient homeostasis and responses to biotic and abiotic stresses. The biogenesis and regulatory functions of miRNAs have been intensively studied in model angiosperms, such as Arabidopsis thaliana, Oryza sativa and Populus trichocarpa. However, global identification of Pinus densata miRNAs has not been reported in previous research.

Results

Here, we report the identification of 34 conserved miRNAs belonging to 25 miRNA families from a P. densata mRNA transcriptome database using local BLAST and MIREAP programs. The primary and/or precursor sequences of 29 miRNAs were further confirmed by RT-PCR amplification and subsequent sequencing. The average value of the minimal folding free energy indexes of the 34 miRNA precursors was 0.92. Nineteen (58%) mature miRNAs began with a 5' terminal uridine residue. Analysis of miRNA precursors showed that 19 mature miRNAs were novel members of 14 conserved miRNA families, of which 17 miRNAs were further validated by subcloning and sequencing. Using real-time quantitative RT-PCR, we found that the expression levels of 7 miRNAs were more than 2-fold higher in needles than in stems. In addition, 72 P. densata mRNAs were predicted to be targets of 25 miRNA families. Four target genes, including a nodal modulator 1-like protein gene, two GRAS family transcription factor protein genes and one histone deacetylase gene, were experimentally verified to be the targets of 3 P. densata miRNAs, pde-miR162a, pde-miR171a and pde-miR482a, respectively.

Conclusions

This study led to the discovery of 34 conserved miRNAs comprising 25 miRNA families from Pinus densata. These results lay a solid foundation for further studying the regulative roles of miRNAs in the development, growth and responses to environmental stresses in P. densata.

Global analysis of tomato gene expression during Potato spindle tuber viroid infection reveals a complex array of changes affecting hormone signaling

Viroids like Potato spindle tuber viroid (PSTVd) are the smallest known agents of infectious disease-small, highly structured, circular RNA molecules that lack detectable messenger RNA activity, yet are able to replicate autonomously in susceptible plant species. To better understand the possible role of RNA silencing in disease induction, a combination of microarray analysis and large-scale RNA sequence analysis was used to compare changes in tomato gene expression and microRNA levels associated with PSTVd infection in two tomato cultivars plus a third transformed line expressing small PSTVd small interfering RNAs in the absence of viroid replication. Changes in messenger (m)RNA levels for the sensitive cultivar 'Rutgers' were extensive, involving more than half of the approximately 10,000 genes present on the array. Chloroplast biogenesis was down-regulated in both sensitive and tolerant cultivars, and effects on mRNAs encoding enzymes involved in the biosynthesis of gibberellin and other hormones were accompanied by numerous changes affecting their respective signaling pathways. In the dwarf cultivar 'MicroTom', a marked upregulation of genes involved in response to stress and other stimuli was observed only when exogenous brassinosteroid was applied to infected plants, thereby providing the first evidence for the involvement of brassinosteroid-mediated signaling in viroid disease induction.

De novo sequencing and characterization of Picrorhiza kurrooa transcriptome at two temperatures showed major transcriptome adjustments

BACKGROUND

Picrorhiza kurrooa Royle ex Benth. is an endangered plant species of medicinal importance. The medicinal property is attributed to monoterpenoids picroside I and II, which are modulated by temperature. The transcriptome information of this species is limited with the availability of few hundreds of expressed sequence tags (ESTs) in the public databases. In order to gain insight into temperature mediated molecular changes, high throughput de novo transcriptome sequencing and analyses were carried out at 15 °C and 25 °C, the temperatures known to modulate picrosides content.

RESULTS

Using paired-end (PE) Illumina sequencing technology, a total of 20,593,412 and 44,229,272 PE reads were obtained after quality filtering for 15 °C and 25 °C, respectively. Available (e.g., De-Bruijn/Eulerian graph) and in-house developed bioinformatics tools were used for assembly and annotation of transcriptome. A total of 74,336 assembled transcript sequences were obtained, with an average coverage of 76.6 and average length of 439.5. Guanine-cytosine (GC) content was observed to be 44.6%, while the transcriptome exhibited abundance of trinucleotide simple sequence repeat (SSR; 45.63%) markers.Large scale expression profiling through "read per exon kilobase per million (RPKM)", showed changes in several biological processes and metabolic pathways including cytochrome P450s (CYPs), UDP-glycosyltransferases (UGTs) and those associated with picrosides biosynthesis. RPKM data were validated by reverse transcriptase-polymerase chain reaction using a set of 19 genes, wherein 11 genes behaved in accordance with the two expression methods.

CONCLUSIONS

Study generated transcriptome of P. kurrooa at two different temperatures. Large scale expression profiling through RPKM showed major transcriptome changes in response to temperature reflecting alterations in major biological processes and metabolic pathways, and provided insight of GC content and SSR markers. Analysis also identified putative CYPs and UGTs that could help in discovering the hitherto unknown genes associated with picrosides biosynthesis.