MicroRNAS and their regulatory roles in plants

Mining NGS transcriptomes for miRNAs and dissecting their role in regulating growth, development, and secondary metabolites production in different organs of a medicinal herb, Picrorhiza kurroa

This study is the first endeavor on mining of miRNAs and analyzing their involvement in development and secondary metabolism of an endangered medicinal herb Picrorhiza kurroa (P. kurroa ). miRNAs are ubiquitous non-coding RNA species that target complementary sequences of mRNA and result in either translational repression or target degradation in eukaryotes. The role of miRNAs has not been investigated in P. kurroa which is a medicinal herb of industrial value due to the presence of secondary metabolites, picroside-I and picroside-II. Computational identification of miRNAs was done in 6 transcriptomes of P. kurroa generated from root, shoot, and stolon organs varying for growth, development, and culture conditions. All available plant miRNA entries were retrieved from miRBase and used as backend datasets to computationally identify conserved miRNAs in transcriptome data sets. Total 18 conserved miRNAs were detected in P. kurroa followed by target prediction and functional annotation which suggested their possible role in controlling various biological processes. Validation of miRNA and expression analysis by qRT-PCR and 5' RACE revealed that miRNA-4995 has a regulatory role in terpenoid biosynthesis ultimately affecting the production of picroside-I. miR-5532 and miR-5368 had negligible expression in field-grown samples as compared to in vitro-cultured samples suggesting their role in regulating P. kurroa growth in culture conditions. The study has thus identified novel functions for existing miRNAs which can be further validated for their potential regulatory role.

Transcriptome-Wide Identification of miRNAs and Their Targets from Typha angustifolia by RNA-Seq and Their Response to Cadmium Stress

MicroRNAs (miRNAs) play important roles in plant responses to environmental stress. In this work, we used high-throughput sequencing to analyze transcriptome and small RNAs (sRNAs) in Typha angustifolia under cadmium (Cd) stress. 57,608,230 raw reads were obtained from deep sequencing of a pooled cDNA library. Sequence assembly and analysis yielded 102,473 unigenes. We subsequently sequenced two sRNA libraries from T. angustifolia with or without Cd exposure respectively. Based on transcriptome data of T. angustifolia, we catalogued and analyzed the sRNAs, resulting in the identification of 114 conserved miRNAs and 41 novel candidate miRNAs in both small RNA libraries. In silico analysis revealed 764 targets for 89 conserved miRNAs and 21 novel miRNAs. Statistical analysis on sequencing reads abundance and experimental validation revealed that 4 conserved and 6 novel miRNAs showed specific expression. Combined with function of target genes, these results suggested that miRNAs might play a role in plant Cd stress response. This study provided the first transcriptome-based analysis of miRNAs and their targets responsive to Cd stress in T. angustifolia, which provide a framework for further analysis of miRNAs and their role in regulating plant responses to Cd stress.

An integrative approach to identify hexaploid wheat miRNAome associated with development and tolerance to abiotic stress

Background

Wheat is a major staple crop with broad adaptability to a wide range of environmental conditions. This adaptability involves several stress and developmentally responsive genes, in which microRNAs (miRNAs) have emerged as important regulatory factors. However, the currently used approaches to identify miRNAs in this polyploid complex system focus on conserved and highly expressed miRNAs avoiding regularly those that are often lineage-specific, condition-specific, or appeared recently in evolution. In addition, many environmental and biological factors affecting miRNA expression were not yet considered, resulting still in an incomplete repertoire of wheat miRNAs.

Results

We developed a conservation-independent technique based on an integrative approach that combines machine learning, bioinformatic tools, biological insights of known miRNA expression profiles and universal criteria of plant miRNAs to identify miRNAs with more confidence. The developed pipeline can potentially identify novel wheat miRNAs that share features common to several species or that are species specific or clade specific. It allowed the discovery of 199 miRNA candidates associated with different abiotic stresses and development stages. We also highlight from the raw data 267 miRNAs conserved with 43 miRBase families. The predicted miRNAs are highly associated with abiotic stress responses, tolerance and development. GO enrichment analysis showed that they may play biological and physiological roles associated with cold, salt and aluminum (Al) through auxin signaling pathways, regulation of gene expression, ubiquitination, transport, carbohydrates, gibberellins, lipid, glutathione and secondary metabolism, photosynthesis, as well as floral transition and flowering.

Conclusion

This approach provides a broad repertoire of hexaploid wheat miRNAs associated with abiotic stress responses, tolerance and development. These valuable resources of expressed wheat miRNAs will help in elucidating the regulatory mechanisms involved in freezing and Al responses and tolerance mechanisms as well as for development and flowering. In the long term, it may help in breeding stress tolerant plants.

Electronic supplementary material

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

PlantMirnaT: miRNA and mRNA integrated analysis fully utilizing characteristics of plant sequencing data

miRNA is known to regulate up to several hundreds coding genes, thus the integrated analysis of miRNA and mRNA expression data is an important problem. Unfortunately, the integrated analysis is challenging since it needs to consider expression data of two different types, miRNA and mRNA, and target relationship between miRNA and mRNA is not clear, especially when microarray data is used. Fortunately, due to the low sequencing cost, small RNA and RNA sequencing are routinely processed and we may be able to infer regulation relationships between miRNAs and mRNAs more accurately by using sequencing data. However, no method is developed specifically for sequencing data. Thus we developed PlantMirnaT, a new miRNA-mRNA integrated analysis system. To fully leverage the power of sequencing data, three major features are developed and implemented in PlantMirnaT. First, we implemented a plant-specific short read mapping tool based on recent discoveries on miRNA target relationship in plant. Second, we designed and implemented an algorithm considering miRNA targets in the full intragenic region, not just 3' UTR. Lastly but most importantly, our algorithm is designed to consider quantity of miRNA expression and its distribution on target mRNAs. The new algorithm was used to characterize rice under drought condition using our proprietary data. Our algorithm successfully discovered that two miRNAs, miRNA1425-5p, miRNA 398b, that are involved in suppression of glucose pathway in a naturally drought resistant rice, Vandana. The system can be downloaded at https://sites.google.com/site/biohealthinformaticslab/resources.

Implications of miR166 and miR159 induction to the basal response mechanisms of an andigena potato (Solanum tuberosum subsp. andigena) to salinity stress, predicted from network models in Arabidopsis

MicroRNA (miRNA) mediated changes in gene expression by post-transcriptional modulation of major regulatory transcription factors is a potent mechanism for integrating growth and stress-related responses. Exotic plants including many traditional varieties of Andean potatoes (Solanum tuberosum subsp. andigena) are known for better adaptation to marginal environments. Stress physiological studies confirmed earlier reports on the salinity tolerance potentials of certain andigena cultivars. Guided by the hypothesis that certain miRNAs play important roles in growth modulation under suboptimal conditions, we identified and characterized salinity stress-responsive miRNA-target gene pairs in the andigena cultivar Sullu by parallel analysis of noncoding and coding RNA transcriptomes. Inverse relationships were established by the reverse co-expression between two salinity stress-regulated miRNAs (miR166, miR159) and their target transcriptional regulators HD-ZIP-Phabulosa/Phavulota and Myb101, respectively. Based on heterologous models in Arabidopsis, the miR166-HD-ZIP-Phabulosa/Phavulota network appears to be involved in modulating growth perhaps by mediating vegetative dormancy, with linkages to defense-related pathways. The miR159-Myb101 network may be important for the modulation of vegetative growth while also controlling stress-induced premature transition to reproductive phase. We postulate that the induction of miR166 and miR159 under salinity stress represents important network hubs for balancing gene expression required for basal growth adjustments.

Small RNAs in plants: recent development and application for crop improvement

The phenomenon of RNA interference (RNAi) which involves sequence-specific gene regulation by small non-coding RNAs, i.e., small interfering RNA (siRNA) and microRNA (miRNA) has emerged as one of most powerful approaches for crop improvement. RNAi based on siRNA is one of the widely used tools of reverse genetics which aid in revealing gene functions in many species. This technology has been extensively applied to alter the gene expression in plants with an aim to achieve desirable traits. RNAi has been used for enhancing the crop yield and productivity by manipulating the gene involved in biomass, grain yield and enhanced shelf life of fruits and vegetables. It has also been applied for developing resistance against various biotic (bacteria, fungi, viruses, nematodes, insects) and abiotic stresses (drought, salinity, cold, etc.). Nutritional improvements of crops have also been achieved by enriching the crops with essential amino acids, fatty acids, antioxidants and other nutrients beneficial for human health or by reducing allergens or anti-nutrients. microRNAs are key regulators of important plant processes like growth, development, and response to various stresses. In spite of similarity in size (20-24 nt), miRNA differ from siRNA in precursor structures, pathway of biogenesis, and modes of action. This review also highlights the miRNA based genetic modification technology where various miRNAs/artificial miRNAs and their targets can be utilized for improving several desirable plant traits. microRNA based strategies are much efficient than siRNA-based RNAi strategies due to its specificity and less undesirable off target effects. As per the FDA guidelines, small RNA (sRNA) based transgenics are much safer for consumption than those over-expressing proteins. This review thereby summarizes the emerging advances and achievement in the field of sRNAs and its application for crop improvement.

