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

Comparative analysis of miRNA expression profiles in flowering and non-flowering tissue of Crocus sativus L

Crocus sativus is a valuable plant due to the presence of apocarotenoids in its stigma. Considerable work has been done in the past to understand the apocarotenoid biosynthetic pathway in saffron. However, the reports on understanding the regulation of flowering at the post-transcriptional level are meagre. The study aimed to discover the candidate miRNAs, target genes, transcription factors (TFs), and apocarotenoid biosynthetic pathway genes associated with the regulation and transition of flowering in C. sativus. In the present investigation, miRNA profiling was performed in flowering and non-flowering corms of saffron, along with expression analysis of apocarotenoid genes and transcription factors involved in the synthesis of secondary metabolites. Significant modulation in the expression of miR156, miR159, miR166, miR172, miR395, miR396, miR399, and miR408 gene families was observed. We obtained 36 known miRNAs (26 in flowering and 10 in non-flowering) and 64 novel miRNAs (40 in flowering and 24 in non-flowering) unique to specific tissues in our analysis. TFs, including CsMADS and CsMYb, showed significant modulation in expression in flowering tissue, followed by CsHB. Additionally, the miRNAs were predicted to be involved in carbohydrate metabolism, phytohormone signalling, regulation of flower development, and response to stress, cold, and defence. The comprehensive study has enhanced our understanding of the regulatory machinery comprising factors like phytohormones, abiotic stress, apocarotenoid genes, transcription factors, and miRNAs responsible for the synthesis of apocarotenoids and developmental processes during and after flowering.

An overview of floral regulatory genes in annual and perennial plants

The floral initiation in angiosperms is a complex process influenced by endogenous and exogenous signals. With this approach, we aim to provide a comprehensive review to integrate this complex floral regulatory process and summarize the regulatory genes and their functions in annuals and perennials. Seven primary paths leading to flowering have been discovered in Arabidopsis under several growth condition that include; photoperiod, ambient temperature, vernalization, gibberellins, autonomous, aging and carbohydrates. These pathways involve a series of interlinked signaling pathways that respond to both internal and external signals, such as light, temperature, hormones, and developmental cues, to coordinate the expression of genes that are involved in flower development. Among them, the photoperiodic pathway was the most important and conserved as some of the fundamental loci and mechanisms are shared even by closely related plant species. The activation of floral regulatory genes such as FLC, FT, LFY, and SOC1 that determine floral meristem identity and the transition to the flowering stage result from the merging of these pathways. Recent studies confirmed that alternative splicing, antisense RNA and epigenetic modification play crucial roles by regulating the expression of genes related to blooming. In this review, we documented recent progress in the floral transition time in annuals and perennials, with emphasis on the specific regulatory mechanisms along with the application of various molecular approaches including overexpression studies, RNA interference and Virus-induced flowering. Furthermore, the similarities and differences between annual and perennial flowering will aid significant contributions to the field by elucidating the mechanisms of perennial plant development and floral initiation regulation.

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

BACKGROUND

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

METHODS AND RESULTS

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

CONCLUSIONS

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

Two citrus KNAT-like genes, CsKN1 and CsKN2, are involved in the regulation of spring shoot development in sweet orange

Shoot-tip abortion is a very common phenomenon in some perennial woody plants and it affects the height, architecture, and branch orientation of trees; however, little is currently known about the underlying mechanisms. In this study, we identified a gene in sweet orange (Citrus sinensis) encoding a KNAT-like protein (CsKN1) and found high expression in the shoot apical meristem (SAM). Overexpression of CsKN1 in transgenic plants prolonged the vegetative growth of SAMs, whilst silencing resulted in either the loss or inhibition of SAMs. Yeast two-hybrid analysis revealed that CsKN1 interacted with another citrus KNAT-like protein (CsKN2), and overexpression of CsKN2 in lemon and tobacco caused an extreme multiple-meristem phenotype. Overexpression of CsKN1 and CsKN2 in transgenic plants resulted in the differential expression of numerous genes related to hormone biosynthesis and signaling. Yeast one-hybrid analysis revealed that the CsKN1-CsKN2 complex can bind to the promoter of citrus floral meristem gene LEAFY (CsLFY) and inhibit its expression. These results indicated that CsKN1 might prolong the vegetative growth period of SAMs by delaying flowering. In addition, an ethylene-responsive factor (CsERF) was found to bind to the CsKN1 promoter and suppresses its transcription. Overexpression of CsERF in Arabidopsis increased the contents of ethylene and reactive oxygen species, which might induce the occurrence of shoot-tip abscission. On the basis of our results, we conclude that CsKN1 and CsKN2 might work cooperatively to regulate the shoot-tip abscission process in spring shoots of sweet orange.

