Integrative analysis of miRNA and mRNA profiles in response to ethylene in rose petals during flower opening

Epigenetics and plant hormone dynamics: a functional and methodological perspective

Plant hormones, pivotal regulators of plant growth, development, and response to environmental cues, have recently emerged as central modulators of epigenetic processes governing gene expression and phenotypic plasticity. This review addresses the complex interplay between plant hormones and epigenetic mechanisms, highlighting the diverse methodologies that have been harnessed to decipher these intricate relationships. We present a comprehensive overview to understand how phytohormones orchestrate epigenetic modifications, shaping plant adaptation and survival strategies. Conversely, we explore how epigenetic regulators ensure hormonal balance and regulate the signalling pathways of key plant hormones. Furthermore, our investigation includes a search for novel genes that are regulated by plant hormones under the control of epigenetic processes. Our review offers a contemporary overview of the epigenetic-plant hormone crosstalk, emphasizing its significance in plant growth, development, and potential agronomical applications.

Petal size is controlled by the MYB73/TPL/HDA19-miR159-CKX6 module regulating cytokinin catabolism in Rosa hybrida

The size of plant lateral organs is determined by well-coordinated cell proliferation and cell expansion. Here, we report that miR159, an evolutionarily conserved microRNA, plays an essential role in regulating cell division in rose (Rosa hybrida) petals by modulating cytokinin catabolism. We uncover that Cytokinin Oxidase/Dehydrogenase6 (CKX6) is a target of miR159 in petals. Knocking down miR159 levels results in the accumulation of CKX6 transcripts and earlier cytokinin clearance, leading to a shortened cell division period and smaller petals. Conversely, knocking down CKX6 causes cytokinin accumulation and a prolonged developmental cell division period, mimicking the effects of exogenous cytokinin application. MYB73, a R2R3-type MYB transcription repressor, recruits a co-repressor (TOPLESS) and a histone deacetylase (HDA19) to form a suppression complex, which regulates MIR159 expression by modulating histone H3 lysine 9 acetylation levels at the MIR159 promoter. Our work sheds light on mechanisms for ensuring the correct timing of the exit from the cell division phase and thus organ size regulation by controlling cytokinin catabolism.

Identification and Profiling Analysis of microRNAs in Guava Fruit (Psidium guajava L.) and Their Role during Ripening

The guava (Psidium guajava L.) is a climacteric fruit with an accelerated post-harvest overripening. miRNAs are small RNA sequences that function as gene regulators in eukaryotes and are essential for their survival and development. In this study, miRNA libraries were constructed, sequenced and analyzed from the breaker and ripe stages of guava fruit cv. Siglo XXI. One hundred and seventy-four mature miRNA sequences from 28 miRNA families were identified. The taxonomic distribution of the guava miRNAs showed a high level of conservation among the dicotyledonous plants. Most of the predicted miRNA target genes were transcription factors and genes involved in the metabolism of phytohormones such as abscisic acid, auxins, and ethylene, as revealed through an ontology enrichment analysis. The miRNA families miR168, miR169, miR396, miR397, and miR482 were classified as being directly associated with maturation, whereas the miRNA families miR160, miR165, miR167, miR3930, miR395, miR398, and miR535 were classified as being indirectly associated. With this study, we intended to increase our knowledge and understanding of the regulatory process involved in the ripening process, thereby providing valuable information for future research on the ripening of guava fruit.

The Integrated mRNA and miRNA Approach Reveals Potential Regulators of Flowering Time in Arundina graminifolia

Orchids are among the most precious flowers in the world. Regulation of flowering time is one of the most important targets to enhance their ornamental value. The beauty of Arundina graminifolia is its year-round flowering, although the molecular mechanism of this flowering ability remains masked. Therefore, we performed a comprehensive assessment to integrate transcriptome and miRNA sequencing to disentangle the genetic regulation of flowering in this valuable species. Clustering analyses provided a set of molecular regulators of floral transition and floral morphogenesis. We mined candidate floral homeotic genes, including FCA, FPA, GI, FT, FLC, AP2, SOC1, SVP, GI, TCP, and CO, which were targeted by a variety of miRNAs. MiR11091 targeted the highest number of genes, including candidate regulators of phase transition and hormonal control. The conserved miR156-miR172 pathway of floral time regulation was evident in our data, and we found important targets of these miRNAs in the transcriptome. Moreover, endogenous hormone levels were determined to decipher the hormonal control of floral buds in A. graminifolia. The qRT-PCR analysis of floral and hormonal integrators validated the transcriptome expression. Therefore, miRNA-mediated mining of candidate genes with hormonal regulation forms the basis for comprehending the complex regulatory network of perpetual flowering in precious orchids. The findings of this study can do a great deal to broaden the breeding programs for flowering time manipulation of orchids.

Glutathione imparts stress tolerance against Alternaria brassicicola infection via miRNA mediated gene regulation

Glutathione (GSH) is well known to play a crucial role in imparting resistance against various pathogen invasions. Nevertheless, the role of GSH in regulating miRNA-mediated defense response is yet to be explored. To decipher the GSH-mediated regulation of miRNA expression during necrotrophic infection in Arabidopsis thaliana, wild-type Col-0 and AtECS1, the transgenic line exhibiting enhanced GSH content, were infected with necrotrophic pathogen Alternaria brassicicola. AtECS1 plants exhibited enhanced resistance as compared to wild-type. MiRNA next-generation sequencing (NGS) was performed to compare the miRNA expression in Col-0 and AtECS1 leaves. Under control condition, differentially expressed 96 known miRNAs and 17 novel miRNAs viz. ath-miR8167f, ath-miR1886.3, ath-miR3932b-5p, etc. were identified. However, under infected condition, 73 known and 43 novel differentially expressed miRNAs viz. ath-miR5652, ath-miR160b, ath-miR865-5p, etc. were identified. Functional annotation and enrichment analysis revealed that several miRNAs that target defense-related genes like leucine-rich repeat protein kinase, MYB transcription factors, TCP8, etc. were down regulated in the AtECS1 line, which, in turn, relieves the repression of their target gene expression, leading to resistance against infection. Together, the present investigation suggests that GSH plays a decisive role in modulating the miRNA-mediated regulation of defense-related genes during pathogen invasion.

Integrated analysis of miRNAome transcriptome and degradome reveals miRNA-target modules governing floral florescence development and senescence across early- and late-flowering genotypes in tree peony

As a candidate national flower of China, tree peony has extremely high ornamental, medicinal and oil value. However, the short florescence and rarity of early-flowering and late-flowering varieties restrict further improvement of the economic value of tree peony. Specific miRNAs and their target genes engaged in tree peony floral florescence, development and senescence remain unknown. This report presents the integrated analysis of the miRNAome, transcriptome and degradome of tree peony petals collected from blooming, initial flowering, full blooming and decay stages in early-flowering variety Paeonia ostii 'Fengdan', an early-flowering mutant line of Paeonia ostii 'Fengdan' and late-flowering variety Paeonia suffruticosa 'Lianhe'. Transcriptome analysis revealed a transcript ('psu.G.00014095') which was annotated as a xyloglucan endotransglycosylase/hydrolase precursor XTH-25 and found to be differentially expressed across flower developmental stages in Paeonia ostii 'Fengdan' and Paeonia suffruticosa 'Lianhe'. The miRNA-mRNA modules were presented significant enrichment in various pathways such as plant hormone signal transduction, indole alkaloid biosynthesis, arachidonic acid metabolism, folate biosynthesis, fatty acid elongation, and the MAPK signaling pathway. Multiple miRNA-mRNA-TF modules demonstrated the potential functions of MYB-related, bHLH, Trihelix, NAC, GRAS and HD-ZIP TF families in floral florescence, development, and senescence of tree peony. Comparative spatio-temporal expression investigation of eight floral-favored miRNA-target modules suggested that transcript 'psu.T.00024044' and microRNA mtr-miR166g-5p are involved in the floral florescence, development and senescence associated agronomic traits of tree peony. The results might accelerate the understanding of the potential regulation mechanism in regards to floral florescence, development and abscission, and supply guidance for tree peony breeding of varieties with later and longer florescence characteristics.

