The functional scope of plant microRNA-mediated silencing

Review: Plant microRNAs in pathogen defense: A panacea or a piece of the puzzle?

Plant microRNAs (miRNAs) control key agronomic traits that are associated with their conserved role(s) in development. However, despite a multitude of studies, the utility of miRNAs in plant-pathogen resistance remains less certain. Reviewing the literature identifies three general classes of miRNAs regarding plant pathogen defense. Firstly, a number of evolutionary dynamic 22 nucleotide miRNA families that repress large numbers of plant immunity genes, either directly, or through triggering the biogenesis of secondary siRNAs. However, understanding of their role in defense and of their manipulation to enhance pathogen resistance are still lacking. Secondly, highly conserved miRNAs that indirectly impact disease resistance through their targets that are primarily regulating development or hormone signaling. Any alteration of these miRNAs usually results in pleiotropic impacts, which may alter disease resistance in some plant species, and against some pathogens. Thirdly, are the comparatively diverse and evolutionary dynamic set of non-conserved miRNAs, some of which contribute to pathogen resistance, but whose narrow evolutionary presence will likely restrict their utility. Therefore, reflecting the diverse and evolving nature of plant-pathogen interactions, a complex interplay of plant miRNAs with pathogen responses exists. Any miRNA-based solution for pathogen resistance will likely be highly specific, rather than a general panacea.

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 role of microRNAs in responses to drought and heat stress in peanut (Arachis hypogaea)

MicroRNAs (miRNAs) are 21-24 nt small RNAs (sRNAs) that negatively regulate protein-coding genes and/or trigger phased small-interfering RNA (phasiRNA) production. Two thousand nine hundred miRNA families, of which ∼40 are deeply conserved, have been identified in ∼80 different plant species genomes. miRNA functions in response to abiotic stresses is less understood than their roles in development. Only seven peanut MIRNA families are documented in miRBase, yet a reference genome assembly is now published and over 480 plant-like MIRNA loci were predicted in the diploid peanut progenitor Arachis duranensis genome. We explored by computational analysis of a leaf sRNA library and publicly available sRNA, degradome, and transcriptome datasets the miRNA and phasiRNA space associated with drought and heat stresses in peanut. We characterized 33 novel candidate and 33 ancient conserved families of MIRNAs and present degradome evidence for their cleavage activities on mRNA targets, including several noncanonical targets and novel phasiRNA-producing noncoding and mRNA loci with validated novel targets such as miR1509 targeting serine/threonine-protein phosphatase7 and miRc20 and ahy-miR3514 targeting penta-tricopeptide repeats (PPRs), in contradistinction to other claims of miR1509/173/7122 superfamily miRNAs indirectly targeting PPRs via TAS-like noncoding RNA loci. We characterized the inverse correlations of significantly differentially expressed drought- and heat-regulated miRNAs, assayed by sRNA blots or transcriptome datasets, with target mRNA expressions in the same datasets. Meta-analysis of an expression atlas and over representation of miRNA target genes in co-expression networks suggest that miRNAs have functions in unique aspects of peanut gynophore development. Genome-wide MIRNA annotation of the published allopolyploid peanut genome can facilitate molecular breeding of value-added traits.

Target Landscape of Conserved Plant MicroRNAs and the Complexities of Their Ancient MicroRNA-Binding Sites

In plants, microRNA (miRNA)-target interactions (MTIs) require high complementarity, a feature from which bioinformatic programs have predicted numerous and diverse targets for any given miRNA, promoting the idea of complex miRNA networks. Opposing this is a hypothesis of constrained miRNA specificity, in which functional MTIs are restricted to the few targets whose required expression output is compatible with the expression of the miRNA. To explore these opposing views, the bioinformatic pipeline Targets Ranked Using Experimental Evidence was applied to strongly conserved miRNAs to identity their high-evidence (HE) targets across species. For each miRNA family, HE targets predominantly consisted of homologs from one conserved target gene family (primary family). These primary families corresponded to the known canonical miRNA-target families, validating the approach. Very few additional HE target families were identified (secondary family), and if so, they were likely functionally related to the primary family. Many primary target families contained highly conserved nucleotide sequences flanking their miRNA-binding sites that were enriched in HE homologs across species. A number of these flanking sequences are predicted to form conserved RNA secondary structures that preferentially base pair with the miRNA-binding site, implying that these sites are highly structured. Our findings support a target landscape view that is dominated by the conserved primary target families, with a minority of either secondary target families or non-conserved targets. This is consistent with the constrained hypothesis of functional miRNA specificity, which potentially in part is being facilitated by features beyond complementarity.

TRUEE; a bioinformatic pipeline to define the functional microRNA targetome of Arabidopsis