Functional plasticity of miR165/166 in plant development revealed by small tandem target mimic

MicroRNA 165 and 166 (miR165/166) is composed of nine members and targets five members (PHB, PHV, REV, ATHB8 and ATHB15) of the HD-ZIP III transcription factor family. Mutants generated by traditional methods could hardly reveal the overall functions of miR165/166 in plant development. In this study, the expressions of all miR165/166 members were simultaneously blocked by over-expressing STTM165/166-31 in Arabidopsis and tomato for functional dissection of miR165/166 family. Following a down-regulation of over 90% endogenous miR165/166, the target HD-ZIP III genes were correspondingly up-regulated in the STTM transgenic Arabidopsis and tomato plants. Notably, the STTM165/166-31 over-expressed Arabidopsis and tomato displayed pleiotropic effects on development which were not frequently observed in previously identified genetic mutants of either individual miR165/166 gene or any of the five target genes. Furthermore, the transgenic Arabidopsis showed increased IAA content and decreased IAA sensitivity accompanied by enhanced expressions of genes responsible for auxin biosynthesis and signaling, suggesting possible roles of auxin in mediation of miR165/166-regulated processes. Importantly, the transgenic Arabidopsis exhibited the improved behavior under salt stress. Overall, such diverse variations in plant development and physiological process revealed by STTM165/166 demonstrate a key role of miR165/166-mediated network in regulating plant development and responses to abiotic stresses.

Application of microRNA gene resources in the improvement of agronomic traits in rice

microRNAs (miRNAs) are important nonprotein-coding genes that are involved in almost all biological processes, including cell differentiation and fate determination, developmental regulation, and immune responses. Investigations have shown that some miRNAs can highly affect plant agricultural traits, including virus resistance, nematode resistance, drought and salinity tolerance, heavy metal detoxification, biomass yield, grain yield, fruit development and flower development. Therefore, these miRNAs are considered a newly identified gene resource for the genetic improvement of crops. In this review, we will summarize the recent findings of the rice miRNA-directed regulatory network, which controls agronomic traits such as yield, quality and stress tolerance, and explore the outlook for the uses of these miRNA-associated traits in rice biotechnology.

Role of microRNAs in plant drought tolerance

Drought is a normal and recurring climate feature in most parts of the world and plays a major role in limiting crop productivity. However, plants have their own defence systems to cope with adverse climatic conditions. One of these defence mechanisms is the reprogramming of gene expression by microRNAs (miRNAs). miRNAs are small noncoding RNAs of approximately 22 nucleotides length, which have emerged as important regulators of genes at post-transcriptional levels in a range of organisms. Some miRNAs are functionally conserved across plant species and are regulated by drought stress. These properties suggest that miRNA-based genetic modifications have the potential to enhance drought tolerance in cereal crops. This review summarizes the current understanding of the regulatory mechanisms of plant miRNAs, involvement of plant miRNAs in drought stress responses in barley (Hordeum vulgare L.), wheat (Triticum spp.) and other plant species, and the involvement of miRNAs in plant-adaptive mechanisms under drought stress. Potential strategies and directions for future miRNA research and the utilization of miRNAs in the improvement of cereal crops for drought tolerance are also discussed.

Floral induction and flower formation--the role and potential applications of miRNAs

The multiple regulatory pathways controlling flowering and flower development are varied and complex, and they require tight control of gene expression and protein levels. MicroRNAs (miRNAs) act at both the transcriptional and post-transcriptional level to regulate key genes involved in flowering-related processes such as the juvenile-adult transition, the induction of floral competence and flower development. Many different miRNA families are involved in these processes and their roles are summarized in this review, along with potential biotechnological applications for miRNAs in controlling processes related to flowering and flower development.

Methylation of microRNA genes regulates gene expression in bisexual flower development in andromonoecious poplar

Previous studies showed sex-specific DNA methylation and expression of candidate genes in bisexual flowers of andromonoecious poplar, but the regulatory relationship between methylation and microRNAs (miRNAs) remains unclear. To investigate whether the methylation of miRNA genes regulates gene expression in bisexual flower development, the methylome, microRNA, and transcriptome were examined in female and male flowers of andromonoecious poplar. 27 636 methylated coding genes and 113 methylated miRNA genes were identified. In the coding genes, 64.5% of the methylated reads mapped to the gene body region; by contrast, 60.7% of methylated reads in miRNA genes mainly mapped in the 5' and 3' flanking regions. CHH methylation showed the highest methylation levels and CHG showed the lowest methylation levels. Correlation analysis showed a significant, negative, strand-specific correlation of methylation and miRNA gene expression (r=0.79, P <0.05). The methylated miRNA genes included eight long miRNAs (lmiRNAs) of 24 nucleotides and 11 miRNAs related to flower development. miRNA172b might play an important role in the regulation of bisexual flower development-related gene expression in andromonoecious poplar, via modification of methylation. Gynomonoecious, female, and male poplars were used to validate the methylation patterns of the miRNA172b gene, implying that hyper-methylation in andromonoecious and gynomonoecious poplar might function as an important regulator in bisexual flower development. Our data provide a useful resource for the study of flower development in poplar and improve our understanding of the effect of epigenetic regulation on genes other than protein-coding genes.

MicroRNA profiling analysis throughout tomato fruit development and ripening reveals potential regulatory role of RIN on microRNAs accumulation

The development and ripening of tomato fruit are complex processes involving many gene regulatory pathways at the transcriptional and post-transcriptional level. Ripening inhibitor (RIN) is a vital transcription factor, which targets numerous ripening-related genes at the transcriptional level during tomato fruit ripening. MicroRNAs (miRNAs) are a class of short noncoding RNAs that play important roles in post-transcriptional gene regulation. To elucidate the potential regulatory relationship between rin and miRNAs during fruit development and ripening, we identified known miRNAs and profiled their expression in wild-type tomato and rin mutant using a deep sequencing approach combined with quantitative RT-PCR. A total of 33 known miRNA families were identified, of which 14 miRNA families were differently accumulated. Subsequent promoter analysis showed that possible RIN-binding motifs (CArG-box) tended to occur frequently in the promoter regions of partial differently expressed miRNAs. In addition, ethylene may participate in the regulation of miRNAs accumulation during tomato fruit ripening. Chromatin immunoprecipitation (ChIP) and electrophoretic mobility shift assay confirmed the direct binding of RIN to the promoter of MIR172a. Collectively, these results showed a close correlation between miRNA expression and RIN as well as ethylene, which further elucidated the regulatory roles of miRNAs during fruit development and ripening and enriched the regulatory network of RIN in tomato fruit.

Intronic regions of plant genes potentially encode RDR (RNA-dependent RNA polymerase)-dependent small RNAs

Recent research has linked the non-coding intronic regions of plant genes to the production of small RNAs (sRNAs). Certain introns, called 'mirtrons' and 'sirtrons', could serve as the single-stranded RNA precursors for the generation of microRNA and small interfering RNA, respectively. However, whether the intronic regions could serve as the template for double-stranded RNA synthesis and then for sRNA biogenesis through an RDR (RNA-dependent RNA polymerase)-dependent pathway remains unclear. In this study, a genome-wide search was made for the RDR-dependent sRNA loci within the intronic regions of the Arabidopsis genes. Hundreds of intronic regions encoding three or more RDR-dependent sRNAs were found to be covered by dsRNA-seq (double-stranded RNA sequencing) reads, indicating that the intron-derived sRNAs were indeed generated from long double-stranded RNA precursors. More interestingly, phase-distributed sRNAs were discovered on some of the dsRNA-seq read-covered intronic regions, and those sRNAs were largely 24 nt in length. Based on these results, the opinion is put forward that the intronic regions might serve as the genomic origins for the RDR-dependent sRNAs. This opinion might add a novel layer to the current biogenesis model of the intron-derived sRNAs.

Identification and functional analysis of microRNAs and their targets in Platanus acerifolia under lead (Pb) stress