Genetic and Biochemical Aspects of Floral Scents in Roses

Floral scents possess high ornamental and economic values to rose production in the floricultural industry. In the past two decades, molecular bases of floral scent production have been studied in the rose as well as their genetic inheritance. Some significant achievements have been acquired, such as the comprehensive rose genome and the finding of a novel geraniol synthase in plants. In this review, we summarize the composition of floral scents in modern roses, focusing on the recent advances in the molecular mechanisms of floral scent production and emission, as well as the latest developments in molecular breeding and metabolic engineering of rose scents. It could provide useful information for both studying and improving the floral scent production in the rose.

Integrative Analysis of miRNAs and Their Targets Involved in Ray Floret Growth in Gerbera hybrida

MicroRNAs (miRNAs) are involved in regulating many aspects of plant growth and development at the post-transcriptional level. Gerbera (Gerbera hybrida) is an important ornamental crop. However, the role of miRNAs in the growth and development of gerbera is still unclear. In this study, we used high-throughput sequencing to analyze the expression profiles of miRNAs in ray floret during inflorescence opening. A total of 164 miRNAs were obtained, comprising 24 conserved miRNAs and 140 novel miRNAs. Ten conserved and 15 novel miRNAs were differentially expressed during ray floret growth, and 607 differentially expressed target genes of these differentially expressed miRNAs were identified using psRNATarget. We performed a comprehensive analysis of the expression profiles of the miRNAs and their targets. The changes in expression of five miRNAs (ghy-miR156, ghy-miR164, ghy-miRn24, ghy-miRn75 and ghy-miRn133) were inversely correlated with the changes in expression of their eight target genes. The miRNA cleavage sites in candidate target gene mRNAs were determined using 5′-RLM-RACE. Several miRNA-mRNA pairs were predicted to regulate ray floret growth and anthocyanin biosynthesis. In conclusion, the results of small RNA sequencing provide valuable information to reveal the mechanisms of miRNA-mediated ray floret growth and anthocyanin accumulation in gerbera.

Unravelling the Role of Epigenetic Modifications in Development and Reproduction of Angiosperms: A Critical Appraisal

Epigenetics are the heritable changes in gene expression patterns which occur without altering DNA sequence. These changes are reversible and do not change the sequence of the DNA but can alter the way in which the DNA sequences are read. Epigenetic modifications are induced by DNA methylation, histone modification, and RNA-mediated mechanisms which alter the gene expression, primarily at the transcriptional level. Such alterations do control genome activity through transcriptional silencing of transposable elements thereby contributing toward genome stability. Plants being sessile in nature are highly susceptible to the extremes of changing environmental conditions. This increases the likelihood of epigenetic modifications within the composite network of genes that affect the developmental changes of a plant species. Genetic and epigenetic reprogramming enhances the growth and development, imparts phenotypic plasticity, and also ensures flowering under stress conditions without changing the genotype for several generations. Epigenetic modifications hold an immense significance during the development of male and female gametophytes, fertilization, embryogenesis, fruit formation, and seed germination. In this review, we focus on the mechanism of epigenetic modifications and their dynamic role in maintaining the genomic integrity during plant development and reproduction.

Effect of aneuploidy of a non-essential chromosome on gene expression in maize

The non-essential supernumerary maize (Zea mays) B chromosome (B) has recently been shown to contain active genes and to be capable of impacting gene expression of the A chromosomes. However, the effect of the B chromosome on gene expression is still unclear. In addition, it is unknown whether the accumulation of the B chromosome has a cumulative effect on gene expression. To examine these questions, the global expression of genes, microRNAs (miRNAs), and transposable elements (TEs) of leaf tissue of maize W22 plants with 0-7 copies of the B chromosome was studied. All experimental genotypes with B chromosomes displayed a trend of upregulated gene expression for a subset of A-located genes compared to the control. Over 3000 A-located genes are significantly differentially expressed in all experimental genotypes with the B chromosome relative to the control. Modulations of these genes are largely determined by the presence rather than the copy number of the B chromosome. By contrast, the expression of most B-located genes is positively correlated with B copy number, showing a proportional gene dosage effect. The B chromosome also causes increased expression of A-located miRNAs. Differentially expressed miRNAs potentially regulate their targets in a cascade of effects. Furthermore, the varied copy number of the B chromosome leads to the differential expression of A-located and B-located TEs. The findings provide novel insights into the function and properties of the B chromosome.

Genome-Wide Identification and Comparative Profiling of MicroRNAs Reveal Flavonoid Biosynthesis in Two Contrasting Flower Color Cultivars of Tree Peony

MicroRNA (miRNA)-mediated gene regulation is involved in various physiological processes in plants. Flower color is one of the vital ornamental traits of tree peony (Paeonia suffruticosa Andr.). However, the yellow-flowered tree peony cultivars are particularly rare. To elucidate the miRNA-mediated gene regulatory mechanism underlying yellow pigmentation in tree peony, we combined pigment assessment, miRNA identification, expression analysis, and gene functional verification in two contrasting flower color cultivars "High Noon" and "Roufurong." Flavones/flavonols and anthocyanins were found to be the main contributors to the coloration of "High Noon" and "Roufurong" petals, respectively. Subsequently, miRNA analysis based on available genome data identified 9 differentially expressed miRNAs and 12 relevant target genes implicated in flavonoid biosynthesis. Their dynamic expression patterns determined the key role of mdm-miR156b-PsSPL2 module in yellow pigmentation of tree peony flowers. The sequence analysis and subcellular localization validated that PsSPL2 might function as a nuclear-localized transcription factor. Overexpression of PsSPL2 in tobacco resulted in a decrease of anthocyanin content and down-regulation of NtF3'H and NtDFR transcripts. PsSPL2-silenced petals exhibited lighter yellow color, and the contents of THC, Ap, and Ch decreased significantly. Meanwhile, expression levels of PsCHS, PsCHI, and PsF3H were significantly decreased in the petals with PsSPL2 silencing, while those of PsF3'H and PsDFR were remarkably increased. This study offers a novel insight into yellow pigmentation-related miRNA regulation network in tree peony, and further provides the valuable information on physiological changes during yellow coloring process of tree peony.

Rose (Rosa hybrida) Ethylene Responsive Factor 3 Promotes Rose Flower Senescence via Direct Activation of the Abscisic Acid Synthesis-Related 9-CIS-EPOXYCAROTENOID DIOXYGENASE Gene

During plant senescence, energy and nutrients are transferred to young leaves, fruits or seeds. However, senescence reduces flower quality, which leads to huge economic losses in flower production. Ethylene is an important factor affecting the quality of cut roses during transportation and storage. Ethylene-responsive factors (ERFs) are key nodes in ethylene signaling, but the molecular mechanism underlying ERFs regulated flower senescence is not well understood. We addressed this issue in the present study by focusing on RhERF3 from Rosa hybrida, an ERF identified in a previous transcriptome analysis of ethylene-treated rose flowers. Expression of RhERF3 was strongly induced by ethylene during rose flower senescence. Transient silencing of RhERF3 delayed flower senescence, whereas overexpression (OE) accelerated the process. RNA sequencing analysis of RhERF3 OE and pSuper vector control samples identified 13,214 differentially expressed genes that were mostly related to metabolic process and plant hormone signal transduction. Transient activation and yeast one-hybrid assays demonstrated that RhERF3 directly bound the promoter of the 9-cis-epoxycarotenoid dioxygenase (RhNCED1) gene and activated gene expression. Thus, a RhERF3/RhNCED1 axis accelerates rose flower senescence.