Central to plant microRNA (miRNA) biology is the identification of functional miRNA-target interactions (MTIs). However, the complementarity basis of bioinformatic target prediction results in mostly false positives, and the degree of complementarity does not equate with regulation. Here, we develop a bioinformatic workflow named TRUEE (Targets Ranked Using Experimental Evidence) that ranks MTIs on the extent to which they are subjected to miRNA-mediated cleavage. It sorts predicted targets into high (HE) and low evidence (LE) groupings based on the frequency and strength of miRNA-guided cleavage degradome signals across multiple degradome experiments. From this, each target is assigned a numerical value, termed a Category Score, ranking the extent to which it is subjected to miRNA-mediated cleavage. As a proof-of-concept, the 428 Arabidopsis miRNAs annotated in miRBase were processed through the TRUEE pipeline to determine the miRNA 'targetome'. The majority of high-ranking Category Score targets corresponded to highly conserved MTIs, validating the workflow. Very few Arabidopsis-specific, Brassicaceae-specific, or Conserved-passenger miRNAs had HE targets with high Category Scores. In total, only several hundred MTIs were found to have Category Scores characteristic of currently known physiologically significance MTIs. Although non-exhaustive, clearly the number of functional MTIs is much narrower than many studies claim. Therefore, using TRUEE to numerically rank targets  directly on experimental evidence has given insights into the scope of the functional miRNA targetome of Arabidopsis.

microRNA172 targets APETALA2 to regulate flavonoid biosynthesis in apple (Malus domestica)

MicroRNA172 (miR172) plays a role in regulating a diverse range of plant developmental processes, including flowering, fruit development and nodulation. However, its role in regulating flavonoid biosynthesis is unclear. In this study, we show that transgenic apple plants over-expressing miR172 show a reduction in red coloration and anthocyanin accumulation in various tissue types. This reduction was consistent with decreased expression of APETALA2 homolog MdAP2_1a (a miR172 target gene), MdMYB10, and targets of MdMYB10, as demonstrated by both RNA-seq and qRT-PCR analyses. The positive role of MdAP2_1a in regulating anthocyanin biosynthesis was supported by the enhanced petal anthocyanin accumulation in transgenic tobacco plants overexpressing MdAP2_1a, and by the reduction in anthocyanin accumulation in apple and cherry fruits transfected with an MdAP2_1a virus-induced-gene-silencing construct. We demonstrated that MdAP2_1a could bind directly to the promoter and protein sequences of MdMYB10 in yeast and tobacco, and enhance MdMYB10 promotor activity. In Arabidopsis, over-expression of miR172 reduced flavonoid (including anthocyanins and flavonols) concentration and RNA transcript abundance of flavonoid genes in plantlets cultured on medium containing 7% sucrose. The anthocyanin content and RNA abundance of anthocyanin genes could be partially restored by using a synonymous mutant of MdAP2_1a, which had lost the miR172 target sequences at mRNA level, but not restored by using a WT MdAP2_1a. These results indicate that miR172 inhibits flavonoid biosynthesis through suppressing the expression of an AP2 transcription factor that positively regulates MdMYB10.

Laccase Directed Lignification Is One of the Major Processes Associated With the Defense Response Against Pythium ultimum Infection in Apple Roots

Apple replant disease (ARD), incited by a pathogen complex including Pythium ultimum, causes stunted growth or death of newly planted trees at replant sites. Development and deployment of resistant or tolerant rootstocks offers a cost-effective, ecologically friendly, and durable approach for ARD management. Maximized exploitation of natural resistance requires integrated efforts to identify key regulatory mechanisms underlying resistance traits in apple. In this study, miRNA profiling and degradome sequencing identified major miRNA pathways and candidate genes using six apple rootstock genotypes with contrasting phenotypes to P. ultimum infection. The comprehensive RNA-seq dataset offered an expansive view of post-transcriptional regulation of apple root defense activation in response to infection from P. ultimum. Several pairs of miRNA families and their corresponding targets were identified for their roles in defense response in apple roots, including miR397-laccase, miR398-superoxide dismutase, miR10986-polyphenol oxidase, miR482-resistance genes, and miR160-auxin response factor. Of these families, the genotype-specific expression patterns of miR397 indicated its fundamental role in developing defense response patterns to P. ultimum infection. Combined with other identified copper proteins, the importance of cellular fortification, such as lignification of root tissues by the action of laccase, may critically contribute to genotype-specific resistance traits. Our findings suggest that quick and enhanced lignification of apple roots may significantly impede pathogen penetration and minimize the disruption of effective defense activation in roots of resistant genotypes. The identified target miRNA species and target genes consist of a valuable resource for subsequent functional analysis of their roles during interaction between apple roots and P. ultimum.

Genome-wide analysis of plant miRNA action clarifies levels of regulatory dynamics across developmental contexts