MicroRNAs (miRNAs) play important regulatory roles in development and stress responses in plants. Lead (Pb) is a non-essential element that is highly toxic to living organisms. Platanus acerifolia is grown as a street tree in cities throughout temperate regions for its importance in improving the urban ecological environment. MiRNAs that respond to abiotic stresses have been identified in plants; however, until now, the influence of Pb stress on P. acerifolia miRNAs has not been reported. To identify miRNAs and predict their target genes under Pb stress, two small RNA and two degradome libraries were constructed from Pb-treated and Pb-free leaves of P.acerifolia seedlings. After sequencing, 55 known miRNAs and 129 novel miRNAs were obtained, and 104 target genes for the miRNAs were identified by degradome sequencing. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed to predict the functions of the targets. The expressions of eight differentially expressed miRNAs were validated by quantitative real-time polymerase chain reaction (qRT-PCR). This is the first report about P. acerifolia miRNAs and their target genes under Pb stress. This study has provided data for further research into molecular mechanisms involved in resistance of P.acerifolia to Pb stress.

Decoding mechanisms by which silent codon changes influence protein biogenesis and function

SCOPE

Synonymous codon usage has been a focus of investigation since the discovery of the genetic code and its redundancy. The occurrences of synonymous codons vary between species and within genes of the same genome, known as codon usage bias. Today, bioinformatics and experimental data allow us to compose a global view of the mechanisms by which the redundancy of the genetic code contributes to the complexity of biological systems from affecting survival in prokaryotes, to fine tuning the structure and function of proteins in higher eukaryotes. Studies analyzing the consequences of synonymous codon changes in different organisms have revealed that they impact nucleic acid stability, protein levels, structure and function without altering amino acid sequence. As such, synonymous mutations inevitably contribute to the pathogenesis of complex human diseases. Yet, fundamental questions remain unresolved regarding the impact of silent mutations in human disorders. In the present review we describe developments in this area concentrating on mechanisms by which synonymous mutations may affect protein function and human health.

PURPOSE

This synopsis illustrates the significance of synonymous mutations in disease pathogenesis. We review the different steps of gene expression affected by silent mutations, and assess the benefits and possible harmful effects of codon optimization applied in the development of therapeutic biologics.

PHYSIOLOGICAL AND MEDICAL RELEVANCE

Understanding mechanisms by which synonymous mutations contribute to complex diseases such as cancer, neurodegeneration and genetic disorders, including the limitations of codon-optimized biologics, provides insight concerning interpretation of silent variants and future molecular therapies.

Identification of reference genes for quantitative expression analysis of microRNAs and mRNAs in barley under various stress conditions

For accurate and reliable gene expression analysis using quantitative real-time reverse transcription PCR (qPCR), the selection of appropriate reference genes as an internal control for normalization is crucial. We hypothesized that non-coding, small nucleolar RNAs (snoRNAs) would be stably expressed in different barley varieties and under different experimental treatments, in different tissues and at different developmental stages of plant growth and therefore might prove to be suitable reference genes for expression analysis of both microRNAs (miRNAs) and mRNAs. In this study, we examined the expression stability of ten candidate reference genes in six barley genotypes under five experimental stresses, drought, fungal infection, boron toxicity, nutrient deficiency and salinity. We compared four commonly used housekeeping genes; Actin (ACT), alpha-Tubulin (α-TUB), Glycolytic glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ADP-ribosylation factor 1-like protein (ADP), four snoRNAs; (U18, U61, snoR14 and snoR23) and two microRNAs (miR168, miR159) as candidate reference genes. We found that ADP, snoR14 and snoR23 were ranked as the best of these candidates across diverse samples. For accurate and reliable gene expression analysis using quantitative real-time reverse transcription PCR (qPCR), the selection of appropriate reference genes as an internal control for normalization is crucial. We hypothesized that non-coding, small nucleolar RNAs (snoRNAs) would be stably expressed in different barley varieties and under different experimental treatments, in different tissues and at different developmental stages of plant growth and therefore might prove to be suitable reference genes for expression analysis of both microRNAs (miRNAs) and mRNAs. In this study, we examined the expression stability of ten candidate reference genes in six barley genotypes under five experimental stresses, drought, fungal infection, boron toxicity, nutrient deficiency and salinity. We compared four commonly used housekeeping genes; Actin (ACT), alpha-Tubulin (α-TUB), Glycolytic glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ADP-ribosylation factor 1-like protein (ADP), four snoRNAs; (U18, U61, snoR14 and snoR23) and two microRNAs (miR168, miR159) as candidate reference genes. We found that ADP, snoR14 and snoR23 were ranked as the best of these candidates across diverse samples. Additionally, we found that miR168 was a suitable reference gene for expression analysis in barley. Finally, we validated the performance of our stable and unstable candidate reference genes for both mRNA and miRNA qPCR data normalization under different stress conditions and demonstrated the superiority of the stable candidates. Our data demonstrate the suitability of barley snoRNAs and miRNAs as potential reference genes for miRNA and mRNA qPCR data normalization under different stress treatments.

The role of plant small RNAs in NB-LRR regulation

Small RNAs constitute a fundamental layer of gene regulation for diverse biological processes in plants, including development, metabolism and stress responses. With the advance of high-throughput sequencing technologies and the rapid accumulation of transcriptomic data, the scope of regulation afforded by small RNAs has expanded to encompass plant innate immune responses. Plants have evolved the capacity to control the infection through intracellular surveillance proteins of the nucleotide binding site-leucine-rich repeat (NB-LRR) family that recognize pathogen-encoded effectors and initiate effector-triggered immunity. Emerging evidence indicates that plants have evolved to use specific microRNAs that target conserved domains of NB-LRR-encoding genes and trigger the production of a phased array of 21-nucleotide secondary small interfering RNAs to amplify the silencing effect. Herein, this review describes recent advances in understanding the roles of small RNAs in NB-LRR regulation that provide new insights into small RNA-mediated arms race between plants and their pathogens and discuss the unresolved questions and the future prospects for research on this topic.

High-Throughput Sequencing Reveals Diverse Sets of Conserved, Nonconserved, and Species-Specific miRNAs in Jute

MicroRNAs play a pivotal role in regulating a broad range of biological processes, acting by cleaving mRNAs or by translational repression. A group of plant microRNAs are evolutionarily conserved; however, others are expressed in a species-specific manner. Jute is an agroeconomically important fibre crop; nonetheless, no practical information is available for microRNAs in jute to date. In this study, Illumina sequencing revealed a total of 227 known microRNAs and 17 potential novel microRNA candidates in jute, of which 164 belong to 23 conserved families and the remaining 63 belong to 58 nonconserved families. Among a total of 81 identified microRNA families, 116 potential target genes were predicted for 39 families and 11 targets were predicted for 4 among the 17 identified novel microRNAs. For understanding better the functions of microRNAs, target genes were analyzed by Gene Ontology and their pathways illustrated by KEGG pathway analyses. The presence of microRNAs identified in jute was validated by stem-loop RT-PCR followed by end point PCR and qPCR for randomly selected 20 known and novel microRNAs. This study exhaustively identifies microRNAs and their target genes in jute which will ultimately pave the way for understanding their role in this crop and other crops.

Targeted genome modifications in soybean with CRISPR/Cas9

BACKGROUND

The ability to selectively alter genomic DNA sequences in vivo is a powerful tool for basic and applied research. The CRISPR/Cas9 system precisely mutates DNA sequences in a number of organisms. Here, the CRISPR/Cas9 system is shown to be effective in soybean by knocking-out a green fluorescent protein (GFP) transgene and modifying nine endogenous loci.

RESULTS

Targeted DNA mutations were detected in 95% of 88 hairy-root transgenic events analyzed. Bi-allelic mutations were detected in events transformed with eight of the nine targeting vectors. Small deletions were the most common type of mutation produced, although SNPs and short insertions were also observed. Homoeologous genes were successfully targeted singly and together, demonstrating that CRISPR/Cas9 can both selectively, and generally, target members of gene families. Somatic embryo cultures were also modified to enable the production of plants with heritable mutations, with the frequency of DNA modifications increasing with culture time. A novel cloning strategy and vector system based on In-Fusion® cloning was developed to simplify the production of CRISPR/Cas9 targeting vectors, which should be applicable for targeting any gene in any organism.

CONCLUSIONS

The CRISPR/Cas9 is a simple, efficient, and highly specific genome editing tool in soybean. Although some vectors are more efficient than others, it is possible to edit duplicated genes relatively easily. The vectors and methods developed here will be useful for the application of CRISPR/Cas9 to soybean and other plant species.

Genome-wide identification and analysis of drought-responsive genes and microRNAs in tobacco