MicroRNAs in Woody Plants

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

Transcriptome-wide identification and characterization of microRNAs in diverse phases of wood formation in Populus trichocarpa

We applied miRNA expression profiling method to Populus trichocarpa stems of the three developmental stages, primary stem (PS), transitional stem (TS), and secondary stem (SS), to investigate miRNA species and their regulation on lignocellulosic synthesis and related processes. We obtained 892, 872, and 882 known miRNAs and 1727, 1723, and 1597 novel miRNAs, from PS, TS, and SS, respectively. Comparisons of these miRNA species among different developmental stages led to the identification of 114, 306, and 152 differentially expressed miRNAs (DE-miRNAs), which had 921, 2639, and 2042 candidate target genes (CTGs) in the three respective stages of the same order. Correlation analysis revealed 47, 439, and 71 DE-miRNA-CTG pairs of high negative correlation in PS, TS, and SS, respectively. Through biological process analysis, we finally identified 34, 6, and 76 miRNA-CTG pairs from PS, TS, and SS, respectively, and the miRNA target genes in these pairs regulate or participate lignocellulosic biosynthesis-related biological processes: cell division and differentiation, cell wall modification, secondary cell wall biosynthesis, lignification, and programmed cell death processes. This is the first report on an integrated analysis of genome-wide mRNA and miRNA profilings during multiple phases of poplar stem development. Our analysis results imply that individual miRNAs modulate secondary growth and lignocellulosic biosynthesis through regulating transcription factors and lignocellulosic biosynthetic pathway genes, resulting in more dynamic promotion, suppression, or regulatory circuits. This study advanced our understanding of many individual miRNAs and their essential, diversified roles in the dynamic regulation of secondary growth in woody tree species.

AUXIN RESPONSE FACTOR 18-HISTONE DEACETYLASE 6 module regulates floral organ identity in rose (Rosa hybrida)

The phytohormone auxin plays a pivotal role in floral meristem initiation and gynoecium development, but whether and how auxin controls floral organ identity remain largely unknown. Here, we found that auxin levels influence organ specification, and changes in auxin levels influence homeotic transformation between petals and stamens in rose (Rosa hybrida). The PIN-FORMED-LIKES (PILS) gene RhPILS1 governs auxin levels in floral buds during floral organogenesis. RhAUXIN RESPONSE FACTOR 18 (RhARF18), whose expression decreases with increasing auxin content, encodes a transcriptional repressor of the C-class gene RhAGAMOUS (RhAG), and controls stamen-petal organ specification in an auxin-dependent manner. Moreover, RhARF18 physically interacts with the histone deacetylase (HDA) RhHDA6. Silencing of RhHDA6 increases H3K9/K14 acetylation levels at the site adjacent to the RhARF18-binding site in the RhAG promoter and reduces petal number, indicating that RhARF18 might recruit RhHDA6 to the RhAG promoter to reinforce the repression of RhAG transcription. We propose a model for how auxin homeostasis controls floral organ identity via regulating transcription of RhAG.

High-throughput miRNA deep sequencing in response to drought stress in sugarcane

Drought is a major factor which reduces cane growth and productivity. In the present study, we sequenced drought susceptible (V1) and drought tolerant (V2) sugarcane varieties using high-throughput miRNA deep sequencing method to study the regulation of gene expression by miRNAs during drought stress in sugarcane. A total of 1224 conserved miRNAs which belong to 89 miRNA families were identified and 38% of the differentially regulated miRNAs were common for both varieties. Additionally 435 novel miRNAs were also identified from four small RNA libraries. We identified 145 miRNAs that were differentially expressed in susceptible variety (V1-31) and 143 miRNAs differentially expressed in the tolerant variety (V2-31). Target prediction revealed that the genes mainly encoded transcription factors, proteins, phosphatase and kinases involved in signal transduction pathways, integral component of membrane and inorganic ion transport metabolism, enzymes involved in carbohydrate transport and metabolism and drought-stress-related proteins involved in defense mechanisms. Pathway analysis of targets revealed that "General function prediction only" was the most significant pathway observed in both tolerant and susceptible genotypes followed by "signal transduction mechanisms". Functional annotation of the transcripts revealed genes like calcium-dependent protein kinase, respiratory burst oxidase, caffeic acid 3-O-methyltransferase, peroxidase, calmodulin, glutathione S-transferase and transcription factors like MYB, WRKY that are involved in drought tolerant pathways. qRT-PCR was used to verify the expression levels of miRNAs and their potential targets obtained from RNA sequencing results.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-021-02857-x.

Novel approaches on identification of conserved miRNAs for broad-spectrum Potyvirus control measures

Potyviridae comprises more than 200 ssRNA viruses, many of which have a broad host range and geographical distributions. Potyvirids (members of Potyviridae) infect several economically important plants such as saffron, cardamom, cucumber, pepper, potato, tomato, yam, etc. Cumulatively, potyvirids cause a substantial economic loss. The major bottleneck in developing an efficient antiviral strategy is that viruses quickly evade host immunity owing to their higher mutation and recombination rates. Due to this reason, the emergence of newer and improved broad-spectrum approaches to combat viral infections is essential. The use of microRNA's (miRNA) to circumvent viral infection against animal viruses has been successfully employed. Fewer studies reported the development of efficient miRNA-based antivirus resistant strategies against plant viruses and none focused on multiple virus resistance. We focused on potyviruses since studies are limited and identification of conserved miRNAs among various host plants would be an initiative to design broad-spectrum antivirus strategies. In this study, we predicted evolutionarily conserved miRNAs by BLAST searching of reported miRNAs from 15 plants against the GSS and EST sequences of banana. A total of nine miRNAs were predicted and screened in nine diverse potyvirids' hosts (Banana, Tomato, Green gram, Jasmine, Chilli, Coriander, Onion, Rose and Colocasia) belonging to eight different orders (Zingiberales, Solanales, Fabales, Lamiales, Apiales, Asperagales, Rosales and Alismatales). Results suggested that miR168 and miR162 are conserved among all the selected plants. This comprehensive study laid the foundations to design broad-spectrum antivirus resistance using miRNAs. To conclude miR168 and miR162 are conserved among many plants and play a crucial role in evading virus infection which could be used as a potential candidate for developing antiviral strategies against potyvirid infections.

Azadirachta indica MicroRNAs: Genome-Wide Identification, Target Transcript Prediction, and Expression Analyses

MicroRNAs are short, endogenous, non-coding RNAs, liable for essential regulatory function. Numerous miRNAs have been identified and studied in plants with known genomic or small RNA resources. Despite the availability of genomic and transcriptomic resources, the miRNAs have not been reported in the medicinal tree Azadirachta indica (Neem) till date. Here for the first time, we report extensive identification of miRNAs and their possible targets in A. indica which might help to unravel their therapeutic potential. A comprehensive search of miRNAs in the A. indica genome by C-mii tool was performed. Overall, 123 miRNAs classified into 63 families and their stem-loop hairpin structures were predicted. The size of the A. indica (ain)-miRNAs ranged between 19 and 23 nt in length, and their corresponding ain-miRNA precursor sequence MFEI value averaged as -1.147 kcal/mol. The targets of ain-miRNAs were predicted in A. indica as well as Arabidopsis thaliana plant. The gene ontology (GO) annotation revealed the involvement of ain-miRNA targets in developmental processes, transport, stress, and metabolic processes including secondary metabolism. Stem-loop qRT-PCR was carried out for 25 randomly selected ain-miRNAs and differential expression patterns were observed in different A. indica tissues. Expression of miRNAs and its targets shows negative correlation in a dependent manner.