Development of complex organisms requires the delicate and dynamic spatiotemporal regulation of gene expression. Central to this are microRNAs (miRNAs). These mobile small RNAs offer specificity in conveying positional information and versatility in patterning the outcomes of gene expression. However, the parameters that shape miRNA output during development are still to be clarified. Here, we address this question on a genome-wide scale, using the maize shoot apex as a model. We show that patterns and levels of miRNA accumulation are largely determined at the transcriptional level, but are finessed post-transcriptionally in a tissue-dependent manner. The stem cell environments of the shoot apical meristem and vasculature appear particularly liable to this. Tissue-specific effects are also apparent at the level of target repression, with target cleavage products in the vasculature exceeding those of other tissues. Our results argue against a clearance mode of regulation purely at the level of transcript cleavage, leading us to propose that transcript cleavage provides a baseline level of target repression, onto which miRNA-driven translational repression can act to toggle the mode of target regulation between clearance and rheostat. Our data show how the inherent complexities of miRNA pathways allow the accumulation and activity of these small RNAs to be tailored in space and time to bring about the gene expression versatility needed during development.

MYB Transcription Factors and Its Regulation in Secondary Cell Wall Formation and Lignin Biosynthesis during Xylem Development

The secondary wall is the main part of wood and is composed of cellulose, xylan, lignin, and small amounts of structural proteins and enzymes. Lignin molecules can interact directly or indirectly with cellulose, xylan and other polysaccharide molecules in the cell wall, increasing the mechanical strength and hydrophobicity of plant cells and tissues and facilitating the long-distance transportation of water in plants. MYBs (v-myb avian myeloblastosis viral oncogene homolog) belong to one of the largest superfamilies of transcription factors, the members of which regulate secondary cell-wall formation by promoting/inhibiting the biosynthesis of lignin, cellulose, and xylan. Among them, MYB46 and MYB83, which comprise the second layer of the main switch of secondary cell-wall biosynthesis, coordinate upstream and downstream secondary wall synthesis-related transcription factors. In addition, MYB transcription factors other than MYB46/83, as well as noncoding RNAs, hormones, and other factors, interact with one another to regulate the biosynthesis of the secondary wall. Here, we discuss the biosynthesis of secondary wall, classification and functions of MYB transcription factors and their regulation of lignin polymerization and secondary cell-wall formation during wood formation.

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.

Base-Pairing Requirements for Small RNA-Mediated Gene Silencing of Recessive Self-Incompatibility Alleles in Arabidopsis halleri

Small noncoding RNAs are central regulators of genome activity and stability. Their regulatory function typically involves sequence similarity with their target sites, but understanding the criteria by which they specifically recognize and regulate their targets across the genome remains a major challenge in the field, especially in the face of the diversity of silencing pathways involved. The dominance hierarchy among self-incompatibility alleles in Brassicaceae is controlled by interactions between a highly diversified set of small noncoding RNAs produced by dominant S-alleles and their corresponding target sites on recessive S-alleles. By controlled crosses, we created numerous heterozygous combinations of S-alleles in Arabidopsis halleri and developed an real-time quantitative PCR assay to compare allele-specific transcript levels for the pollen determinant of self-incompatibility (SCR). This provides the unique opportunity to evaluate the precise base-pairing requirements for effective transcriptional regulation of this target gene. We found strong transcriptional silencing of recessive SCR alleles in all heterozygote combinations examined. A simple threshold model of base pairing for the small RNA-target interaction captures most of the variation in SCR transcript levels. For a subset of S-alleles, we also measured allele-specific transcript levels of the determinant of pistil specificity (SRK), and found sharply distinct expression dynamics throughout flower development between SCR and SRK In contrast to SCR, both SRK alleles were expressed at similar levels in the heterozygote genotypes examined, suggesting no transcriptional control of dominance for this gene. We discuss the implications for the evolutionary processes associated with the origin and maintenance of the dominance hierarchy among self-incompatibility alleles.

An Early Arising Role of the MicroRNA156/529-SPL Module in Reproductive Development Revealed by the Liverwort Marchantia polymorpha

In angiosperms, the phase transition from vegetative to reproductive growth involves the de-repression of the squamosa promoter-binding-protein-like (SPL) class of transcription factors, which is negatively regulated by the specific microRNAs (miRNAs/miRs) miR156/529 [1]. Non-vascular land plants also undergo growth-phase transition to the reproductive state, but knowledge regarding the controlling mechanisms is limited. Here, we investigate the reproductive transition in the liverwort Marchantia polymorpha, focusing on the roles of miR529c [2-4] and MpSPL2. First, we established mir529c-null mutants using CRISPR/Cas9. Even in the absence of far-red light-supplemented long-day condition, which is usually needed to induce reproductive development [5, 6], the mutant thalli developed sexual reproductive organs (gametangia) and produced gametes. Transgenic plants expressing a miR529-resistant MpSPL2 transgene also showed a similar phenotype of reproductive transition in the absence of inductive far-red light signals. In these mutants and transgenic plants, the MpSPL2 mRNA abundance was elevated. Mpspl2^ko mutant plants showed successful gamete development and fertilization, which suggests that MpSPL2 is involved in, but not essential for, sexual reproduction in M. polymorpha. Furthermore, analysis of Mpspl2^ko mutant and its complemented lines suggests that MpSPL2 may have a role in promotion of reproductive transition. These findings support the notion that the transition to reproductive development in liverworts is controlled by a system similar to that in angiosperms, and the miR156/529-SPL module has common significance in the control of the vegetative-to-reproductive transition during development in many land plants, including liverworts.