Drought stress response is a complex trait regulated at transcriptional and post-transcriptional levels in tobacco. Since the 1990s, many studies have shown that miRNAs act in many ways to regulate target expression in plant growth, development and stress response. The recent draft genome sequence of Nicotiana benthamiana has provided a framework for Digital Gene Expression (DGE) and small RNA sequencing to understand patterns of transcription in the context of plant response to environmental stress. We sequenced and analyzed three Digital Gene Expression (DGE) libraries from roots of normal and drought-stressed tobacco plants, and four small RNA populations from roots, stems and leaves of control or drought-treated tobacco plants, respectively. We identified 276 candidate drought responsive genes (DRGs) with sequence similarities to 64 known DRGs from other model plant crops, 82 were transcription factors (TFs) including WRKY, NAC, ERF and bZIP families. Of these tobacco DRGs, 54 differentially expressed DRGs included 21 TFs, which belonged to 4 TF families such as NAC (6), MYB (4), ERF (10), and bZIP (1). Additionally, we confirmed expression of 39 known miRNA families (122 members) and five conserved miRNA families, which showed differential regulation under drought stress. Targets of miRNAs were further surveyed based on a recently published study, of which ten targets were DRGs. An integrated gene regulatory network is proposed for the molecular mechanisms of tobacco root response to drought stress using differentially expressed DRGs, the changed expression profiles of miRNAs and their target transcripts. This network analysis serves as a reference for future studies on tobacco response stresses such as drought, cold and heavy metals.

To bloom or not to bloom: role of microRNAs in plant flowering

During the course of their life cycles, plants undergo various morphological and physiological changes underlying juvenile-to-adult and adult-to-flowering phase transitions. To flower or not to flower is a key step of plasticity of a plant toward the start of its new life cycle. In addition to the previously revealed intrinsic genetic programs, exogenous cues, and endogenous cues, a class of small non-coding RNAs, microRNAs (miRNAs), plays a key role in plants making the decision to flower by integrating into the known flowering pathways. This review highlights the age-dependent flowering pathway with a focus on a number of timing miRNAs in determining such a key process. The contributions of other miRNAs which exist mainly outside the age pathway are also discussed. Approaches to study the flowering-determining miRNAs, their interactions, and applications are presented.

Identification and characterization of microRNAs from tree peony (Paeonia ostii) and their response to copper stress

MicroRNAs (miRNAs) are a class of non-coding, small RNAs recognized as important regulators of gene expression. Although plant miRNAs have been extensively studied in model systems, less is known in other plants with limited genome sequence data, including Paeonia ostii. In this work, we used high-throughput sequencing to identify conserved and nonconserved miRNAs and other short RNAs in Paeonia ostii under control and copper stressed condition. 102 previously known plant miRNAs were identified and classified into 89 families according to their gene sequence identity. Some miRNAs were highly conserved in the plant kingdom suggesting that these miRNA play important and conserved roles. Combined our transcriptome sequencing data of Paeonia ostii under same conditions, 34 novel potential miRNAs were identified. The potential targets of the identified known and novel miRNAs were also predicted based on sequence homology search. Comparing the two libraries, it was observed that 12 conserved miRNAs and 18 novel miRNAs showed significantly changes in response to copper stress. Some of the new identified potential miRNAs might be involved in Paeonia ostii-specific regulating mechanisms under copper stress. These results provide a framework for further analysis of miRNAs and their role in regulating Paeonia ostii response to copper stress.

Duplicate gene divergence by changes in microRNA binding sites in Arabidopsis and Brassica

Gene duplication provides large numbers of new genes that can lead to the evolution of new functions. Duplicated genes can diverge by changes in sequences, expression patterns, and functions. MicroRNAs play an important role in the regulation of gene expression in many eukaryotes. After duplication, two paralogs may diverge in their microRNA binding sites, which might impact their expression and function. Little is known about conservation and divergence of microRNA binding sites in duplicated genes in plants. We analyzed microRNA binding sites in duplicated genes in Arabidopsis thaliana and Brassica rapa. We found that duplicates are more often targeted by microRNAs than singletons. The vast majority of duplicated genes in A. thaliana with microRNA binding sites show divergence in those sites between paralogs. Analysis of microRNA binding sites in genes derived from the ancient whole-genome triplication in B. rapa also revealed extensive divergence. Paralog pairs with divergent microRNA binding sites show more divergence in expression patterns compared with paralog pairs with the same microRNA binding sites in Arabidopsis. Close to half of the cases of binding site divergence are caused by microRNAs that are specific to the Arabidopsis genus, indicating evolutionarily recent gain of binding sites after target gene duplication. We also show rapid evolution of microRNA binding sites in a jacalin gene family. Our analyses reveal a dynamic process of changes in microRNA binding sites after gene duplication in Arabidopsis and highlight the role of microRNA regulation in the divergence and contrasting evolutionary fates of duplicated genes.

Dissection of miRNA pathways using arabidopsis mesophyll protoplasts

MicroRNAs (miRNAs) control gene expression mostly post-transcriptionally by guiding transcript cleavage and/or translational repression of complementary mRNA targets, thereby regulating developmental processes and stress responses. Despite the remarkable expansion of the field, the mechanisms underlying miRNA activity are not fully understood. In this article, we describe a transient expression system in Arabidopsis mesophyll protoplasts, which is highly amenable for the dissection of miRNA pathways. We show that by transiently overexpressing primary miRNAs and target mimics, we can manipulate miRNA levels and consequently impact on their targets. Furthermore, we developed a set of luciferase-based sensors for quantifying miRNA activity that respond specifically to both endogenous and overexpressed miRNAs and target mimics. We demonstrate that these miRNA sensors can be used to test the impact of putative components of the miRNA pathway on miRNA activity, as well as the impact of specific mutations, by either overexpression or the use of protoplasts from the corresponding mutants. We further show that our miRNA sensors can be used for investigating the effect of chemicals on miRNA activity. Our cell-based transient expression system is fast and easy to set up, and generates quantitative results, being a powerful tool for assaying miRNA activity in vivo.

Pathogenesis of Soybean mosaic virus in soybean carrying Rsv1 gene is associated with miRNA and siRNA pathways, and breakdown of AGO1 homeostasis

Profiling small RNAs in soybean Williams 82 (rsv), susceptible to Soybean mosaic virus (SMV, the genus Potyvirus, family Potyviridae) strains G2 and G7, and soybean PI96983 (Rsv1), resistant to G2 but susceptible to G7, identified the microRNA miR168 that was highly overexpressed only in G7-infected PI96983 showing a lethal systemic hypersensitive response (LSHR). Overexpression of miR168 was in parallel with the high-level expression of AGO1 mRNA, high-level accumulation of miR168-mediated AGO1 mRNA cleavage products but with severely repressed AGO1 protein. In contrast, AGO1 mRNA, degradation products and protein remained without significant changes in G2- and G7-infected Williams 82. Moreover, knock-down of SGS3, an essential component in RNA silencing, suppressed AGO1 siRNA, partially recovered repressed AGO1 protein, and alleviated LSHR severity in G7-infected Rsv1 soybean. These results suggest that both miRNA and siRNA pathways are involved in G7 infection of Rsv1 soybean, and LSHR is associated with breakdown of AGO1 homeostasis.

Involvement of miR160/miR393 and their targets in cassava responses to anthracnose disease

Cassava is a starchy root crop for food and industrial applications in many countries around the world. Among the factors that affect cassava production, diseases remain the major cause of yield loss. Cassava anthracnose disease is caused by the fungus Colletotrichum gloeosporioides. Severe anthracnose attacks can cause tip die-backs and stem cankers, which can affect the availability of planting materials especially in large-scale production systems. Recent studies indicate that plants over- or under-express certain microRNAs (miRNAs) to cope with various stresses. Understanding how a disease-resistant plant protects itself from pathogens should help to uncover the role of miRNAs in the plant immune system. In this study, the disease severity assay revealed different response to C. gloeosporioides infection in two cassava cultivars. Quantitative RT-PCR analysis uncovered the differential expression of the two miRNAs and their target genes in the two cassava cultivars that were subjected to fungal infection. The more resistant cultivar revealed the up-regulation of miR160 and miR393, and consequently led to low transcript levels in their targets, ARF10 and TIR1, respectively. The more susceptible cultivar exhibited the opposite pattern. The cis-regulatory elements relevant to defense and stress responsiveness, fungal elicitor responsiveness and hormonal responses were the most prevalent present in the miRNAs gene promoter regions. The possible dual role of these specific miRNAs and their target genes associated with cassava responses to C. gloeosporioides is discussed. This is the first study to address the molecular events by which miRNAs which might play a role in fungal-infected cassava. A better understanding of the functions of miRNAs target genes should greatly increase our knowledge of the mechanism underlying susceptibility and lead to new strategies to enhance disease tolerance in this economically important crop.

miRNA expression during prickly pear cactus fruit development