The regulatory framework of developmentally programmed cell death in floral organs: A review

Developmentally programmed cell death (dPCD) is a tightly controlled biological process. In recent years, vital roles of dPCD on regulating floral organ growth and development have been reported. It is well known that flower is an essential organ for reproduction and a turning point of plants' life cycle. Hence, uncovering the complex molecular networks which regulates dPCD processes in floral organs is utmost important. So far, our understanding of dPCD on floral organ growth and development is just starting. Herein, we summarize the important factors that involved in the tapetal degeneration, pollen tube rupture, receptive synergid cell death, nucellar degradation, and antipodal cell degradation. Meanwhile, the known factors that involved in transmitting tract formation and self-incompatibility-induced PCD were also introduced. Furthermore, the genes that associated with anther dehiscence and petal senescence and abscission were reviewed as well. The functions of various types of factors involved in floral dPCD processes are highlighted principally. The regulatory panorama described here can provide us some insights about flower-specific dPCD process.

Research Tools for the Functional Genomics of Plant miRNAs During Zygotic and Somatic Embryogenesis

During early plant embryogenesis, some of the most fundamental decisions on fate and identity are taken making it a fascinating process to study. It is no surprise that higher plant embryogenesis was intensively analysed during the last century, while somatic embryogenesis is probably the most studied regeneration model. Encoded by the MIRNA, short, single-stranded, non-coding miRNAs, are commonly present in all Eukaryotic genomes and are involved in the regulation of the gene expression during the essential developmental processes such as plant morphogenesis, hormone signaling, and developmental phase transition. During the last few years dedicated to miRNAs, analytical methods and tools have been developed, which have afforded new opportunities in functional analyses of plant miRNAs, including (i) databases for in silico analysis; (ii) miRNAs detection and expression approaches; (iii) reporter and sensor lines for a spatio-temporal analysis of the miRNA-target interactions; (iv) in situ hybridisation protocols; (v) artificial miRNAs; (vi) MIM and STTM lines to inhibit miRNA activity, and (vii) the target genes resistant to miRNA. Here, we attempted to summarise the toolbox for functional analysis of miRNAs during plant embryogenesis. In addition to characterising the described tools/methods, examples of the applications have been presented.

RhMYB108, an R2R3-MYB transcription factor, is involved in ethylene- and JA-induced petal senescence in rose plants

Rose (Rosa hybrida) plants are major ornamental species worldwide, and their commercial value greatly depends on their open flowers, as both the quality of fully open petals and long vase life are important. Petal senescence can be started and accelerated by various hormone signals, and ethylene is considered an accelerator of petal senescence in rose. To date, however, the underlying mechanism of signaling crosstalk between ethylene and other hormones such as JA in petal senescence remains largely unknown. Here, we isolated RhMYB108, an R2R3-MYB transcription factor, which is highly expressed in senescing petals as well as in petals treated with exogenous ethylene and JA. Applications of exogenous ethylene and JA markedly accelerated petal senescence, while the process was delayed in response to applications of 1-MCP, an ethylene action inhibitor. In addition, silencing of RhMYB108 alter the expression of SAGs such as RhNAC029, RhNAC053, RhNAC092, RhSAG12, and RhSAG113, and finally block ethylene- and JA-induced petal senescence. Furthermore, RhMYB108 was identified to target the promoters of RhNAC053, RhNAC092, and RhSAG113. Our results reveal a model in which RhMYB108 functions as a receptor of ethylene and JA signals to modulate the onset of petal senescence by targeting and enhancing senescence-associated gene expression.

Integration of Transcriptomes, Small RNAs, and Degradome Sequencing to Identify Putative miRNAs and Their Targets Related to Eu-Rubber Biosynthesis in Eucommia ulmoides

Eucommia ulmoides has attracted much attention as a valuable natural rubber (Eu-rubber) production tree. As a strategic material, Eu-rubber plays a vital role in general and defence industries. However, the study of Eu-rubber biosynthesis at a molecular level is scarce, and the regulatory network between microRNAs (miRNAs) and messenger RNAs (mRNAs) in Eu-rubber biosynthesis has not been assessed. In this study, we comprehensively analyzed the transcriptomes, small RNAs (sRNAs) and degradome to reveal the regulatory network of Eu-rubber biosynthesis in E. ulmoides. A total of 82,065 unigenes and 221 miRNAs were identified using high-throughput sequencing; 20,815 targets were predicted using psRNATarget software. Of these targets, 779 miRNA-target pairs were identified via degradome sequencing. Thirty-one miRNAs were differentially expressed; 22 targets of 34 miRNAs were annotated in the terpenoid backbone biosynthesis pathway (ko00900) based on the Kyoto Encyclopedia of Genes and Genomes (KEGG). These miRNAs were putatively related to Eu-rubber biosynthesis. A regulatory network was constructed according to the expression profiles of miRNAs and their targets. These results provide a comprehensive analysis of transcriptomics, sRNAs and degradome to reveal the Eu-rubber accumulation, and provide new insights into genetic engineering techniques which may improve the content of Eu-rubber in E. ulmoides.

Hybrid-Transcriptome Sequencing and Associated Metabolite Analysis Reveal Putative Genes Involved in Flower Color Difference in Rose Mutants

Gene mutation is a common phenomenon in nature that often leads to phenotype differences, such as the variations in flower color that frequently occur in roses. With the aim of revealing the genomic information and inner mechanisms, the differences in the levels of both transcription and secondary metabolism between a pair of natural rose mutants were investigated by using hybrid RNA-sequencing and metabolite analysis. Metabolite analysis showed that glycosylated derivatives of pelargonidin, e.g., pelargonidin 3,5 diglucoside and pelargonidin 3-glucoside, which were not detected in white flowers (Rosa 'Whilte Mrago Koster'), constituted the major pigments in pink flowers. Conversely, the flavonol contents of petal, such as kaempferol-3-glucoside, quercetin 3-glucoside, and rutin, were higher in white flowers. Hybrid RNA-sequencing obtained a total of 107,280 full-length transcripts in rose petal which were annotated in major databases. Differentially expressed gene (DEG) analysis showed that the expression of genes involved in the flavonoid biosynthesis pathway was significantly different, e.g., CHS, FLS, DFR, LDOX, which was verified by qRT-PCR during flowering. Additionally, two MYB transcription factors were found and named RmMYBAN2 and RmMYBPA1, and their expression patterns during flowering were also analyzed. These findings indicate that these genes may be involved in the flower color difference in the rose mutants, and competition between anthocyanin and flavonol biosynthesis is a primary cause of flower color variation, with its regulation reflected by transcriptional and secondary metabolite levels.

Transcriptome-wide identification and characterization of microRNAs responsive to phosphate starvation in Populus tomentosa

miRNAs (microRNAs) are ~ 21-nt non-coding small RNAs (sRNAs) that play crucial regulatory roles in plant biotic and abiotic stress responses. Phosphorus (Pi) deficiency constrains plant growth and reduces yields worldwide. To identify tree miRNAs and evaluate their functions in the response to low Pi, we identified 261 known and 31 candidate novel miRNA families from three sRNA libraries constructed from Populus tomentosa subjected to sufficient or Pi deficiency condition or to restoration of a sufficient Pi level after Pi deficiency. Pi deficiency resulted in significant changes in the abundance of TPM (transcript per million) of 65 known and 3 novel miRNAs. Interestingly, four miRNAs responsive to low N-miR167, miR394, miR171, and miR857-were found to be involved in the response to low Pi. Thirty-five known and one novel miRNAs responded dynamically to Pi fluctuations, suggesting their involvement in the response to Pi deficiency. miRNA clusters comprising 36 miRNAs were identified in 10 chromosomes. Intriguingly, nine pairs of sense and antisense miRNAs transcribed from the same loci were detected in P. tomentosa, which is the first such report in woody plants. Moreover, target genes of the known miRNAs and novel miRNA candidates with significantly changed abundance were predicted, and their functions were annotated. Degradome sequencing supported the identified targets of miRNAs in P. tomentosa. These findings will enhance our understanding of universal and specific molecular regulatory mechanisms of trees under nutrition stress and may facilitate improvement of the Pi utilization efficiency of woody plants.