A potential endogenous gibberellin-mediated signaling cascade regulated floral transition in Magnolia × soulangeana 'Changchun'

The floral transition is a critical developmental switch in plants, and has profound effects on the flower production and yield. Magnolia × soulangeana 'Changchun' is known as a woody ornamental plant, which can bloom in spring and summer, respectively. In this study, anatomical observation, physiological measurement, transcriptome, and small RNA sequencing were performed to investigate potential endogenous regulatory mechanisms underlying floral transition in 'Changchun'. Transition of the shoot apical meristem from vegetative to reproductive growth occurred between late April and early May. During this specific developmental process, a total of 161,645 unigenes were identified, of which 73,257 were significantly differentially expressed, while a number of these two categories of miRNAs were 299 and 148, respectively. Further analysis of differentially expressed genes (DEGs) revealed that gibberellin signaling could regulate floral transition in 'Changchun' in a DELLA-dependent manner. In addition, prediction and analysis of miRNA targeted genes suggested that another potential molecular regulatory module was mediated by the miR172 family and other several novel miRNAs (Ms-novel_miR139, Ms-novel_miR229, and Ms-novel_miR232), with the participation of up- or down-regulating genes, including MsSVP, MsAP2, MsTOE3, MsAP1, MsGATA6, MsE2FA, and MsMDS6. Through the integrated analysis of mRNA and miRNA, our research results will facilitate the understanding of the potential molecular mechanism underlying floral transition in 'Changchun', and also provide basic experimental data for the plant germplasm resources innovation in Magnolia.

Genome-Wide Identification and Characterization of SQUAMOSA-Promoter-Binding Protein (SBP) Genes Involved in the Flowering Development of Citrus Clementina

SQUAMOSA-promoter binding protein (SBP)-box genes encode a family of plant-specific transcription factors that play vital roles in plant growth and development. In this study, 15 SBP-box genes were identified and isolated from Citrus clementina (CclSBPs), where 10 of these genes were predicted to be putative targets of Citrus clementina microRNA156 (CclmiR156). The 15 CclSBP genes could be classified into six groups based on phylogenetic analysis, diverse intron–exon structure, and motif prediction, similar to the SQUAMOSA promoter binding protein-like (SPL) gene family of Populus trichocarpa and Arabidopsis thaliana. Furthermore, CclSBPs classified into a group/subgroup have similar gene structures and conserved motifs, implying their functional redundancy. Tissue-specific expression analysis of CclSBPs demonstrated their diversified expression patterns. To further explore the potential role of CclSBPs during floral inductive water deficits, the dynamic changes of the 15 CclSBPs were investigated during floral inductive water deficits, and the results showed that some CclSBPs were associated with floral induction. Among these genes, CclSBP6 was not homologous to the Arabidopsis SBP-box gene family, and CclSBP7 was regulated by being alternatively spliced. Therefore, CclSBP6 and CclSBP7 were genetically transformed in Arabidopsis. Overexpression of the two genes changed the flowering time of Arabidopsis.

Transcriptome-wide identification and functional prediction of novel and flowering-related circular RNAs from trifoliate orange (Poncirus trifoliata L. Raf.)

A total of 558 potential circular RNAs (circRNAs) were identified in citrus, and these were analyzed and compared. One hundred seventy-six differentially expressed circRNAs were identified in two genotypes of trifoliate orange. Circular RNAs (circRNAs) play diverse roles in transcriptional control and microRNA (miRNA) function. However, little information is known about circRNAs in citrus. To identify citrus circRNAs and investigate their functional roles, high-throughput sequencing of precocious trifoliate orange (an early-flowering trifoliate orange mutant, Poncirus trifoliata L. Raf.) and its wild type was performed. A total of 558 potential circRNAs were identified by bioinformatic analysis, and 86.02% of these were sense-overlapping circRNAs. Their sequence features, alternative circularization, and other characteristics were investigated in this study. Compared with the wild type, 176 circRNAs were identified as differentially expressed circRNAs, 61 were significantly up-regulated and 115 were down-regulated in precocious trifoliate orange, indicating that they may play an important role in the early flowering process. Alternative circularization and differential expression of some circRNAs were verified by Sanger sequencing and real-time polymerase chain reaction. The functions of differentially expressed circRNAs and their host genes were predicted by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. We found that many differentially expressed circRNAs had abundant miRNA binding sites: 29 circRNAs were found to act as the 16 miRNA targets. Overall, these results will help to reveal the biological functions of circRNAs in growth and development of citrus.