miRNAs are a class of small non-coding RNAs that regulate gene expression. They are involved in the control of many developmental processes, including fruit development. The increasing amount of information on miRNAs, on their expression, abundance, and conservation between various species, provides a new opportunity to study the role of miRNAs in non-model plant species. In this work, we used a combination of Northern blot and tissue print hybridization analysis to identify conserved miRNAs expressed during prickly pear cactus (Opuntia ficus indica) fruit development. Comparative profiling detected the expression of 34 miRNAs, which were clustered in three different groups that were associated with the different phases of fruit development. Variation in the level of miRNA expression was observed. Gradual expression increase of several miRNAs was observed during fruit development, including miR164. miR164 was selected for stem-loop RT-PCR and for a detailed spatial-temporal expression analysis. At early floral stages, miR164 was mainly localized in meristematic tissues, boundaries and fusion zones, while it was more homogenously expressed in fruit tissues. Our results provide the first evidence of miRNA expression in the prickly pear cactus and provide the basis for future research on miRNAs in Opuntia. Moreover, our analyses suggest that miR164 plays different roles during prickly pear cactus fruit development.

HITS-CLIP and PAR-CLIP advance viral miRNA targetome analysis

MiRNAs regulate gene expression by binding predominantly to the 3' untranslated region (UTR) of target transcripts to prevent their translation and/or induce target degradation. In addition to the more than 1200 human miRNAs, human DNA tumor viruses such as Kaposi's sarcoma-associated herpesvirus (KSHV) and Epstein-Barr virus (EBV) encode miRNAs. Target predictions indicate that each miRNA targets hundreds of transcripts, many of which are regulated by multiple miRNAs. Thus, target identification is a big challenge for the field. Most methods used currently investigate single miRNA-target interactions and are not able to analyze complex miRNA-target networks. To overcome these challenges, cross-linking and immunoprecipitation (CLIP), a recently developed method to study direct RNA-protein interactions in living cells, has been successfully applied to miRNA target analysis. It utilizes Argonaute (Ago)-immunoprecipitation to isolate native Ago-miRNA-mRNA complexes. In four recent publications, two variants of the CLIP method (HITS-CLIP and PAR-CLIP) were utilized to determine the targetomes of human and viral miRNAs in cells infected with the gamma-herpesviruses KSHV and EBV, which are associated with a number of human cancers. Here, we briefly introduce herpesvirus-encoded miRNAs and then focus on how CLIP technology has largely impacted our understanding of viral miRNAs in viral biology and pathogenesis.

Tuning LeSPL-CNR expression by SlymiR157 affects tomato fruit ripening

In plants, microRNAs (miRNAs) play essential roles in growth, development, yield, stress response and interactions with pathogens. However no miRNA has been experimentally documented to be functionally involved in fruit ripening although many miRNAs have been profiled in fruits. Here we show that SlymiR157 and SlymiR156 differentially modulate ripening and softening in tomato (Solanum lycopersicum). SlymiR157 is expressed and developmentally regulated in normal tomato fruits and in those of the Colourless non-ripening (Cnr) epimutant. It regulates expression of the key ripening gene LeSPL-CNR in a likely dose-dependent manner through miRNA-induced mRNA degradation and translation repression. Viral delivery of either pre-SlymiR157 or mature SlymiR157 results in delayed ripening. Furthermore, qRT-PCR profiling of key ripening regulatory genes indicates that the SlymiR157-target LeSPL-CNR may affect expression of LeMADS-RIN, LeHB1, SlAP2a and SlTAGL1. However SlymiR156 does not affect the onset of ripening, but it impacts fruit softening after the red ripe stage. Our findings reveal that working together with a ripening network of transcription factors, SlymiR157 and SlymiR156 form a critical additional layer of regulatory control over the fruit ripening process in tomato.

Characterization of miRNAs associated with Botrytis cinerea infection of tomato leaves

BACKGROUND

Botrytis cinerea Pers. Fr. is an important pathogen causing stem rot in tomatoes grown indoors for extended periods. MicroRNAs (miRNAs) have been reported as gene expression regulators related to several stress responses and B. cinerea infection in tomato. However, the function of miRNAs in the resistance to B. cinerea remains unclear.

RESULTS

The miRNA expression patterns in tomato in response to B. cinerea stress were investigated by high-throughput sequencing. In total, 143 known miRNAs and seven novel miRNAs were identified and their corresponding expression was detected in mock- and B. cinerea-inoculated leaves. Among those, one novel and 57 known miRNAs were differentially expressed in B. cinerea-infected leaves, and 8 of these were further confirmed by quantitative reverse-transcription PCR (qRT-PCR). Moreover, five of these eight differentially expressed miRNAs could hit 10 coding sequences (CDSs) via CleaveLand pipeline and psRNAtarget program. In addition, qRT-PCR revealed that four targets were negatively correlated with their corresponding miRNAs (miR319, miR394, and miRn1).

CONCLUSION

Results of sRNA high-throughput sequencing revealed that the upregulation of miRNAs may be implicated in the mechanism by which tomato respond to B. cinerea stress. Analysis of the expression profiles of B. cinerea-responsive miRNAs and their targets strongly suggested that miR319, miR394, and miRn1 may be involved in the tomato leaves' response to B. cinerea infection.

Identification of microRNAs and target genes in the fruit and shoot tip of Lycium chinense: a traditional Chinese medicinal plant

Although Lycium chinense (goji berry) is an important traditional Chinese medicinal plant, little genome information is available for this plant, particularly at the small-RNA level. Recent findings indicate that the evolutionary role of miRNAs is very important for a better understanding of gene regulation in different plant species. To elucidate small RNAs and their potential target genes in fruit and shoot tissues, high-throughput RNA sequencing technology was used followed by qRT-PCR and RLM 5’-RACE experiments. A total of 60 conserved miRNAs belonging to 31 families and 30 putative novel miRNAs were identified. A total of 62 significantly differentially expressed miRNAs were identified, of which 15 (14 known and 1 novel) were shoot-specific, and 12 (7 known and 5 novel) were fruit-specific. Additionally, 28 differentially expressed miRNAs were recorded as up-regulated in fruit tissues. The predicted potential targets were involved in a wide range of metabolic and regulatory pathways. GO (Gene Ontology) enrichment analysis and the KEGG (Kyoto Encyclopedia of Genes and Genomes) database revealed that “metabolic pathways” is the most significant pathway with respect to the rich factor and gene numbers. Moreover, five miRNAs were related to fruit maturation, lycopene biosynthesis and signaling pathways, which might be important for the further study of fruit molecular biology. This study is the first, to detect known and novel miRNAs, and their potential targets, of L. chinense. The data and findings that are presented here might be a good source for the functional genomic study of medicinal plants and for understanding the links among diversified biological pathways.

Genome-Wide Characterization of miRNAs Involved in N Gene-Mediated Immunity in Response to Tobacco Mosaic Virus in Nicotiana benthamiana

microRNAs (miRNAs) are a class of endogenous small RNAs (sRNAs) that play pivotal roles in plant development, abiotic stress response, and pathogen response. miRNAs have been extensively studied in plants, but rarely in Nicotiana benthamiana, despite its wide use in plant virology studies, particularly for studying N protein-tobacco mosaic virus (TMV) interactions. We report an efficient method using high-throughput sequencing and bioinformatics to identify genome-wide miRNAs in N. benthamiana. A total of 30 conserved miRNA families and 113 novel miRNAs belonging to 93 families were identified. Some miRNAs were clustered on chromosomes, and some were embedded in host gene introns. The predicted miRNA targets were involved in diverse biological processes, such as metabolism, signaling, and responses to stimuli. miRNA expression profiling revealed that most of them were differentially expressed during N-mediated immunity to TMV. This study provides a framework for further analysis of miRNA functions in plant immunity.

Hsa-miR-1 downregulates long non-coding RNA urothelial cancer associated 1 in bladder cancer

MicroRNAs (miRNAs) are known to mainly target protein-coding genes at post-transcriptional level, resulting in mRNA destabilization and/or translational repression. Long non-coding RNAs (lncRNAs) are emerging as a novel set of targets for miRNAs. Here, we report that downregulated hsa-miR-1 and upregulated lncRNA urothelial cancer associated 1 (UCA1) were inversely expressed in bladder cancer. Hsa-miR-1 decreased the expression of UCA1 in bladder cancer cells in an Ago2-slicer-dependent manner. The binding site between UCA1 and hsa-miR-1 was confirmed. Overexpression of hsa-miR-1 inhibited bladder cancer cell growth, induced apoptosis, and decreased cell motility. Knockdown of UCA1 expression phenocopied the effects of upregulation of hsa-miR-1. Transfection of UCA1 expression vector partly reversed the changes caused by transfection of pre-miR-1 plasmids. This study provides evidence for hsa-miR-1 to play tumor suppressive roles via downregulating lncRNA UCA1 in bladder cancer, which may have potential therapeutic significance.