Biology and Function of miR159 in Plants

MicroR159 (miR159) is ancient, being present in the majority of land plants where it targets a class of regulatory genes called GAMYB or GAMYB-like via highly conserved miR159-binding sites. These GAMYB genes encode R2R3 MYB domain transcription factors that transduce the gibberellin (GA) signal in the seed aleurone and the anther tapetum. Here, GAMYB plays a conserved role in promoting the programmed cell death of these tissues, where miR159 function appears weak. By contrast, GAMYB is not involved in GA-signaling in vegetative tissues, but rather its expression is deleterious, leading to the inhibition of growth and development. Here, the major function of miR159 is to mediate strong silencing of GAMYB to enable normal growth. Highlighting this requirement of strong silencing are conserved RNA secondary structures associated with the miR159-binding site in GAMYB mRNA that promotes miR159-mediated repression. Although the miR159-GAMYB pathway in vegetative tissues has been implicated in a number of different functions, presently no conserved role for this pathway has emerged. We will review the current knowledge of the different proposed functions of miR159, and how this ancient pathway has been used as a model to help form our understanding of miRNA biology in plants.

Expression of miR159 Is Altered in Tomato Plants Undergoing Drought Stress

In a scenario of global climate change, water scarcity is a major threat for agriculture, severely limiting crop yields. Therefore, alternatives are urgently needed for improving plant adaptation to drought stress. Among them, gene expression reprogramming by microRNAs (miRNAs) might offer a biotechnologically sound strategy. Drought-responsive miRNAs have been reported in many plant species, and some of them are known to participate in complex regulatory networks via their regulation of transcription factors involved in water stress signaling. We explored the role of miR159 in the response of Solanum lycopersicum Mill. plants to drought stress by analyzing the expression of sly-miR159 and its target SlMYB transcription factor genes in tomato plants of cv. Ailsa Craig grown in deprived water conditions or in response to mechanical damage caused by the Colorado potato beetle, a devastating insect pest of Solanaceae plants. Results showed that sly-miR159 regulatory function in the tomato plants response to distinct stresses might be mediated by differential stress-specific MYB transcription factor targeting. sly-miR159 targeting of SlMYB33 transcription factor transcript correlated with accumulation of the osmoprotective compounds proline and putrescine, which promote drought tolerance. This highlights the potential role of sly-miR159 in tomato plants’ adaptation to water deficit conditions.

Cloning and characterization of a tyrosine decarboxylase involved in the biosynthesis of galanthamine in Lycoris aurea

Background

Galanthamine, one kind of Amaryllidaceae alkaloid extracted from the Lycoris species, is used in the treatment of Alzheimer’s disease. In regards to medical and economic importance, the biosynthesis and regulatory mechanism of the secondary metabolites in Lycoris remain uninvestigated.

Methods

BLAST was used to identify the sequence of tyrosine decarboxylase in the transcriptome of Lycoris aurea (L’Hér) Herb. The enzyme activity of this TYDC was determined by using heterologous expressed protein in the Escherichia coli cells. The related productive contents of tyramine were detected using High Performance Liquid Chromatography (HPLC). According to the available micro RNA sequencing profiles and degradome database of L. aurea, microRNA396 were isolated, which targets to LaTYDC1 and RNA Ligase-Mediated-Rapid Amplification of cDNA Ends (RLM-RACE) were used to confirm the cleavage. The expression levels of miR396 and LaTYDC1 were measured using a quantitative real-time polymerase chain reaction (qRT-PCR).

Results

LaTYDC1 was mainly expressed in root, bulb, leaf and flower fitting the models for galanthamine accumulation. This decarboxylase efficiently catalyzes tyrosine to tyramine conversion. Under methyl jasmonate (MeJA) treatment, the expression of LaTYDC1 and the content of tyramine sharply increase. The use of RLM-RACE confirms that miR396 promotes the degradation of LaTYDC1 mRNA. Under MeJA treatment, the expression of miR396 was suppressed while the expression level of LaTYDC1 sharply increased. Following the increase of the miR396 transcriptional level, LaTYDC1 was significantly repressed.

Conclusion

LaTYDC1 participates in the biosynthesis of galanthamine, and is regulated by miR396. This finding also provides genetic strategy for improving the yield of galanthamine in the future.