Identification of Viscum album L. miRNAs and prediction of their medicinal values

MicroRNAs (miRNAs) are a class of approximately 22 nucleotides single-stranded non-coding RNA molecules that play crucial roles in gene expression. It has been reported that the plant miRNAs might enter mammalian bloodstream and have a functional role in human metabolism, indicating that miRNAs might be one of the hidden bioactive ingredients in medicinal plants. Viscum album L. (Loranthaceae, European mistletoe) has been widely used for the treatment of cancer and cardiovascular diseases, but its functional compounds have not been well characterized. We considered that miRNAs might be involved in the pharmacological activities of V. album. High-throughput Illumina sequencing was performed to identify the novel and conserved miRNAs of V. album. The putative human targets were predicted. In total, 699 conserved miRNAs and 1373 novel miRNAs have been identified from V. album. Based on the combined use of TargetScan, miRanda, PITA, and RNAhybrid methods, the intersection of 30697 potential human genes have been predicted as putative targets of 29 novel miRNAs, while 14559 putative targets were highly enriched in 33 KEGG pathways. Interestingly, these highly enriched KEGG pathways were associated with some human diseases, especially cancer, cardiovascular diseases and neurological disorders, which might explain the clinical use as well as folk medicine use of mistletoe. However, further experimental validation is necessary to confirm these human targets of mistletoe miRNAs. Additionally, target genes involved in bioactive components synthesis in V. album were predicted as well. A total of 68 miRNAs were predicted to be involved in terpenoid biosynthesis, while two miRNAs including val-miR152 and miR9738 were predicted to target viscotoxins and lectins, respectively, which increased the knowledge regarding miRNA-based regulation of terpenoid biosynthesis, lectin and viscotoxin expressions in V. album.

Small RNA sequencing reveals a role for sugarcane miRNAs and their targets in response to Sporisorium scitamineum infection

BACKGROUND

Sugarcane smut caused by Sporisorium scitamineum leads to a significant reduction in cane yield and sucrose content. MicroRNAs (miRNAs) play an important role in regulating plant responses to biotic stress. The present study was the first to use two sugarcane genotypes, YA05-179 (smut-resistant) and ROC22 (smut-susceptible), to identify differentially expressed miRNAs in sugarcane challenged with S. scitamineum by using high-throughput sequencing.

RESULTS

The predicted target gene number corresponding to known differentially expressed miRNAs in YA05-179 was less than that in ROC22, however most of them were in common. Expression of differential miRNAs under S. scitamineum challenge was mostly downregulated, with similar trends in the two varieties. Gene ontology (GO) analysis showed that the target gene classification of known miRNAs was similar to that of the newly identified miRNAs. These were mainly associated with cellular processes and metabolic processes in the biological process category, as well as combination and catalytic activity in the molecular function category. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed that these predicted target genes involved in a series of physiological and biochemical pathways or disease resistance-related physiological metabolism and signal transduction pathways, suggesting that the molecular interaction mechanism between sugarcane and S. scitamineum was a complex network system. These findings also showed certain predicted target genes of miR5671, miR5054, miR5783, miR5221, and miR6478 play roles in the mitogen-activated protein kinase (MAPK) signaling pathway, plant hormone signal transduction, and plant-pathogen interaction. Quantitative real-time PCR (qRT-PCR) analysis showed that majority of the known miRNAs and its predicted target genes followed a negatively regulated mode. Seven out of eight predicted target genes showed identical expression after 12 h treatment and reached the highest degree of matching at 48 h, indicating that the regulatory role of miRNAs on the target genes in sugarcane was maximized at 48 h after S. scitamineum challenge.

CONCLUSIONS

Taken together, our findings serve as evidence for the association of miRNA expression with the molecular mechanism underlying the pathogenesis of sugarcane smut, particularly on the significance of miRNA levels in relation to the cultivation of smut-resistant sugarcane varieties.

Identification and comparative profiling of miRNAs in herbaceous peony (Paeonia lactiflora Pall.) with red/yellow bicoloured flowers