Identification of cold-inducible microRNAs in grapevine

Low temperature is one of the most important environmental factors that limits the geographical distribution and productivity of grapevine. However, the molecular mechanisms on how grapevine responds to cold stress remains to be elucidated. MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that play an essential role during plant development and stress responses. Although miRNAs and their targets have been identified in several Vitis species, their participation during cold accumulation in grapevine remains unknown. In this study, two small RNA libraries were generated from micropropagated 'Muscat Hamburg' (V. vinifera) plantlets under normal and low temperatures (4°C). A total of 163 known miRNAs and 67 putative novel miRNAs were detected from two small RNA libraries by Solexa sequencing. Forty-four cold-inducible miRNAs were identified through differentially expressed miRNAs (DEMs) analysis; among which, 13 belonged to upregulated DEMs while 31 belonged downregulated DEMs. The expression patterns of the 13 DEMs were verified by real-time RT-PCR analysis. The prediction of the target genes for DEMs indicated that miRNA may regulate transcription factors, including AP2, SBP, MYB, bHLH, GRAS, and bZIP under cold stress. The 5'-RLM RACE were conducted to verify the cleavage site of predicted targets. Seven predicted target genes for four known and three novel vvi-miRNAs showed specific cleavage sites corresponding to their miRNA complementary sequences. The expression pattern of these seven target genes revealed negative correlation with the expression level of the corresponding vvi-miRNAs. Our results indicated that a diverse set of miRNAs in V. vinifera are cold-inducible and may play an important role in cold stress response.

Genome-wide identification of microRNAs and their targets in wild type and phyB mutant provides a key link between microRNAs and the phyB-mediated light signaling pathway in rice

Phytochrome B (phyB), a member of the phytochrome family in rice, plays important roles in regulating a range of developmental processes and stress responses. However, little information about the mechanisms involved in the phyB-mediated light signaling pathway has been reported in rice. MicroRNAs (miRNAs) also perform important roles in plant development and stress responses. Thus, it is intriguing to explore the role of miRNAs in the phyB-mediated light signaling pathway in rice. In this study, comparative high-throughput sequencing and degradome analysis were used to identify candidate miRNAs and their targets that participate in the phyB-mediated light signaling pathway. A total of 720 known miRNAs, 704 novel miRNAs and 1957 target genes were identified from the fourth leaves of wild-type (WT) and phyB mutant rice at the five-leaf stage. Among them, 135 miRNAs showed differential expression, suggesting that the expression of these miRNAs is directly or indirectly under the control of phyB. In addition, 32 out of the 135 differentially expressed miRNAs were found to slice 70 genes in the rice genome. Analysis of these target genes showed that members of various transcription factor families constituted the largest proportion, indicating miRNAs are probably involved in the phyB-mediated light signaling pathway mainly by regulating the expression of transcription factors. Our results provide new clues for functional characterization of miRNAs in the phyB-mediated light signaling pathway, which should be helpful in comprehensively uncovering the molecular mechanisms of phytochrome-mediated photomorphogenesis and stress responses in plants.

Identification of miRNAs and their targets in transgenic Brassica napus and its acceptor (Westar) by high-throughput sequencing and degradome analysis

MicroRNAs (miRNAs) are a class of noncoding small RNAs (sRNAs) that play many roles in plant growth, development, and the stress response.

Differential expression of microRNAs and other small RNAs in barley between water and drought conditions

Drought is a major constraint to crop production, and microRNAs (miRNAs) play an important role in plant drought tolerance. Analysis of miRNAs and other classes of small RNAs (sRNAs) in barley grown under water and drought conditions reveals that drought selectively regulates expression of miRNAs and other classes of sRNAs. Low-expressed miRNAs and all repeat-associated siRNAs (rasiRNAs) tended towards down-regulation, while tRNA-derived sRNAs (tsRNAs) had the tendency to be up-regulated, under drought. Antisense sRNAs (putative siRNAs) did not have such a tendency under drought. In drought-tolerant transgenic barley overexpressing DREB transcription factor, most of the low-expressed miRNAs were also down-regulated. In contrast, tsRNAs, rasiRNAs and other classes of sRNAs were not consistently expressed between the drought-treated and transgenic plants. The differential expression of miRNAs and siRNAs was further confirmed by Northern hybridization and quantitative real-time PCR (qRT-PCR). Targets of the drought-regulated miRNAs and siRNAs were predicted, identified by degradome libraries and confirmed by qRT-PCR. Their functions are diverse, but most are involved in transcriptional regulation. Our data provide insight into the expression profiles of miRNAs and other sRNAs, and their relationship under drought, thereby helping understand how miRNAs and sRNAs respond to drought stress in cereal crops.

Contribution of bioinformatics prediction in microRNA-based cancer therapeutics

Despite enormous efforts, cancer remains one of the most lethal diseases in the world. With the advancement of high throughput technologies massive amounts of cancer data can be accessed and analyzed. Bioinformatics provides a platform to assist biologists in developing minimally invasive biomarkers to detect cancer, and in designing effective personalized therapies to treat cancer patients. Still, the early diagnosis, prognosis, and treatment of cancer are an open challenge for the research community. MicroRNAs (miRNAs) are small non-coding RNAs that serve to regulate gene expression. The discovery of deregulated miRNAs in cancer cells and tissues has led many to investigate the use of miRNAs as potential biomarkers for early detection, and as a therapeutic agent to treat cancer. Here we describe advancements in computational approaches to predict miRNAs and their targets, and discuss the role of bioinformatics in studying miRNAs in the context of human cancer.

A genome-wide identification of the miRNAome in response to salinity stress in date palm (Phoenix dactylifera L.)

Although date palm is relatively salt-tolerant, little is known about the underlying molecular mechanisms that contribute to its salt tolerance. Only recently, investigators have uncovered microRNA-mediated post-transcriptional gene regulation, which is critical for typical plant development and adaptation to stress conditions such as salinity. To identify conserved and novel miRNAs in date palm and to characterize miRNAs that could play a role in salt tolerance, we have generated sRNA libraries from the leaves and roots of NaCl-treated and untreated seedlings of date palm. Deep sequencing of these four sRNA libraries yielded approximately 251 million reads. The bioinformatics analysis has identified 153 homologs of conserved miRNAs, 89 miRNA variants, and 180 putative novel miRNAs in date palm. Expression profiles under salinity revealed differential regulation of some miRNAs in date palm. In leaves, 54 of the identified miRNAs were significantly affected and the majority (70%) of them were upregulated, whereas in roots, 25 of the identified miRNAs were significantly affected and 76% of them were upregulated by the salinity stress. The salt-responsiveness of some of these miRNAs was further validated using semi-quantitative PCR (qPCR). Some of the predicted targets for the identified miRNA include genes with known functions in plant salt tolerance, such as potassium channel AKT2-like proteins, vacuolar protein sorting-associated protein, calcium-dependent and mitogen-activated proteins. As one of the first cultivated trees in the world that can tolerate a wide range of abiotic stresses, date palm contains a large population of conserved and non-conserved miRNAs that function at the post-transcriptional level. This study provided insights into miRNA-mediated gene expression that are important for adaptation to salinity in date palms.

Targeting Non-Coding RNAs in Plants with the CRISPR-Cas Technology is a Challenge yet Worth Accepting

Non-coding RNAs (ncRNAs) have emerged as versatile master regulator of biological functions in recent years. MicroRNAs (miRNAs) are small endogenous ncRNAs of 18-24 nucleotides in length that originates from long self-complementary precursors. Besides their direct involvement in developmental processes, plant miRNAs play key roles in gene regulatory networks and varied biological processes. Alternatively, long ncRNAs (lncRNAs) are a large and diverse class of transcribed ncRNAs whose length exceed that of 200 nucleotides. Plant lncRNAs are transcribed by different RNA polymerases, showing diverse structural features. Plant lncRNAs also are important regulators of gene expression in diverse biological processes. There has been a breakthrough in the technology of genome editing, the CRISPR-Cas9 (clustered regulatory interspaced short palindromic repeats/CRISPR-associated protein 9) technology, in the last decade. CRISPR loci are transcribed into ncRNA and eventually form a functional complex with Cas9 and further guide the complex to cleave complementary invading DNA. The CRISPR-Cas technology has been successfully applied in model plants such as Arabidopsis and tobacco and important crops like wheat, maize, and rice. However, all these studies are focused on protein coding genes. Information about targeting non-coding genes is scarce. Hitherto, the CRISPR-Cas technology has been exclusively used in vertebrate systems to engineer miRNA/lncRNAs, but it is still relatively unexplored in plants. While briefing miRNAs, lncRNAs and applications of the CRISPR-Cas technology in human and animals, this review essentially elaborates several strategies to overcome the challenges of applying the CRISPR-Cas technology in editing ncRNAs in plants and the future perspective of this field.