Analysis of microRNAs, phased small interfering RNAs and their potential targets in Rosarugosa Thunb

BACKGROUND

MicroRNAs (miRNAs) are small non-coding RNAs that play important roles by regulating other genes. Rosa rugosa Thunb. is an important ornamental and edible plant, yet there are only a few studies on the miRNAs and their functions in R. rugosa.

RESULTS

We sequenced 10 samll RNA profiles from the roots, petals, pollens, stamens, and leaves and 4 RNA-seq profiles in leaves and petals to analysis miRNA, phasiRNAs and mRNAs in R. rugosa. In addition, we acquired a degradome sequencing profile from leaf of R. rugosa to identify miRNA and phasiRNA targets using the SeqTar algorithm. We have identified 321 conserved miRNA homologs including primary transcripts for 25 conserved miRNAs, and 22 novel miRNAs. We identified 592 putative targets of the conserved miRNAs or tasiRNAs that showed significant accumulations of degradome reads. We found differential expression patterns of conserved miRNAs in five different tissues of R. rugosa. We identified three hundred and thirty nine 21 nucleotide (nt) PHAS loci, and forty nine 24 nt PHAS loci, respectively. Our results suggest that miR482 triggers generations of phasiRNAs by targeting nucleotide-binding, leucine-rich repeat (NB-LRR) disease resistance genes in R. rugosa. Our results also suggest that the deregulated genes in leaves and petals are significantly enriched in GO terms and KEGG pathways related to metabolic processes and photosynthesis.

CONCLUSIONS

These results significantly enhanced our knowledge of the miRNAs and phasiRNAs, as well as their potential functions, in R. rugosa.

Identification of Ethylene Responsive miRNAs and Their Targets from Newly Harvested Banana Fruits Using High-Throughput Sequencing

The roles of microRNAs (miRNAs) related to ethylene response in banana fruits remain unknown because many miRNAs are differentially expressed as the fruit ripens, making the identification of ethylene-responsive miRNAs difficult. Using newly harvested banana fruits (within 5 h after harvest) as material, we found that these fruit did not ripen when treated with 5 μL/L of ethylene for 12 h at 22 °C. Two miRNA libraries were generated from newly harvested banana fruits with and without ethylene treatment and sequenced. In total, 128 known miRNAs belonging to 42 miRNA families were obtained, and 12 novel miRNAs were identified. Among them, 22 were differentially expressed in response to ethylene treatment, among which 6 known miRNAs and their putative targets were validated using qRT-PCR. These putative targets encoded proteins including GATA, ARF, DLC, and AGO, etc. KEGG and GO analyses showed that miRNAs differentially expressed in response to ethylene mainly function in the molecular and biological processes.

Analysis of microRNA reveals cleistogamous and chasmogamous floret divergence in dimorphic plant

Cleistogenes songorica, a grass species that exhibits two spatially different type of inflorescence, chastogamy (CH), flowers localized at the top, and cleistogamy (CL) flowers embedded in leaf sheath. This study aimed at dissecting reasons underlying these distinct floral development patterns at morphological and microRNA level. Phenotyping for CH and CL was conducted and four small RNA libraries were constructed from the CH and CL flowers for high-throughput sequencing to identify the differentiated miRNAs. As results, spikelet, stigma, anther, lemma and lodicule length of CH flowers were found larger than that of CL, and so was seed setting. Also, 17 flower-related differential expression miRNAs were identified which were associated with floral organ development and morphogenesis, and the flower development. Further results showed that miR159a.1-CL3996.Contig2 pair was related to anther development, miR156a-5p-CL1954.Contig2 was linked to response to high light intensity, miR408-3p/miR408d-Unigene429 was related to pollination and Unigene429 positively regulated flower development. To our knowledge, this is the first study on differential miRNA accumulation between CH and CL flowers and our study serves as a foundation to the future elucidation of regulatory mechanisms of miRNAs in the divergent development of CL and CH flowers in a single plant.

Identification and characterization of microRNAs in tree peony during chilling induced dormancy release by high-throughput sequencing

Tree peony, one of the most valuable horticultural and medicinal plants in the world, has to go through winter to break dormancy. Growing studies from molecular aspects on dormancy release process have been reported, but inadequate study has been done on miRNA-guided regulation in tree peony. In this study, high-throughput sequencing was employed to identify and characterize miRNAs in three libraries (6 d, 18 d and 24 d chilling treatments). There were 7,122, 10,076 and 9,097 unique miRNA sequences belonging to 52, 87 and 68 miRNA families, respectively. A total of 32 conserved miRNAs and 17 putative novel miRNAs were identified during dormancy release. There were 771 unigenes as potential targets of 62 miRNA families. Total 112 known miRNAs were differentially expressed, of which 55 miRNAs were shared among three libraries and 28 miRNAs were only found in 18 d chilling duration library. The expression patterns of 15 conserved miRNAs were validated and classified into four types by RT-qPCR. Combining with our microarray data under same treatments, five miRNAs (miR156k, miR159a, miR167a, miR169a and miR172a) were inversely correlated to those of their target genes. Our results would provide new molecular basis about dormancy release in tree peony.

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

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

Petal senescence: a hormone view

Flowers are highly complex organs that have evolved to enhance the reproductive success of angiosperms. As a key component of flowers, petals play a vital role in attracting pollinators and ensuring successful pollination. Having fulfilled this function, petals senesce through a process that involves many physiological and biochemical changes that also occur during leaf senescence. However, petal senescence is distinct, due to the abundance of secondary metabolites in petals and the fact that petal senescence is irreversible. Various phytohormones are involved in regulating petal senescence, and are thought to act both synergistically and antagonistically. In this regard, there appears to be developmental point during which such regulatory signals are sensed and senescence is initiated. Here, we review current understanding of petal senescence, and discuss associated regulatory mechanisms involving hormone interactions and epigenetic regulation.

Small RNAs, emerging regulators critical for the development of horticultural traits

Small RNAs (sRNAs) have been recently recognized as key genetic and epigenetic regulators in various organisms, ranging from the modification of DNA and histone methylations to the modulation of the abundance of coding or non-coding RNAs. In plants, major regulatory sRNAs are classified as respective microRNA (miRNA) and small interfering RNA (siRNA) species, with the former primarily engaging in posttranscriptional regulation while the latter in transcriptional one. Many of these characterized sRNAs are involved in regulation of diverse biological programs, processes, and pathways in response to developmental cues, environmental signals/stresses, pathogen infection, and pest attacks. Recently, sRNAs-mediated regulations have also been extensively investigated in horticultural plants, with many novel mechanisms unveiled, which display far more mechanistic complexity and unique regulatory features compared to those studied in model species. Here, we review the recent progress of sRNA research in horticultural plants, with emphasis on mechanistic aspects as well as their relevance to trait regulation. Given that major and pioneered sRNA research has been carried out in the model and other plants, we also discuss ongoing sRNA research on these plants. Because miRNAs and phased siRNAs (phasiRNAs) are the most studied sRNA regulators, this review focuses on their biogenesis, conservation, function, and targeted genes and traits as well as the mechanistic relation between them, aiming at providing readers comprehensive information instrumental for future sRNA research in horticulture crops.