Herbaceous peony (Paeonia lactiflora Pall.) is popular worldwide because of its gorgeous flower colour, and the yellow flower is the rarest. However, its mechanism of yellow formation is still unexplored from the post-translational level. In this study, the anatomy of the petal, cell sap pH and metal elements were investigated in bicoloured flower cultivar 'Jinhui' with red outer-petal and yellow inner-petal, and the yellow formation was influenced by the anatomy of petal, while not by the cell sap pH and metal elements. Subsequently, microRNAs sequencing (miRNA-seq) was used to identify small RNAs (sRNAs). A total of 4,172,810 and 3,565,152 specific unique sRNAs were obtained, 207 and 204 conserved miRNAs and 38 and 42 novel miRNAs were identified from red outer-petal and yellow inner-petal, respectively, which were confirmed by subcloning. Among these miRNAs, 163 conserved and 28 novel miRNAs were differentially expressed in two wheel of petals. And 5 differentially expressed miRNAs and their corresponding target genes related to yellow formation were screened, and their dynamic expression patterns confirmed that the yellow formation might be under the regulation of miR156e-3p-targeted squamosa promoter binding protein-like gene (SPL1). These results improve the understanding of miRNA regulation of the yellow formation in P. lactiflora.

Genome-wide screening and characterization of long non-coding RNAs involved in flowering development of trifoliate orange (Poncirus trifoliata L. Raf.)

Long non-coding RNAs (lncRNAs) have been demonstrated to play critical regulatory roles in post-transcriptional and transcriptional regulation in Arabidopsis. However, lncRNAs and their functional roles remain poorly characterized in woody plants, including citrus. To identify lncRNAs and investigate their role in citrus flowering, paired-end strand-specific RNA sequencing was performed for precocious trifoliate orange and its wild-type counterpart. A total of 6,584 potential lncRNAs were identified, 51.6% of which were from intergenic regions. Additionally, 555 lncRNAs were significantly up-regulated and 276 lncRNAs were down-regulated in precocious trifoliate orange, indicating that lncRNAs could be involved in the regulation of trifoliate orange flowering. Comparisons between lncRNAs and coding genes indicated that lncRNAs tend to have shorter transcripts and lower expression levels and that they display significant expression specificity. More importantly, 59 and 7 lncRNAs were identified as putative targets and target mimics of citrus miRNAs, respectively. In addition, the targets of Pt-miR156 and Pt-miR396 were confirmed using the regional amplification reverse-transcription polymerase chain reaction method. Furthermore, overexpression of Pt-miR156a1 and Pt-miR156a1 in Arabidopsis resulted in an extended juvenile phase, short siliques, and smaller leaves in transgenic plants compared with control plants. These findings provide important insight regarding citrus lncRNAs, thus enabling in-depth functional analyses.

Application of Genomic Technologies to the Breeding of Trees

The recent introduction of next generation sequencing (NGS) technologies represents a major revolution in providing new tools for identifying the genes and/or genomic intervals controlling important traits for selection in breeding programs. In perennial fruit trees with long generation times and large sizes of adult plants, the impact of these techniques is even more important. High-throughput DNA sequencing technologies have provided complete annotated sequences in many important tree species. Most of the high-throughput genotyping platforms described are being used for studies of genetic diversity and population structure. Dissection of complex traits became possible through the availability of genome sequences along with phenotypic variation data, which allow to elucidate the causative genetic differences that give rise to observed phenotypic variation. Association mapping facilitates the association between genetic markers and phenotype in unstructured and complex populations, identifying molecular markers for assisted selection and breeding. Also, genomic data provide in silico identification and characterization of genes and gene families related to important traits, enabling new tools for molecular marker assisted selection in tree breeding. Deep sequencing of transcriptomes is also a powerful tool for the analysis of precise expression levels of each gene in a sample. It consists in quantifying short cDNA reads, obtained by NGS technologies, in order to compare the entire transcriptomes between genotypes and environmental conditions. The miRNAs are non-coding short RNAs involved in the regulation of different physiological processes, which can be identified by high-throughput sequencing of RNA libraries obtained by reverse transcription of purified short RNAs, and by in silico comparison with known miRNAs from other species. All together, NGS techniques and their applications have increased the resources for plant breeding in tree species, closing the former gap of genetic tools between trees and annual species.

Over-Expression of GmGIa-Regulated Soybean miR172a Confers Early Flowering in Transgenic Arabidopsis thaliana