Transcriptome analysis of a subtropical deciduous tree: autumn leaf senescence gene expression profile of Formosan gum

Autumn leaf senescence is a spectacular natural phenomenon; however, the regulation networks controlling autumnal colors and the leaf senescence program remain largely unelucidated. Whether regulation of leaf senescence is similar in subtropical deciduous plants and temperate deciduous plants is also unknown. In this study, the gene expression of a subtropical deciduous tree, Formosan gum (Liquidambar formosana Hance), was profiled. The transcriptomes of April leaves (green leaves, 'G') and December leaves (red leaves, 'R') were investigated by next-generation gene sequencing. Out of 58,402 de novo assembled contigs, 32,637 were annotated as putative genes. Furthermore, the L. formosana-specific microarray designed based on total contigs was used to extend the observation period throughout the growing seasons of 2011-2013. Network analysis from the gene expression profile focused on the genes up-regulated when autumn leaf senescence occurred. LfWRKY70, LfWRKY75, LfWRKY65, LfNAC1, LfSPL14, LfNAC100 and LfMYB113 were shown to be key regulators of leaf senescnece, and the genes regulated by LfWRKY75, LfNAC1 and LfMYB113 are candidates to link chlorophyll degradation and anthocyanin biosynthesis to senescence. In summary, the gene expression profiles over the entire year of the developing leaf from subtropical deciduous trees were used for in silico analysis and the putative gene regulation in autumn coloration and leaf senescence is discussed in this study.

Endogenous target mimics down-regulate miR160 mediation of ARF10, -16, and -17 cleavage during somatic embryogenesis in Dimocarpus longan Lour

MicroRNA160 plays a critical role in plant development by negatively regulating the auxin response factors ARF10, -16, and -17. However, the ways in which miR160 expression is regulated at the transcriptional level, and how miR160 interacts with its targets during plant embryo development, remain unknown. Here, we studied the regulatory relationships among endogenous target mimics (eTMs), and miR160 and its targets, and their involvement in hormone signaling and somatic embryogenesis (SE) in Dimocarpus longan. We identified miR160 family members and isolated the miR160 precursor, primary transcript, and promoter. The promoter contained cis-acting elements responsive to stimuli such as light, abscisic acid, salicylic acid (SA) and heat stress. The pri-miR160 was down-regulated in response to SA but up-regulated by gibberellic acid, ethylene, and methyl jasmonate treatment, suggesting that pri-miR160 was associated with hormone transduction. Dlo-miR160a, -a(∗) and -d(∗) reached expression peaks in torpedo-shaped embryos, globular embryos and cotyledonary embryos, respectively, but were barely detectable in friable-embryogenic callus. This suggests that they have expression-related and functional diversity, especially during the middle and later developmental stages of SE. Four potential eTMs for miR160 were identified. Two of them, glucan endo-1,3-beta- glucosidase-like protein 2-like and calpain-type cysteine protease DEK1, were confirmed to control the corresponding dlo-miR160a(∗) expression level. This suggests that they may function to abolish the binding between dlo-miR160a(∗) and its targets. These two eTMs also participated in 2,4-D and ABA signal transduction. DlARF10, -16, and -17 targeting by dlo-miR160a was confirmed; their expression levels were higher in friable-embryogenic callus and incomplete compact pro-embryogenic cultures and responded to 2,4-D, suggesting they may play a major role in the early stages of longan SE dependent on 2,4-D. The eTMs, miR160, and ARF10, -16, and -17 exhibited tissue specificity in 'Sijimi' longan vegetative and reproductive organs, but were not significant negatively correlated. These results provide insights into the possible role of the eTM-miR160-ARF10-16-17 pathway in longan somatic embryo development.

Endogenous target mimics down-regulate miR160 mediation of ARF10, -16, and -17 cleavage during somatic embryogenesis in Dimocarpus longan Lour

MicroRNA160 plays a critical role in plant development by negatively regulating the auxin response factors ARF10, -16, and -17. However, the ways in which miR160 expression is regulated at the transcriptional level, and how miR160 interacts with its targets during plant embryo development, remain unknown. Here, we studied the regulatory relationships among endogenous target mimics (eTMs), and miR160 and its targets, and their involvement in hormone signaling and somatic embryogenesis (SE) in Dimocarpus longan. We identified miR160 family members and isolated the miR160 precursor, primary transcript, and promoter. The promoter contained cis-acting elements responsive to stimuli such as light, abscisic acid, salicylic acid (SA) and heat stress. The pri-miR160 was down-regulated in response to SA but up-regulated by gibberellic acid, ethylene, and methyl jasmonate treatment, suggesting that pri-miR160 was associated with hormone transduction. Dlo-miR160a, -a(∗) and -d(∗) reached expression peaks in torpedo-shaped embryos, globular embryos and cotyledonary embryos, respectively, but were barely detectable in friable-embryogenic callus. This suggests that they have expression-related and functional diversity, especially during the middle and later developmental stages of SE. Four potential eTMs for miR160 were identified. Two of them, glucan endo-1,3-beta- glucosidase-like protein 2-like and calpain-type cysteine protease DEK1, were confirmed to control the corresponding dlo-miR160a(∗) expression level. This suggests that they may function to abolish the binding between dlo-miR160a(∗) and its targets. These two eTMs also participated in 2,4-D and ABA signal transduction. DlARF10, -16, and -17 targeting by dlo-miR160a was confirmed; their expression levels were higher in friable-embryogenic callus and incomplete compact pro-embryogenic cultures and responded to 2,4-D, suggesting they may play a major role in the early stages of longan SE dependent on 2,4-D. The eTMs, miR160, and ARF10, -16, and -17 exhibited tissue specificity in 'Sijimi' longan vegetative and reproductive organs, but were not significant negatively correlated. These results provide insights into the possible role of the eTM-miR160-ARF10-16-17 pathway in longan somatic embryo development.

Comparative expression profiling of miRNAs between the cytoplasmic male sterile line MeixiangA and its maintainer line MeixiangB during rice anther development

miRNAs are involved in the pollen development during the CMS occurrence in rice. miRNAs are 20-24 nt endogenously expressed small RNAs that play key roles in the regulation of many growth and developmental processes in plants. The knowledge on cytoplasmic male sterility (CMS) regulation by miRNAs in rice is rather limited. In this study, Illumina sequencing was employed to examine the expression profiles of rice anther miRNAs from the CMS line MeixiangA (MxA) and its maintainer line MeixiangB (MxB). A total of 518 known miRNAs and 144 novel miRNAs were identified during rice anther development. Based on the number of sequencing reads, a total of 24 miRNAs were discovered to be differentially expressed between MxA and MxB, and the results were partially validated by qRT-PCR. Among these, 16 miRNAs were decreased and 8 miRNAs were increased in MxA compared with MxB. Target prediction showed that they target genes encoding EF-hand family proteins, F-box domain-containing proteins, MYB transcription factors, PPR-containing proteins and transposons. The expression patterns for targets of osa-miR528, osa-miR5793, osa-miR1432, osa-miR159, osa-miR812d, osa-miR2118c, osa-miR172d and osa-miR5498 were selectively examined, and the results showed that there was a negative correlation on the expression patterns between miRNAs and their targets. These targets have previously been reported to be related with pollen development and male sterility, suggesting that miRNAs might act as regulators of CMS occurrence in rice anthers. Furthermore, miRNA editing events were observed. The U → C and U → A editing phenomenon was validated by molecular cloning and sequencing. These findings contribute to our understanding of the roles of miRNAs during anther development and CMS occurrence in rice.

Identification of chilling stress-responsive tomato microRNAs and their target genes by high-throughput sequencing and degradome analysis

BACKGROUND

MicroRNAs (miRNAs) are a class of noncoding small RNAs (sRNAs) that are 20-24 nucleotides (nt) in length. Extensive studies have indicated that miRNAs play versatile roles in plants, functioning in processes such as growth, development and stress responses. Chilling is a common abiotic stress that seriously affects plants growth and development. Recently, chilling-responsive miRNAs have been detected in several plant species. However, little is known about the miRNAs in the model plant tomato. 'LA1777' (Solanum habrochaites) has been shown to survive chilling stress due to its various characteristics.

RESULTS