Dynamic mRNA and miRNA expression analysis in response to hypoxia and reoxygenation in the blunt snout bream (Megalobrama amblycephala)

Adaptation to hypoxia is a complex process involving various pathways and regulation mechanisms. A better understanding of the genetic influence on these mechanisms could permit selection for hypoxia-sensitive fish. To aid this understanding, an integrated analysis of miRNA and mRNA expression was performed in Megalobrama amblycephala under four acute hypoxia and reoxygenation stages. A number of significantly differentially-expressed miRNAs and genes associated with oxidative stress were identified, and their functional characteristics were revealed by GO function and KEGG pathway analysis. They were found to be involved in HIF-1 pathways known to affect energy metabolism and apoptosis. MiRNA-mRNA interaction pairs were detected from comparison of expression between the four different stages. The function annotation results also showed that many miRNA-mRNA interaction pairs were likely to be involved in regulating hypoxia stress. As a unique resource for gene expression and regulation during hypoxia and reoxygenation, this study could provide a starting point for further studies to better understand the genetic background of hypoxia stress.

Uncovering Male Fertility Transition Responsive miRNA in a Wheat Photo-Thermosensitive Genic Male Sterile Line by Deep Sequencing and Degradome Analysis

MicroRNAs (miRNAs) are endogenous small RNAs which play important negative regulatory roles at both the transcriptional and post-transcriptional levels in plants. Wheat is the most commonly cultivated plant species worldwide. In this study, RNA-seq analysis was used to examine the expression profiles of miRNA in the spikelets of photo-thermosenisitive genic male sterile (PTGMS) wheat line BS366 during male fertility transition. Through mapping on their corresponding precursors, 917-7,762 novel miRNAs were found in six libraries. Six novel miRNAs were selected for examination of their secondary structures and confirmation by stem-loop RT-PCR. In a differential expression analysis, 20, 22, and 58 known miRNAs exhibited significant differential expression between developmental stages 1 (secondary sporogenous cells had formed), 2 (all cells layers were present and mitosis had ceased), and 3 (meiotic division stage), respectively, of fertile and sterile plants. Some of these differential expressed miRNAs, such as tae-miR156, tae-miR164, tae-miR171, and tae-miR172, were shown to be associated with their targets. These targets were previously reported to be related to pollen development and/or male sterility, indicating that these miRNAs and their targets may be involved in the regulation of male fertility transition in the PTGMS wheat line BS366. Furthermore, target genes of miRNA cleavage sites were validated by degradome sequencing. In this study, a possible signal model for the miRNA-mediated signaling pathway during the process of male fertility transition in the PTGMS wheat line BS366 was developed. This study provides a new perspective for understanding the roles of miRNAs in male fertility in PTGMS lines of wheat.

Plant hormone signaling in flowering: An epigenetic point of view

Reproduction is one of the most important phases in an organism's lifecycle. In the case of angiosperm plants, flowering provides the major developmental transition from the vegetative to the reproductive stage, and requires genetic and epigenetic reprogramming to ensure the success of seed production. Flowering is regulated by a complex network of genes that integrate multiple environmental cues and endogenous signals so that flowering occurs at the right time; hormone regulation, signaling and homeostasis are very important in this process. Working alone or in combination, hormones are able to promote flowering by epigenetic regulation. Some plant hormones, such as gibberellins, jasmonic acid, abscisic acid and auxins, have important effects on chromatin compaction mediated by DNA methylation and histone posttranslational modifications, which hints at the role that epigenetic regulation may play in flowering through hormone action. miRNAs have been viewed as acting independently from DNA methylation and histone modification, ignoring their potential to interact with hormone signaling - including the signaling of auxins, gibberellins, ethylene, jasmonic acid, salicylic acid and others - to regulate flowering. Therefore, in this review we examine new findings about interactions between epigenetic mechanisms and key players in hormone signaling to coordinate flowering.

Genetic control of flowering time in woody plants: Roses as an emerging model

Genetic control of the timing of flowering in woody plants is complex and has yet to be adequately investigated due to their long life-cycle and difficulties in genetic modification. Studies in Populus, one of the best woody plant models, have revealed a highly conserved genetic network for flowering timing in annuals. However, traits like continuous flowering cannot be addressed with Populus. Roses and strawberries have relatively small, diploid genomes and feature enormous natural variation. With the development of new genetic populations and genomic tools, roses and strawberries have become good models for studying the molecular mechanisms underpinning the regulation of flowering in woody plants. Here, we review findings on the molecular and genetic factors controlling continuous flowering in roses and woodland strawberries. Natural variation at TFL1 orthologous genes in both roses and strawberries seems be the key plausible factor that regulates continuous flowering. However, recent efforts suggest that a two-recessive-loci model may explain the controlling of continuous flowering in roses. We propose that epigenetic factors, including non-coding RNAs or chromatin-related factors, might also play a role. Insights into the genetic control of flowering time variation in roses should benefit the development of new germplasm for woody crops and shed light on the molecular genetic bases for the production and maintenance of plant biodiversity.

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.

Differentially Expressed microRNAs and Target Genes Associated with Plastic Internode Elongation in Alternanthera philoxeroides in Contrasting Hydrological Habitats

Phenotypic plasticity is crucial for plants to survive in changing environments. Discovering microRNAs, identifying their targets and further inferring microRNA functions in mediating plastic developmental responses to environmental changes have been a critical strategy for understanding the underlying molecular mechanisms of phenotypic plasticity. In this study, the dynamic expression patterns of microRNAs under contrasting hydrological habitats in the amphibious species Alternanthera philoxeroides were identified by time course expression profiling using high-throughput sequencing technology. A total of 128 known and 18 novel microRNAs were found to be differentially expressed under contrasting hydrological habitats. The microRNA:mRNA pairs potentially associated with plastic internode elongation were identified by integrative analysis of microRNA and mRNA expression profiles, and were validated by qRT-PCR and 5' RLM-RACE. The results showed that both the universal microRNAs conserved across different plants and the unique microRNAs novelly identified in A. philoxeroides were involved in the responses to varied water regimes. The results also showed that most of the differentially expressed microRNAs were transiently up-/down-regulated at certain time points during the treatments. The fine-scale temporal changes in microRNA expression highlighted the importance of time-series sampling in identifying stress-responsive microRNAs and analyzing their role in stress response/tolerance.

Characterization of miRNAs responsive to exogenous ethylene in grapevine berries at whole genome level

MicroRNAs (miRNAs) are critical regulators of various biological and metabolic processes of plants. Numerous miRNAs and their functions have been identified and analyzed in many plants. However, till now, the involvement of miRNAs in the response of grapevine berries to ethylene has not been reported yet. Here, Solexa technology was employed to deeply sequence small RNA libraries constructed from grapevine berries treated with and without ethylene. A total of 124 known and 78 novel miRNAs were identified. Among these miRNAs, 162 miRNAs were clearly responsive to ethylene, with 55 downregulated, 59 upregulated, and 14 unchanged miRNAs detected only in the control. The other 35 miRNAs responsive to ethylene were induced by ethylene and detected only in the ethylene-treated grapevine materials. Expression analysis of 27 conserved and 26 novel miRNAs revealed that 13 conserved and 18 novel ones were regulated by ethylene during the whole development of grapevine berries. High-throughput sequencing and qRT-PCR assays revealed consistent results on the expression results of ethylene-responsive miRNAs. Moreover, 90 target genes for 34 novel miRNAs were predicted, most of which were involved in responses to various stresses, especially like exogenous ethylene treatment. The identified miRNAs may be mainly involved in grapevine berry development and response to various environmental conditions.