Flowering is a pivotal event in the life cycle of plants. miR172 has been widely confirmed to play critical roles in flowering time control by regulating its target gene expression in Arabidopsis. However, the role of its counterpart in soybean remains largely unclear. In the present study, we found that the gma-miR172a was regulated by a GIGANTEA ortholog, GmGIa, in soybean through miRNA metabolism. The expression analysis revealed that gma-miR172a has a pattern of diurnal rhythm expression and its abundance increased rapidly as plants grew until the initiation of flowering phase in soybean. One target gene of gma-miR172a, Glyma03g33470, was predicted and verified using a modified RLM 5'-RACE (RNA ligase-mediated rapid amplification of 5' cDNA ends) assay. Overexpression of gma-miR172a exhibited an early flowering phenotype and the expression of FT, AP1 and LFY were simultaneously increased in gma-miR172a-transgenic Arabidopsis plants, suggesting that the early flowering phenotype was associated with up-regulation of these genes. The overexpression of the gma-miR172a-resistant version of Glyma03g33470 weakened early flowering phenotype in the toe1 mutant of Arabidopsis. Taken together, our results suggested that gma-miR172a played an important role in GmGIa-mediated flowering by repressing Glyma03g33470, which in turn increased the expression of FT, AP1 and LFY to promote flowering in soybean.

Omics studies of citrus, grape and rosaceae fruit trees

Recent advance of bioinformatics and analytical apparatuses such as next generation DNA sequencer (NGS) and mass spectrometer (MS) has brought a big wave of comprehensive study to biology. Comprehensive study targeting all genes, transcripts (RNAs), proteins, metabolites, hormones, ions or phenotypes is called genomics, transcriptomics, proteomics, metabolomics, hormonomics, ionomics or phenomics, respectively. These omics are powerful approaches to identify key genes for important traits, to clarify events of physiological mechanisms and to reveal unknown metabolic pathways in crops. Recently, the use of omics approach has increased dramatically in fruit tree research. Although the most reported omics studies on fruit trees are transcriptomics, proteomics and metabolomics, and a few is reported on hormonomics and ionomics. In this article, we reviewed recent omics studies of major fruit trees, i.e. citrus, grapevine and rosaceae fruit trees. The effectiveness and prospects of omics in fruit tree research will as well be highlighted.

Identification of bolting-related microRNAs and their targets reveals complex miRNA-mediated flowering-time regulatory networks in radish (Raphanus sativus L.)

MicroRNAs (miRNAs) play vital regulatory roles in plant growth and development. The phase transition from vegetative growth to flowering is crucial in the life cycle of plants. To date, miRNA-mediated flowering regulatory networks remain largely unexplored in radish. In this study, two small RNA libraries from radish leaves at vegetative and reproductive stages were constructed and sequenced by Solexa sequencing. A total of 94 known miRNAs representing 21 conserved and 13 non-conserved miRNA families, and 44 potential novel miRNAs, were identified from the two libraries. In addition, 42 known and 17 novel miRNAs were significantly differentially expressed and identified as bolting-related miRNAs. RT-qPCR analysis revealed that some miRNAs exhibited tissue- or developmental stage-specific expression patterns. Moreover, 154 target transcripts were identified for 50 bolting-related miRNAs, which were predominately involved in plant development, signal transduction and transcriptional regulation. Based on the characterization of bolting-related miRNAs and their target genes, a putative schematic model of miRNA-mediated bolting and flowering regulatory network was proposed. These results could provide insights into bolting and flowering regulatory networks in radish, and facilitate dissecting the molecular mechanisms underlying bolting and flowering time regulation in vegetable crops.

Identification and Characterization of MicroRNAs in Ginkgo biloba var. epiphylla Mak

Ginkgo biloba, a dioecious plant known as a living fossil, is an ancient gymnosperm that stands distinct from other gymnosperms and angiosperms. Ginkgo biloba var. epiphylla (G. biloba var. epiphylla), with ovules borne on the leaf blade, is an unusual germplasm derived from G. biloba. MicroRNAs (miRNAs) are post-transcriptional gene regulators that play critical roles in diverse biological and metabolic processes. Currently, little is known about the miRNAs involved in the key stage of partly epiphyllous ovule germination in G. biloba var. epiphylla. Two small RNA libraries constructed from epiphyllous ovule leaves and normal leaves of G. biloba var. epiphylla were sequenced on an Illumina/Solexa platform. A total of 82 miRNA sequences belonging to 23 families and 53 putative novel miRNAs were identified in the two libraries. Differential expression analysis showed that 25 conserved and 21 novel miRNAs were differentially expressed between epiphyllous ovule leaves and normal leaves. The expression patterns of partially differentially expressed miRNAs and the transcript levels of their predicted target genes were validated by quantitative real time RT-PCR. All the expression profiles of the 21 selected miRNAs were similar to those detected by Solexa deep sequencing. Additionally, the transcript levels of almost all the putative target genes of 9 selected miRNAs were opposite to those of the corresponding miRNAs. The putative target genes of the differentially expressed miRNAs were annotated with Gene Ontology terms related to reproductive process, metabolic process and responding to stimulus. This work presents a broad range of small RNA transcriptome data obtained from epiphyllous ovule and normal leaves of G. biloba var. epiphylla, which may provide insights into the miRNA-mediated regulation in the epiphyllous ovule germination process.