Here, two small RNA libraries and two degradome libraries were produced from chilling-treated (CT) and non-chilling-treated (NT) leaves of S. habrochaites seedlings. Following high-throughput sequencing and filtering, 161 conserved and 236 novel miRNAs were identified in the two libraries. Of these miRNAs, 192 increased in the response to chilling stress while 205 decreased. Furthermore, the target genes of the miRNAs were predicted using a degradome sequencing approach. It was found that 62 target genes were cleaved by 42 conserved miRNAs, while nine target genes were cleaved by nine novel miRNAs. Additionally, nine miRNAs and six target genes were validated by quantitative real-time PCR (qRT-PCR). Target gene functional analysis showed that most target genes played positive roles in the chilling response, primarily by regulating the expression of anti-stress proteins, antioxidant enzyme and genes involved in cell wall formation.

CONCLUSIONS

Tomato is an important model plant for basic biological research. In this study, numerous conserved and novel miRNAs involved in the chilling response were identified using high-throughput sequencing, and the target genes were analyzed by degradome sequencing. The work helps identify chilling-responsive miRNAs in tomato and increases the number of identified miRNAs involved in chilling stress. Furthermore, the work provides a foundation for further study of the regulation of miRNAs in the plant response to chilling stress.

The miR9863 family regulates distinct Mla alleles in barley to attenuate NLR receptor-triggered disease resistance and cell-death signaling

Barley (Hordeum vulgare L.) Mla alleles encode coiled-coil (CC), nucleotide binding, leucine-rich repeat (NB-LRR) receptors that trigger isolate-specific immune responses against the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh). How Mla or NB-LRR genes in grass species are regulated at post-transcriptional level is not clear. The microRNA family, miR9863, comprises four members that differentially regulate distinct Mla alleles in barley. We show that miR9863 members guide the cleavage of Mla1 transcripts in barley, and block or reduce the accumulation of MLA1 protein in the heterologous Nicotiana benthamiana expression system. Regulation specificity is determined by variation in a unique single-nucleotide-polymorphism (SNP) in mature miR9863 family members and two SNPs in the Mla miR9863-binding site that separates these alleles into three groups. Further, we demonstrate that 22-nt miR9863s trigger the biogenesis of 21-nt phased siRNAs (phasiRNAs) and together these sRNAs form a feed-forward regulation network for repressing the expression of group I Mla alleles. Overexpression of miR9863 members specifically attenuates MLA1, but not MLA10-triggered disease resistance and cell-death signaling. We propose a key role of the miR9863 family in dampening immune response signaling triggered by a group of MLA immune receptors in barley.

Plants rely on cell-surface and intracellular immune receptors to sense pathogen invasion and to mediate defense responses. However, uncontrolled activation of immune responses is harmful to plant growth and development. Small RNAs have recently been shown to fine-tune the expression of intracellular immune receptors and contribute to the regulation of defense signaling in dicot plants, while similar processes have not been well documented in monocot grain crops, such as barley and wheat. Here, we show that, in barley, some members of the miR9863 family target a subset of Mla alleles that confer race-specific disease resistance to the powdery mildew fungus. These miRNAs act on Mla transcripts by cleavage and translational repression. Production of a type of trans-acting small RNAs, designated as phasiRNAs, enhances the effects of miRNA regulation on Mla targets. We propose that Mla-mediated immune signaling is fine-tuned by the miRNAs at later stage of MLA activation to avoid overloading of immune responses in barley cells.

Small RNA profiling reveals regulation of Arabidopsis miR168 and heterochromatic siRNA415 in response to fungal elicitors

BACKGROUND

Small RNAs (sRNAs), including small interfering RNAs (siRNAs) and microRNAs (miRNAs), have emerged as important regulators of eukaryotic gene expression. In plants, miRNAs play critical roles in development, nutrient homeostasis and abiotic stress responses. Accumulating evidence also reveals that sRNAs are involved in plant immunity. Most studies on pathogen-regulated sRNAs have been conducted in Arabidopsis plants infected with the bacterial pathogen Pseudomonas syringae, or treated with the flagelin-derived elicitor peptide flg22 from P. syringae. This work investigates sRNAs that are regulated by elicitors from the fungus Fusarium oxysporum in Arabidopsis.

RESULTS

Microarray analysis revealed alterations on the accumulation of a set of sRNAs in response to elicitor treatment, including miRNAs and small RNA sequences derived from massively parallel signature sequencing. Among the elicitor-regulated miRNAs was miR168 which regulates ARGONAUTE1, the core component of the RNA-induced silencing complex involved in miRNA functioning. Promoter analysis in transgenic Arabidopsis plants revealed transcriptional activation of MIR168 by fungal elicitors. Furthermore, transgenic plants expressing a GFP-miR168 sensor gene confirmed that the elicitor-induced miR168 is active. MiR823, targeting Chromomethylase3 (CMT3) involved in RNA-directed DNA methylation (RdDM) was also found to be regulated by fungal elicitors. In addition to known miRNAs, microarray analysis allowed the identification of an elicitor-inducible small RNA that was incorrectly annotated as a miRNA. Studies on Arabidopsis mutants impaired in small RNA biogenesis demonstrated that this sRNA, is a heterochromatic-siRNA (hc-siRNA) named as siRNA415. Hc-siRNAs are known to be involved in RNA-directed DNA methylation (RdDM). SiRNA415 is detected in several plant species.

CONCLUSION

Results here presented support a transcriptional regulatory mechanism underlying MIR168 expression. This finding highlights the importance of miRNA functioning in adaptive processes of Arabidopsis plants to fungal infection. The results of this study also lay a foundation for the involvement of RdDM processes through the activity of siRNA415 and miR823 in mediating regulation of immune responses in Arabidopsis plants.

Evidence of purifying selection and co-evolution at the fold-back arm of the novel precursor microRNA159 gene in Phalaenopsis Species (Orchidaceae)

BACKGROUND

MicroRNAs (miRNAs) are small, endogenously transcribed, non-protein-coding RNAs that play important roles in regulation of gene expression in animals and plants. Here, selective constraints on the novel precursor microRNA159 (pre-miR159) gene were investigated in 42 Phalaenopsis species (Orchidaceae).

METHODS/RESULTS

A novel precursor microRNA159 gene was isolated from 42 Phalaenopsis species using a new microRNA-PCR (miR-PCR) approach. Sequencing of pre-miR159 genes revealed differences from the canonical pre-miR159 gene in Phalaenopsis species and other plants. Results demonstrated that the 5' and 3' fold-back arms and the terminal loop of the novel pre-miR159 gene have undergone purifying selection and selective constraint for stabilizing the secondary hairpin structure. Two conserved motifs within the 5' fold-back arm had the highest purifying selective pressure within the novel pre-miR159 gene. Evidence of sequence co-evolution between the 5' and 3' fold-back regions was observed.

CONCLUSIONS

Functional selective pressure might arise from the constraint of forming a hairpin structure and demonstrate co-evolution of sequences between the 5' and 3' fold-back regions of the novel pre-miR159 gene in Phalaenopsis species.

Genome-wide identification of different dormant Medicago sativa L. MicroRNAs in response to fall dormancy

BACKGROUND

MicroRNAs (miRNAs) are a class of regulatory small RNAs (sRNAs) that regulate gene post-transcriptional expression in plants and animals. High-throughput sequencing technology is capable of identifying small RNAs in plant species. Alfalfa (Medicago sativa L.) is one of the most widely cultivated perennial forage legumes worldwide, and fall dormancy is an adaptive characteristic related to the biomass production and winter survival in alfalfa. Here, we applied high-throughput sRNA sequencing to identify some miRNAs that were responsive to fall dormancy in standard variety (Maverick and CUF101) of alfalfa.

RESULTS

Four sRNA libraries were generated and sequenced from alfalfa leaves in two typical varieties at distinct seasons. Through integrative analysis, we identified 51 novel miRNA candidates of 206 families. Additionally, we identified 28 miRNAs associated with fall dormancy in standard variety (Maverick and CUF101), including 20 known miRNAs and eight novel miRNAs. Both high-throughput sequencing and RT-qPCR confirmed that eight known miRNA members were up-regulated and six known miRNA members were down-regulated in response to fall dormancy in standard variety (Maverick and CUF101). Among the 51 novel miRNA candidates, five miRNAs were up-regulated and three miRNAs were down-regulated in response to fall dormancy in standard variety (Maverick and CUF101), and five of them were confirmed by Northern blot analysis.

CONCLUSION

We identified 20 known miRNAs and eight new miRNA candidates that were responsive to fall dormancy in standard variety (Maverick and CUF101) by high-throughput sequencing of small RNAs from Medicago sativa. Our data provide a useful resource for investigating miRNA-mediated regulatory mechanisms of fall dormancy in alfalfa, and these findings are important for our understanding of the roles played by miRNAs in the response of plants to abiotic stress in general and fall dormancy in alfalfa.