Dynamic mRNA and miRNA expression analysis in response to intermuscular bone development of blunt snout bream (Megalobrama amblycephala)

Intermuscular bone (IB), which occurs only in the myosepta of lower teleosts, is attracting more attention because they are difficult to remove and make the fish unpleasant to eat. By gaining a better understanding of the genetic regulation of IB development, an integrated analysis of miRNAs and mRNAs expression profiling was performed on Megalobrama amblycephala. Four key development stages were selected for transcriptome and small RNA sequencing. A number of significantly differentially expressed miRNAs/genes associated with bone formation and differentiation were identified and the functional characteristics of these miRNAs/genes were revealed by GO function and KEGG pathway analysis. These were involved in TGF-β, ERK and osteoclast differentiation pathways known in the literature to affect bone formation and differentiation. MiRNA-mRNA interaction pairs were detected from comparison of expression between different stages. The function annotation results also showed that many miRNA-mRNA interaction pairs were likely to be involved in regulating bone development and differentiation. A negative regulation effect of two miRNAs was verified through dual luciferase reporter assay. As a unique public resource for gene expression and regulation during the IB development, this study is expected to provide forwards ideas and resources for further biological researches to understand the IBs’ development.

Deep-sequence profiling of miRNAs and their target prediction in Monotropa hypopitys

Myco-heterotroph Monotropa hypopitys is a widely spread perennial herb used to study symbiotic interactions and physiological mechanisms underlying the development of non-photosynthetic plant. Here, we performed, for the first time, transcriptome-wide characterization of M. hypopitys miRNA profile using high throughput Illumina sequencing. As a result of small RNA library sequencing and bioinformatic analysis, we identified 55 members belonging to 40 families of known miRNAs and 17 putative novel miRNAs unique for M. hypopitys. Computational screening revealed 206 potential mRNA targets for known miRNAs and 31 potential mRNA targets for novel miRNAs. The predicted target genes were described in Gene Ontology terms and were found to be involved in a broad range of metabolic and regulatory pathways. The identification of novel M. hypopitys-specific miRNAs, some with few target genes and low abundances, suggests their recent evolutionary origin and participation in highly specialized regulatory mechanisms fundamental for non-photosynthetic biology of M. hypopitys. This global analysis of miRNAs and their potential targets in M. hypopitys provides a framework for further investigation of miRNA role in the evolution and establishment of non-photosynthetic myco-heterotrophs.

Functions and Regulation of Programmed Cell Death in Plant Development

Programmed cell death (PCD) is a collective term for diverse processes causing an actively induced, tightly controlled cellular suicide. PCD has a multitude of functions in the development and health of multicellular organisms. In comparison to intensively studied forms of animal PCD such as apoptosis, our knowledge of the regulation of PCD in plants remains limited. Despite the importance of PCD in plant development and as a response to biotic and abiotic stresses, the complex molecular networks controlling different forms of plant PCD are only just beginning to emerge. With this review, we provide an update on the considerable progress that has been made over the last decade in our understanding of PCD as an inherent part of plant development. We highlight both functions of developmental PCD and central aspects of its molecular regulation.

Integration of small RNAs, degradome and transcriptome sequencing in hyperaccumulator Sedum alfredii uncovers a complex regulatory network and provides insights into cadmium phytoremediation

The hyperaccumulating ecotype of Sedum alfredii Hance is a cadmium (Cd)/zinc/lead co-hyperaccumulating species of Crassulaceae. It is a promising phytoremediation candidate accumulating substantial heavy metal ions without obvious signs of poisoning. However, few studies have focused on the regulatory roles of miRNAs and their targets in the hyperaccumulating ecotype of S. alfredii. Here, we combined analyses of the transcriptomics, sRNAs and the degradome to generate a comprehensive resource focused on identifying key regulatory miRNA-target circuits under Cd stress. A total of 87 721 unigenes and 356 miRNAs were identified by deep sequencing, and 79 miRNAs were differentially expressed under Cd stress. Furthermore, 754 target genes of 194 miRNAs were validated by degradome sequencing. A gene ontology (GO) enrichment analysis of differential miRNA targets revealed that auxin, redox-related secondary metabolism and metal transport pathways responded to Cd stress. An integrated analysis uncovered 39 pairs of miRNA targets that displayed negatively correlated expression profiles. Ten miRNA-target pairs also exhibited negative correlations according to a real-time quantitative PCR analysis. Moreover, a coexpression regulatory network was constructed based on profiles of differentially expressed genes. Two hub genes, ARF4 (auxin response factor 4) and AAP3 (amino acid permease 3), which might play central roles in the regulation of Cd-responsive genes, were uncovered. These results suggest that comprehensive analyses of the transcriptomics, sRNAs and the degradome provided a useful platform for investigating Cd hyperaccumulation in S. alfredii, and may provide new insights into the genetic engineering of phytoremediation.

Identification of chilling-responsive microRNAs and their targets in vegetable soybean (Glycine max L.)

Chilling stress is a major factor limiting the yield and quality of vegetable soybean (Glycine max L.) on a global scale. In the present study, systematic identification and functional analysis of miRNAs under chilling stress were carried out to clarify the molecular mechanism of chilling resistance. Two independent small RNA libraries from leaves of soybean were constructed and sequenced with the high-throughput Illumina Solexa system. A total of 434 known miRNAs and 3 novel miRNAs were identified. Thirty-five miRNAs were verified by qRT-PCR analysis. Furthermore, their gene targets were identified via high-throughput degradome sequencing. A total of 898 transcripts were targeted by 54 miRNA families attributed to five categories. More importantly, we identified 51 miRNAs differentially expressed between chilling stress and control conditions. The targets of these miRNAs were enriched in oxidation-reduction, signal transduction, and metabolic process functional categories. Our qRT-PCR analysis confirmed a negative relationship among the miRNAs and their targets under chilling stress. Our work thus provides comprehensive molecular evidence supporting the involvement of miRNAs in chilling-stress responses in vegetable soybean.

High throughput sequencing of small RNAs transcriptomes in two Crassostrea oysters identifies microRNAs involved in osmotic stress response

Increasing evidence suggests that microRNAs post-transcriptionally regulate gene expression and are involved in responses to biotic and abiotic stress. However, the role of miRNAs involved in osmotic plasticity remains largely unknown in marine bivalves. In the present study, we performed low salinity challenge with two Crassostrea species (C. gigas and C. hongkongensis), and conducted high-throughput sequencing of four small RNA libraries constructed from the gill tissues. A total of 202 and 87 miRNAs were identified from C. gigas and C. hongkongensis, respectively. Six miRNAs in C. gigas and two in C. hongkongensis were differentially expressed in response to osmotic stress. The expression profiles of these eight miRNAs were validated by qRT-PCR. Based on GO enrichment and KEGG pathway analysis, genes associated with microtubule-based process and cellular component movement were enriched in both species. In addition, five miRNA-mRNA interaction pairs that showed opposite expression patterns were identified in the C. hongkongensis, Differential expression analysis identified the miRNAs that play important regulatory roles in response to low salinity stress, providing insights into molecular mechanisms that are essential for salinity tolerance in marine bivalves.

Small RNA and degradome sequencing reveals microRNAs and their targets involved in tomato pedicel abscission