Genomewide analysis of small RNAs in nonembryogenic and embryogenic tissues of citrus: microRNA- and siRNA-mediated transcript cleavage involved in somatic embryogenesis

Somatic embryogenesis (SE) is a process of somatic cells becoming dedifferentiated and generating embryos. SE has been widely used in biotechnology as a powerful way of regeneration and a model system for studying plant embryogenesis, but the controlling mechanisms of SE are far from clear. Here, we show the genomewide profiles of miRNAs/siRNAs and their target genes in nonembryogenic and embryogenic tissues of 'Valencia' sweet orange. By high-throughput sequencing (HTS) of small RNAs and RNA degradome tags, we identified 50 known and 45 novel miRNAs, 130 miniature inverted-repeat transposable elements (MITEs) derived, 94 other and 235 phased small interfering RNAs (siRNAs), as well as 203 target genes. The majority of the abundantly expressed miRNAs/siRNAs exhibit lower expression levels in embryogenic callus (EC) or during SE process than in nonembryogenic callus (NEC), which is supposed to derepress the target genes that are involved in development and stress response, thus to activate the biological processes required for cell differentiation. However, the conserved csi-miR156a/b, miR164b and 171c directed suppression of specific transcription factors (TFs) are supposed to inactivate the postembryonic growth thus to maintain normal SE. In this study, miRNA- and siRNA-mediated silencing of target genes was found under sophisticated regulation in citrus SE system; the enhancement effect of specific conserved miRNAs on SE was discussed, providing new clues for future investigation of mechanisms that control SE.

Identification and characterization of microRNAs from Chinese pollination constant non-astringent persimmon using high-throughput sequencing

BACKGROUND

microRNAs (miRNAs) have been shown to play key roles in regulating gene expression at post-transcriptional level, but miRNAs associated with natural deastringency of Chinese pollination-constant nonastringent persimmon (CPCNA) have never been identified.

RESULTS

In this study, two small RNA libraries established using 'Eshi No. 1' persimmon (Diospyros kaki Thunb.; CPCNA) fruits collected at 15 and 20 weeks after flowering (WAF) were sequenced through Solexa platform in order to identify miRNAs involved in deastringency of persimmon. A total of 6,258,487 and 7,634,169 reads were generated for the libraries at 15 and 20 WAF, respectively. Based on sequence similarity and hairpin structure prediction, 236 known miRNAs belonging to 65 miRNA families and 33 novel miRNAs were identified using persimmon transcriptome data. Sixty one of the characterized miRNAs exhibited pronounced difference in the expression levels between 15 and 20 WAF, 17 up-regulated and 44 down-regulated. Expression profiles of 12 conserved and 10 novel miRNAs were validated by stem loop qRT-PCR. A total of 198 target genes were predicted for the differentially expressed miRNAs, including several genes that have been reported to be implicated in proanthocyanidins (PAs, or called tannin) accumulation. In addition, two transcription factors, a GRF and a bHLH, were experimentally confirmed as the targets of dka-miR396 and dka-miR395, respectively.

CONCLUSIONS

Taken together, the present data unraveled several important miRNAs in persimmon. Among them, miR395p-3p and miR858b may regulate bHLH and MYB, respectively, which are influenced by SPL under the control of miR156j-5p and in turn regulate the structural genes involved in PA biosynthesis. In addition, dka-miR396g and miR2911a may regulate their target genes associated with glucosylation and insolubilization of tannin precursors. All of these miRNAs might play key roles in the regulation of (de)astringency in persimmon fruits under normal development conditions.

In silico identification of plant miRNAs in mammalian breast milk exosomes--a small step forward?