A Resource for Inactivation of MicroRNAs Using Short Tandem Target Mimic Technology in Model and Crop Plants

microRNAs (miRNAs) are endogenous small non-coding RNAs that bind to mRNAs and target them for cleavage and/or translational repression, leading to gene silencing. We previously developed short tandem target mimic (STTM) technology to deactivate endogenous miRNAs in Arabidopsis. Here, we created hundreds of STTMs that target both conserved and species-specific miRNAs in Arabidopsis, tomato, rice, and maize, providing a resource for the functional interrogation of miRNAs. We not only revealed the functions of several miRNAs in plant development, but also demonstrated that tissue-specific inactivation of a few miRNAs in rice leads to an increase in grain size without adversely affecting overall plant growth and development. RNA-seq and small RNA-seq analyses of STTM156/157 and STTM165/166 transgenic plants revealed the roles of these miRNAs in plant hormone biosynthesis and activation, secondary metabolism, and ion-channel activity-associated electrophysiology, demonstrating that STTM technology is an effective approach for studying miRNA functions. To facilitate the study and application of STTM transgenic plants and to provide a useful platform for storing and sharing of information about miRNA-regulated gene networks, we have established an online Genome Browser (https://blossom.ffr.mtu.edu/designindex2.php) to display the transcriptomic and miRNAomic changes in STTM-induced miRNA knockdown plants.

MicroRNA MIMIC binding sites: Minor flanking nucleotide alterations can strongly impact MIMIC silencing efficacy in Arabidopsis

In plants, microRNA (miRNA) target MIMICs (MIMs) have been widely used to inhibit miRNA function. They are based on the Arabidopsis INSENSITIVE TO PHOSPATE STARVATION 1 (IPS1) gene that corresponds to a non-coding RNA containing a miR399 binding site that can be modified to sequester and inhibit any miRNA of interest. However, the efficacy of miRNA inhibition of these different MIMs can vary greatly. Using MIMs that have strong efficacy (MIM159) and poor efficacy (MIM165), we investigate the underlying cause of this variation. Firstly, sequence alignments of IPS1 homologs from the Brassicaceae identified a highly conserved sequence immediately downstream of the miRNA binding site. Mutating this sequence in the context of the MIM159 attenuates its strong efficacy. This conserved flanking region contains a predicted stem-loop structure that is also predicted to be present in most modified MIMs that appear to have a strong efficacy, but not in MIM165 that has a poor efficacy. Restoring this predicted stem-loop in MIM165 via mutation of only three or five nucleotides within the conserved flanking region resulted in MIM165 variants that have very strong efficacies of miRNA inhibition. However, specifically mutating this predicted stem-loop in the MIM159 context failed to significantly reduce efficacy, and additional mutations to restore this predicted stem-loop weakened efficacy further. Although this shows there is no simple correlation between this predicted stem-loop and efficacy, these results add to the growing evidence that the sequence context of miRNA binding sites is important, and that minor nucleotide substitutions to flanking sequences of miRNA binding sites can strongly enhance or attenuate the miRNA-target interaction.

Tamarix microRNA Profiling Reveals New Insight into Salt Tolerance

The halophyte tamarisk (Tamarix) is extremely salt tolerant, making it an ideal material for salt tolerance-related studies. Although many salt-responsive genes of Tamarix were identified in previous studies, there are no reports on the role of post-transcriptional regulation in its salt tolerance. We constructed six small RNA libraries of Tamarix chinensis roots with NaCl treatments. High-throughput sequencing of the six libraries was performed and microRNA expression profiles were constructed. We investigated salt-responsive microRNAs to uncover the microRNA-mediated genes regulation. From these analyses, 251 conserved and 18 novel microRNA were identified from all small RNAs. From 191 differentially expressed microRNAs, 74 co-expressed microRNAs were identified as salt-responsive candidate microRNAs. The most enriched GO (gene ontology) terms for the 157 genes targeted by differentially expressed microRNAs suggested that transcriptions factors were highly active. Two hub microRNAs (miR414, miR5658), which connected by several target genes into an organic microRNA regulatory network, appeared to be the key regulators of post-transcriptional salt-stress responses. As the first survey on the tamarisk small RNAome, this study improves the understanding of tamarisk salt-tolerance mechanisms and will contribute to the molecular-assisted resistance breeding.

Evolution of microRNA827 targeting in the plant kingdom

Unlike most ancient microRNAs, which conservatively target homologous genes across species, microRNA827 (miR827) targets two different types of SPX (SYG1/PHO81/XPR1)-domain-containing genes, NITROGEN LIMITATION ADAPTATION (NLA) and PHOSPHATE TRANSPORTER 5 (PHT5), in Arabidopsis thaliana and Oryza sativa to regulate phosphate (Pi) transport and storage, respectively. However, how miR827 shifted its target preference and its evolutionary history are unknown. Based on target prediction analysis, we found that in most angiosperms, miR827 conservatively targets PHT5 homologs, but in Brassicaceae and Cleomaceae it preferentially targets NLA homologs, and we provide evidence for the transition of target preference during Brassicales evolution. Intriguingly, we found a lineage-specific loss of the miR827-regulatory module in legumes. Analysis of miR827-mediated cleavage efficiency and the expression of PHT5 in A. thaliana indicated that accumulation of mutations in the target site and the exclusion of the target site by alternative transcriptional initiation eliminated PHT5 targeting by miR827. Here, we identified a transition of miR827 target preference during plant evolution and revealed the uniqueness of miR827-mediated regulation among conserved plant miRNAs. Despite the change in its target preference, upregulation of miR827 by Pi starvation and its role in regulating cellular Pi homeostasis were retained.

Effects of blue light on flavonoid accumulation linked to the expression of miR393, miR394 and miR395 in longan embryogenic calli

While flavonoid metabolism’s regulation under light conditions by structural genes and transcription factors is understood, the roles of microRNAs (miRNAs) in this pathway have been rarely reported. In this paper, the accurate control of light was firstly enabled through the specially designed plant growth chamber which ensures consistency and accuracy of the cultivation of longan ECs and the repeatability of the experiments. Then, longan ECs were cultured in this chamber for 25 days. The change of growth rate of longan ECs was compared under different light qualities (dark, blue, green, white, green), intensities (16, 32, 64, 128, 256 μmol ·m^-2 ·s^-1), and durations (8 h, 12 h, 16 h, 20h, 24h). Results indicated that longan ECs had a high growth rate in the condition of blue or green light, at intensity ranged from 16 μmol·m^-2·s^-1 to 64 μmol·m^-2·s^-1, and duration from 8 h to 16 h. In addition, the contents of total flavonoids, rutin, and epicatechin were determined. Results indicated that flavonoid contents of longan ECs reached the highest value under blue light, at 32 μmol·m^-2·s^-1 and 12h/d. Blue light promoted the accumulation of epicatechin, but inhibited the synthesis of rutin. Finally, the expressions of flavonoid pathway genes, miRNAs and target genes were analyzed by qPCR. These results indicated that miR393 and its target gene DlTIR1-3, miR394 and its target gene DlAlMT12, and miR395 and its target gene DlAPS1 had a negative regulating relationship under blue light in longan ECs. Furthermore, miR393, miR394, and miR395 acted on target genes, which negatively regulated flavonoid key genes DlFLS and positively regulated key genes DlCHS, DlCHI, DlF3′H, DlDFR, DlLAR, and finally affected the accumulation of flavonoids. The treatment of longan ECs under the blue light at the intensity of 32 μmol·m^-2·s^-1 for 12 h/d inhibited the expression of miR393, miR394 and miR395, which promoted the expression of target genes and the accumulation of flavonoids and epicatechin, but inhibited the synthesis of rutin.

In Silico Analysis of Small RNAs Suggest Roles for Novel and Conserved miRNAs in the Formation of Epigenetic Memory in Somatic Embryos of Norway Spruce

Epigenetic memory in Norway spruce affects the timing of bud burst and bud set, vitally important adaptive traits for this long-lived forest species. Epigenetic memory is established in response to the temperature conditions during embryogenesis. Somatic embryogenesis at different epitype inducing (EpI) temperatures closely mimics the natural processes of epigenetic memory formation in seeds, giving rise to epigenetically different clonal plants in a reproducible and predictable manner, with respect to altered bud phenology. MicroRNAs (miRNAs) and other small non-coding RNAs (sRNAs) play an essential role in the regulation of plant gene expression and may affect this epigenetic mechanism. We used NGS sequencing and computational in silico methods to identify and profile conserved and novel miRNAs among small RNAs in embryogenic tissues of Norway spruce at three EpI temperatures (18, 23 and 28°C). We detected three predominant classes of sRNAs related to a length of 24 nt, followed by a 21–22 nt class and a third 31 nt class of sRNAs. More than 2100 different miRNAs within the prevailing length 21–22 nt were identified. Profiling these putative miRNAs allowed identification of 1053 highly expressed miRNAs, including 523 conserved and 530 novels. 654 of these miRNAs were found to be differentially expressed (DEM) depending on EpI temperature. For most DEMs, we defined their putative mRNA targets. The targets represented mostly by transcripts of multiple-repeats proteins, like TIR, NBS-LRR, PPR and TPR repeat, Clathrin/VPS proteins, Myb-like, AP2, etc. Notably, 124 DE miRNAs targeted 203 differentially expressed epigenetic regulators. Developing Norway spruce embryos possess a more complex sRNA structure than that reported for somatic tissues. A variety of the predicted miRNAs showed distinct EpI temperature dependent expression patterns. These putative EpI miRNAs target spruce genes with a wide range of functions, including genes known to be involved in epigenetic regulation, which in turn could provide a feedback process leading to the formation of epigenetic marks. We suggest that TIR, NBS and LRR domain containing proteins could fulfill more general functions for signal transduction from external environmental stimuli and conversion them into molecular response. Fine-tuning of the miRNA production likely participates in both developmental regulation and epigenetic memory formation in Norway spruce.

New technologies accelerate the exploration of non-coding RNAs in horticultural plants

Non-coding RNAs (ncRNAs), that is, RNAs not translated into proteins, are crucial regulators of a variety of biological processes in plants. While protein-encoding genes have been relatively well-annotated in sequenced genomes, accounting for a small portion of the genome space in plants, the universe of plant ncRNAs is rapidly expanding. Recent advances in experimental and computational technologies have generated a great momentum for discovery and functional characterization of ncRNAs. Here we summarize the classification and known biological functions of plant ncRNAs, review the application of next-generation sequencing (NGS) technology and ribosome profiling technology to ncRNA discovery in horticultural plants and discuss the application of new technologies, especially the new genome-editing tool clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 (Cas9) systems, to functional characterization of plant ncRNAs.

Definition and identification of small RNA sponges: Focus on miRNA sequestration

Targeting RNAs appears as an important opportunity to modulate biological processes. Here, we overviewed critical parameters implied in RNAs competition to bind small RNAs. These competitions influence small RNA availability and thereby gene expression and cell fate. We focused on the ability of RNAs to sequester small RNA, mainly the microRNAs (miRNAs) and proposed experimental workflows to demonstrate the existence and activity of RNA-sponge. From this basic science, we detailed tailored oligonucleotides, developed to challenge the binding of small RNA. In vitro and in vivo, these tailored oligonucleotides efficiently restore small RNA activity by preventing their sequestration on RNA-sponges.

Identification and differential regulation of microRNAs in response to methyl jasmonate treatment in Lycoris aurea by deep sequencing

Background

Lycoris aurea is a medicine-valuable and ornamental herb widely distributed in China. Former studied have showed that methyl jasmonate (MJ) treatment could increase the content of glanthamine-a worldwide medicine for symptomatic treatment of Alzheimer’s disease in genus Lycoris plants. To explore the possible role of miRNAs in the regulation of jasmonic acid signaling pathway and uncover their potential correlations, we investigated the expression profiles of small RNAs (sRNAs) and their targets in Lycoris aurea, with MJ treatment by using next-generation deep sequencing.

Results

A total of 365 miRNAs were identified, comprising 342 known miRNAs (representing 60 miRNA families) and 23 novel miRNAs. Among them, 143 known and 11 novel miRNAs were expressed differently under MJ treatment. Quantitative real-time PCR of eight selected miRNAs validated the expression pattern of these loci in response to MJ treatment. In addition, degradome sequencing analysis showed that 32 target genes were validated to be targeted by the 49 miRNAs, respectively. Gene function and pathway analyses showed that these targets such as auxin response factors (ARFs), squamosa promoter-binding like (SPL) proteins, basic helix-loop-helix (bHLH) proteins, and ubiquitin-conjugating enzyme E2 are involved in different plant processes, indicating miRNAs mediated regulation might play important roles in L. aurea response to MJ treatment. Furthermore, several L. aurea miRNAs associated with their target genes that might be involved in Amaryllidaceae alkloids biosynthehsis were also analyzed.

Conclusions

A number of miRNAs with diverse expression patterns, and complex relationships between expression of miRNAs and targets were identified. This study represents the first transcriptome-based analysis of miRNAs in Lycoris and will contribute to understanding the potential roles of miRNAs involved in regulation of MJ response.

Electronic supplementary material

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

Integrated RNA-seq and sRNA-seq analysis reveals miRNA effects on secondary metabolism in Solanum tuberosum L

Light is a major environmental factor that affects metabolic pathways and stimulates the production of secondary metabolites in potato. However, adaptive changes in potato metabolic pathways and physiological functions triggered by light are partly explained by gene expression changes. Regulation of secondary metabolic pathways in potato has been extensively studied at transcriptional level, but little is known about the mechanisms of post-transcriptional regulation by miRNAs. To identify light-responsive miRNAs/mRNAs and construct putative metabolism pathways regulated by the miRNA-mRNA pairs, an integrated omics (sRNAome and transcriptome) analysis was performed to potato under light stimulus. A total of 31 and 48 miRNAs were identified to be differentially expressed in the leaves and tubers, respectively. Among the DEGs, 1353 genes in the leaves and 1841 genes in the tubers were upregulated, while 1595 genes in the leaves and 897 genes in the tubers were downregulated by light. Mapman enrichment analyses showed that genes related to MVA pathway, alkaloids-like, phenylpropanoids, flavonoids, and carotenoids metabolism were significantly upregulated, while genes associated with major CHO metabolism were repressed in the leaves and tubers. Integrated miRNA and mRNA profiles revealed that light-responsive miRNAs are important regulators in alkaloids metabolism, UMP salvage, lipid biosynthesis, and cellulose catabolism. Moreover, several miRNAs may participate in glycoalkaloids metabolism via JA signaling pathway, UDP-glucose biosynthesis and hydroxylation reaction. This study provides a global view of miRNA and mRNA expression profiles in potato response to light, our results suggest that miRNAs might play important roles in secondary metabolic pathways, especially in glycoalkaloid biosynthesis. The findings will enlighten us on the genetic regulation of secondary metabolite pathways and pave the way for future application of genetically engineered potato.

MicroRNA858 Is a Potential Regulator of Phenylpropanoid Pathway and Plant Development

MicroRNAs (miRNAs) are endogenous, noncoding small RNAs that function as critical regulators of gene expression. In plants, miRNAs have shown their potential as regulators of growth, development, signal transduction, and stress tolerance. Although the miRNA-mediated regulation of several processes is known, the involvement of miRNAs in regulating secondary plant product biosynthesis is poorly understood. In this study, we functionally characterized Arabidopsis (Arabidopsis thaliana) miR858a, which putatively targets R2R3-MYB transcription factors involved in flavonoid biosynthesis. Overexpression of miR858a in Arabidopsis led to the down-regulation of several MYB transcription factors regulating flavonoid biosynthesis. In contrast to the robust growth and early flowering of miR858OX plants, reduction of plant growth and delayed flowering were observed in Arabidopsis transgenic lines expressing an artificial miRNA target mimic (MIM858). Genome-wide expression analysis using transgenic lines suggested that miR858a targets a number of regulatory factors that modulate the expression of downstream genes involved in plant development and hormonal and stress responses. Furthermore, higher expression of MYBs in MIM858 lines leads to redirection of the metabolic flux towards the synthesis of flavonoids at the cost of lignin synthesis. Altogether, our study has established the potential role of light-regulated miR858a in flavonoid biosynthesis and plant growth and development.

Identification of microRNAs in rice root in response to nitrate and ammonium

Nitrate and ammonium are two major nitrogen (N) sources for higher plants, but they differ in utilization and signaling. MicroRNAs (miRNAs) play an essential role in N signal transduction; however, knowledge remains limited about the regulatory role of miRNAs responsive to different N sources, especially in crop plants. To get global overview on miRNAs involved in N response in rice, we performed high-throughput small RNA-sequencing under different nitrate and ammonium treatments. The results demonstrated that only 16 and 11 miRNAs were significantly induced by nitrate and ammonium under short-term treatment, respectively. However, 60 differentially expressed miRNAs were found between nitrate and ammonium under long-term cultivation. These results suggested that miRNA response greatly differentiates between nitrate and ammonium treatments. Furthermore, 44 miRNAs were found to be differentially expressed between high- and low-N conditions. Our study reveals comprehensive expression profiling of miRNAs responsive to different N sources and different N treatments, which advances our understanding on the regulation of different N signaling and homeostasis mediated by miRNAs.

The role of small RNAs in vegetative shoot development

Shoot development consists of the production of lateral organs in predictable spatial and temporal patterns at the shoot apex. To properly integrate such programs of growth across different cell and tissue types, plants require highly complex and robust genetic networks. Over the last twenty years, the roles of small, non-coding RNAs (sRNAs) in these networks have become increasingly apparent, not least in vegetative shoot growth. In this review, we describe recent progress in understanding the contribution of sRNAs to the regulation of vegetative shoot growth, and outline persisting experimental limitations in the field.

Identification of the Abundant Hydroxyproline-Rich Glycoproteins in the Root Walls of Wild-Type Arabidopsis, an ext3 Mutant Line, and Its Phenotypic Revertant

Extensins are members of the cell wall hydroxyproline-rich glycoprotein (HRGP) superfamily that form covalently cross-linked networks in primary cell walls. A knockout mutation in EXT3 (AT1G21310), the gene coding EXTENSIN 3 (EXT3) in Arabidopsis Landsberg erecta resulted in a lethal phenotype, although about 20% of the knockout plants have an apparently normal phenotype (ANP). In this study the root cell wall HRGP components of wild-type, ANP and the ext3 mutant seedlings were characterized by peptide fractionation of trypsin digested anhydrous hydrogen fluoride deglycosylated wall residues and by sequencing using LC-MS/MS. Several HRGPs, including EXT3, were identified in the wild-type root walls but not in walls of the ANP and lethal mutant. Indeed the ANP walls and walls of mutants displaying the lethal phenotype possessed HRGPs, but the profiles suggest that changes in the amount and perhaps type may account for the corresponding phenotypes.