We constructed small RNA and degradome sequencing libraries to identify miRNAs and targets involved in tomato pedicel abscission, and confirmed their roles via quantitative real-time PCR. MicroRNAs (miRNAs) are endogenous small RNAs which play crucial negatively regulatory roles at both the transcriptional and post-transcriptional levels in plants; however, limited knowledge is available on the expression profiles of miRNAs and their target genes during tomato pedicel abscission. Taking advantage of small RNA (sRNA) and degradome sequencing technology, a total of 56 known and 11 novel candidate miRNAs targeting 223 mRNA genes were confirmed during pedicel abscission. Gene ontology annotation and KEGG pathway analysis showed that these target genes were significantly enriched in intracellular, membrane-bounded organelle-related biological processes as well as in metabolic, plant-pathogen interaction and hormone signaling pathways. We screened 17 miRNA/target pairs for further analysis and performed quantitative real-time PCR to identify the roles. Cluster analysis of selected miRNAs revealed that the expression profiles of miRNAs varied in different stages of abscission and could be impacted by ethylene treatment. In the present study, the correlations between miRNAs and targets suggested a complex regulatory network of miRNA-mediated target interaction during pedicel abscission. Additionally, the expression profiles of miRNAs and their targets changed by ethylene might be a considerable reason why ethylene promotes pedicel abscission. Our study provides new insights into the expression and regulatory profiles of miRNAs during tomato pedicel abscission.

Identification and expression profiles of sRNAs and their biogenesis and action-related genes in male and female cones of Pinus tabuliformis

Background

Small RNA (sRNA) play pivotal roles in reproductive development, and their biogenesis and action mechanisms are well characterised in angiosperm plants; however, corresponding studies in conifers are very limited. To improve our understanding of the roles of sRNA pathways in the reproductive development of conifers, the genes associated with sRNA biogenesis and action pathways were identified and analysed, and sRNA sequencing and parallel analysis of RNA ends (PARE) were performed in male and female cones of the Chinese pine (Pinus tabuliformis).

Results

Based on high-quality reference transcriptomic sequences, 21 high-confidence homologues involved in sRNA biogenesis and action in P. tabuliformis were identified, including two different DCL3 genes and one AGO4 gene. More than 75 % of genes involved in sRNA biogenesis and action have higher expression levels in female than in male cones. Twenty-six microRNA (miRNA) families and 74 targets, including 46 24-nt sRNAs with a 5’ A, which are specifically expressed in male cones or female cones and probably bind to AGO4, were identified.

Conclusions

The sRNA pathways have higher activity in female than in male cones, and the miRNA pathways are the main sRNA pathways in P. tabuliformis. The low level of 24-nt short-interfering RNAs in conifers is not caused by the absence of biogenesis-related genes or AGO-binding proteins, but most likely caused by the low accumulation of these key components. The identification of sRNAs and their targets, as well as genes associated with sRNA biogenesis and action, will provide a good starting point for investigations into the roles of sRNA pathways in cone development in conifers.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1885-6) contains supplementary material, which is available to authorized users.

Combined Small RNA and Degradome Sequencing Reveals Novel MiRNAs and Their Targets in the High-Yield Mutant Wheat Strain Yunong 3114

Wheat is one of the main food sources worldwide; large amount studies have been conducted to improve wheat production. MicroRNAs (miRNAs) with about 20-30 nucleotide are a class of regulatory small RNAs (sRNAs), which could regulate gene expression through sequence-specific base pairing with target mRNAs, playing important roles in plant growth. An ideal plant architecture (IPA) is crucial to enhance yield in bread wheat. In this study, the high-yield wheat strain Yunong 3114 was EMS-mutagenesis from the wild-type strain Yunong 201, exhibiting a preferable plant structure compared with the wild-type strain. We constructed small RNA and degradome libraries from Yunong 201 and Yunong 3114, and performed small RNA sequencing of these libraries in order identify miRNAs and their targets related to IPA in wheat. Totally, we identified 488 known and 837 novel miRNAs from Yunong 3114 and 391 known and 533 novel miRNAs from Yunong 201. The number of miRNAs in the mutant increased. A total of 37 known and 432 putative novel miRNAs were specifically expressed in the mutant strain; furthermore, 23 known and 159 putative novel miRNAs were specifically expressed in the wild-type strain. A total of 150 known and 100 novel miRNAs were differentially expressed between mutant and wild-type strains. Among these differentially expressed novel miRNAs, 4 and 8 predict novel miRNAs were evidenced by degradome sequencing and showed up-regulated and down-regulated expressions in the mutant strain Yunong 3114, respectively. Targeted gene annotation and previous results indicated that this set of miRNAs is related to plant structure. Our results further suggested that miRNAs may be necessary to obtain an optimal wheat structure.

Ethylene resistance in flowering ornamental plants - improvements and future perspectives

Various strategies of plant breeding have been attempted in order to improve the ethylene resistance of flowering ornamental plants. These approaches span from conventional techniques such as simple cross-pollination to new breeding techniques which modify the plants genetically such as precise genome-editing. The main strategies target the ethylene pathway directly; others focus on changing the ethylene pathway indirectly via pathways that are known to be antagonistic to the ethylene pathway, e.g. increasing cytokinin levels. Many of the known elements of the ethylene pathway have been addressed experimentally with the aim of modulating the overall response of the plant to ethylene. Elements of the ethylene pathway that appear particularly promising in this respect include ethylene receptors as ETR1, and transcription factors such as EIN3. Both direct and indirect approaches seem to be successful, nevertheless, although genetic transformation using recombinant DNA has the ability to save much time in the breeding process, they are not readily used by breeders yet. This is primarily due to legislative issues, economic issues, difficulties of implementing this technology in some ornamental plants, as well as how these techniques are publically perceived, particularly in Europe. Recently, newer and more precise genome-editing techniques have become available and they are already being implemented in some crops. New breeding techniques may help change the current situation and pave the way toward a legal and public acceptance if products of these technologies are indistinguishable from plants obtained by conventional techniques.

Characterization and comparative profiling of the small RNA transcriptomes in two phases of flowering in Cymbidium ensifolium

Background

Cymbidium ensifolium is one of the most important ornamental flowers in China, with an elegant shape, beautiful appearance, and a fragrant aroma. Its unique flower shape has long attracted scientists. MicroRNAs (miRNAs) are critical regulators in plant development and physiology, including floral development. However, to date, few studies have examined miRNAs in C. ensifolium.

Results

In this study, we employed Solexa technology to sequence four small RNA libraries from two flowering phases to identify miRNAs related to floral development. We identified 48 mature conserved miRNA and 71 precursors. These conserved miRNA belonged to 20 families. We also identified 45 novel miRNA which includes 21 putative novel miRNAs*, and 28 hairpin forming precursors. Two trans-acting small interfering RNAs (ta-siRNAs) were identified, one of which was homologous to TAS3a1. TAS3a1 belongs to the TAS3 family, which has been previously reported to target auxin response factors (ARF) and be involved in plant growth and floral development. Moreover, we built a C. ensifolium transctriptome database to identify genes targeted by miRNA, which resulted in 790 transcriptomic target unigenes. The target unigenes were annotated with information from the non-redundant (Nr), gene ontology database (GO), eukaryotic orthologous groups (KOGs) and Kyoto encyclopedia of genes and genomes (KEGG) database. The unigenes included MADS-box transcription factors targeted by miR156, miR172 and miR5179, and various hormone responding factors targeted by miR159. The MADS-box transcription factors are well known to determine the identity of flower organs and hormone responding factors involved in floral development. In expression analysis, three novel and four conserved miRNA were differentially expressed between two phases of flowering. The results were confirmed by RNA-seq and qRT-PCR. The differential expression of two miRNA, miR160 and miR396, targeted ARFs and growth regulating factor (GRF), respectively. However, most of these small RNA were clustered in the uncharacterized group, which suggests there may be many novel small non-coding RNAs yet to be discovered.

Conclusion

Our study provides a diverse set of miRNAs related to cymbidium floral development and serves as a useful resource for investigating miRNA-mediated regulatory mechanisms of floral development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1764-1) contains supplementary material, which is available to authorized users.