MicroRNAs (miRNAs) are a class of small RNA molecules that regulate gene expression by inhibiting the protein translation or targeting the mRNA cleavage. They play many important roles in living organism cells; however, the knowledge on miRNAs functions has become more extensive upon their identification in biological fluids and recent reports on plant-origin miRNAs abundance in human plasma and serum. Considering these findings, we performed a rigorous bioinformatics analysis of publicly available, raw data from high-throughput sequencing studies on miRNAs composition in human and porcine breast milk exosomes to identify the fraction of food-derived miRNAs. Several processing and filtering steps were applied to increase the accuracy, and to avoid false positives. Through aforementioned analysis, 35 and 17 miRNA species, belonging to 25 and 11 MIR families, were identified, respectively. In the human samples the highest abundance levels yielded the ath-miR166a, pab-miR951, ptc-miR472a and bdi-miR168, while in the porcine breast milk exosomes, the zma-miR168a, zma-miR156a and ath-miR166a have been identified in the largest amounts. The consensus prediction and annotation of potential human targets for select plant miRNAs suggest that the aforementioned molecules may interact with mRNAs coding several transcription factors, protein receptors, transporters and immune-related proteins, thus potentially influencing human organism. Taken together, the presented analysis shows proof of abundant plant miRNAs in mammal breast milk exosomes, pointing at the same time to the new possibilities arising from this discovery.

A comparison of performance of plant miRNA target prediction tools and the characterization of features for genome-wide target prediction

BACKGROUND

Deep-sequencing has enabled the identification of large numbers of miRNAs and siRNAs, making the high-throughput target identification a main limiting factor in defining their function. In plants, several tools have been developed to predict targets, majority of them being trained on Arabidopsis datasets. An extensive and systematic evaluation has not been made for their suitability for predicting targets in species other than Arabidopsis. Nor, these have not been evaluated for their suitability for high-throughput target prediction at genome level.

RESULTS

We evaluated the performance of 11 computational tools in identifying genome-wide targets in Arabidopsis and other plants with procedures that optimized score-cutoffs for estimating targets. Targetfinder was most efficient [89% 'precision' (accuracy of prediction), 97% 'recall' (sensitivity)] in predicting 'true-positive' targets in Arabidopsis miRNA-mRNA interactions. In contrast, only 46% of true positive interactions from non-Arabidopsis species were detected, indicating low 'recall' values. Score optimizations increased the 'recall' to only 70% (corresponding 'precision': 65%) for datasets of true miRNA-mRNA interactions in species other than Arabidopsis. Combining the results of Targetfinder and psRNATarget delivers high true positive coverage, whereas the intersection of psRNATarget and Tapirhybrid outputs deliver highly 'precise' predictions. The large number of 'false negative' predictions delivered from non-Arabidopsis datasets by all the available tools indicate the diversity in miRNAs-mRNA interaction features between Arabidopsis and other species. A subset of miRNA-mRNA interactions differed significantly for features in seed regions as well as the total number of matches/mismatches.

CONCLUSION

Although, many plant miRNA target prediction tools may be optimized to predict targets with high specificity in Arabidopsis, such optimized thresholds may not be suitable for many targets in non-Arabidopsis species. More importantly, non-conventional features of miRNA-mRNA interaction may exist in plants indicating alternate mode of miRNA target recognition. Incorporation of these divergent features would enable next-generation of algorithms to better identify target interactions.

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.

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

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

Identification and functional analysis of flowering related microRNAs in common wild rice (Oryza rufipogon Griff.)

Background

MicroRNAs (miRNAs) is a class of non-coding RNAs involved in post- transcriptional control of gene expression, via degradation and/or translational inhibition. Six-hundred sixty-one rice miRNAs are known that are important in plant development. However, flowering-related miRNAs have not been characterized in Oryza rufipogon Griff. It was approved by supervision department of Guangdong wild rice protection. We analyzed flowering-related miRNAs in O. rufipogon using high-throughput sequencing (deep sequencing) to understand the changes that occurred during rice domestication, and to elucidate their functions in flowering.

Results

Three O. rufipogon sRNA libraries, two vegetative stage (CWR-V1 and CWR-V2) and one flowering stage (CWR-F2) were sequenced using Illumina deep sequencing. A total of 20,156,098, 21,531,511 and 20,995,942 high quality sRNA reads were obtained from CWR-V1, CWR-V2 and CWR-F2, respectively, of which 3,448,185, 4,265,048 and 2,833,527 reads matched known miRNAs. We identified 512 known rice miRNAs in 214 miRNA families and predicted 290 new miRNAs. Targeted functional annotation, GO and KEGG pathway analyses predicted that 187 miRNAs regulate expression of flowering-related genes. Differential expression analysis of flowering-related miRNAs showed that: expression of 95 miRNAs varied significantly between the libraries, 66 are flowering-related miRNAs, such as oru-miR97, oru-miR117, oru-miR135, oru-miR137, et al. 17 are early-flowering -related miRNAs, including osa-miR160f, osa-miR164d, osa-miR167d, osa-miR169a, osa-miR172b, oru-miR4, et al., induced during the floral transition. Real-time PCR revealed the same expression patterns as deep sequencing. miRNAs targets were confirmed for cleavage by 5′-RACE in vivo, and were negatively regulated by miRNAs.

Conclusions

This is the first investigation of flowering miRNAs in wild rice. The result indicates that variation in miRNAs occurred during rice domestication and lays a foundation for further study of phase change and flowering in O. rufipogon. Complicated regulatory networks mediated by multiple miRNAs regulate the expression of flowering genes that control the induction of flowering.

High-throughput sequencing identification of novel and conserved miRNAs in the Brassica oleracea leaves

Background

Plant microRNAs are short (~21 nt) non-coding molecules that regulate gene expression by targeting the mRNA cleavage or protein translation inhibition. In this manner, they play many important roles in the cells of living organisms. One of the plant species in which the entire set of miRNAs has not been yet completely identified is Brassica oleracea var. capitata (cabbage). For this reason and for the economic and nutritional importance of this food crop, high-throughput small RNAs sequencing has been performed to discover the novel and conserved miRNAs in mature cabbage leaves.

Results

In this study, raw reads generated from three small RNA libraries were bioinformatically processed and further analyzed to select sequences homologous to known B. oleracea and other plant miRNAs. As a result of this analysis, 261 conserved miRNAs (belonging to 62 families) have been discovered. MIR169, MIR167 and MIR166 were the largest miRNA families, while the highest abundance molecules were miR167, miR166, miR168c and miR157a. Among the generated sequencing reads, miRNAs* were also found, such as the miR162c*, miR160a* and miR157a*. The unannotated tags were used in the prediction and evaluation of novel miRNAs, which resulted in the 26 potential miRNAs proposal. The expressions of 13 selected miRNAs were analyzed by northern blot hybridization. The target prediction and annotation for identified miRNAs were performed, according to which discovered molecules may target mRNAs encoding several potential proteins – e.g., transcription factors, polypeptides that regulate hormone stimuli and abiotic stress response, and molecules participating in transport and cell communication. Additionally, KEGG maps analysis suggested that the miRNAs in cabbage are involved in important processing pathways, including glycolysis, glycerolipid metabolism, flavonoid biosynthesis and oxidative phosphorylation.

Conclusions

Conclusively, for the first time, the large set of miRNAs was identified in mature cabbage leaves. Potential targets designation for these miRNAs may suggest their essential role in many plants primary biological processes. Presented study not only supplements the knowledge about B. oleracea miRNAs, but additionally it may be used in other research concerning the improvement of the cabbage cultivation.

Identification and expression analysis of microRNAs at the grain filling stage in rice(Oryza sativa L.)via deep sequencing

MicroRNAs (miRNAs) have been shown to play crucial roles in the regulation of plant development. In this study, high-throughput RNA-sequencing technology was used to identify novel miRNAs, and to reveal miRNAs expression patterns at different developmental stages during rice (Oryza sativa L.) grain filling. A total of 434 known miRNAs (380, 402, 390 and 392 at 5, 7, 12 and 17 days after fertilization, respectively.) were obtained from rice grain. The expression profiles of these identified miRNAs were analyzed and the results showed that 161 known miRNAs were differentially expressed during grain development, a high proportion of which were up-regulated from 5 to 7 days after fertilization. In addition, sixty novel miRNAs were identified, and five of these were further validated experimentally. Additional analysis showed that the predicted targets of the differentially expressed miRNAs may participate in signal transduction, carbohydrate and nitrogen metabolism, the response to stimuli and epigenetic regulation. In this study, differences were revealed in the composition and expression profiles of miRNAs among individual developmental stages during the rice grain filling process, and miRNA editing events were also observed, analyzed and validated during this process. The results provide novel insight into the dynamic profiles of miRNAs in developing rice grain and contribute to the understanding of the regulatory roles of miRNAs in grain filling.