The small RNA world of plants

miR156b-targeted VvSBP8/13 functions downstream of the abscisic acid signal to regulate anthocyanins biosynthesis in grapevine fruit under drought

Anthocyanins are the primary color components of grapevine berries and wines. In cultivation practices, a moderate water deficit can promote anthocyanin accumulation in red grape skins. Our previous study showed that abscisic acid (ABA) plays a key role in this process. Herein, we identified a microRNA, vv-miR156b, that is generated in grapevine berries in response to drought stress, along with increasing anthocyanin content and biosynthetic structural gene transcripts. In contrast, vv-miR156b short tandem target mimic (STTM) function-loss callus exhibits the opposite phenotype. Results from in vivo and in vitro experiments revealed that the ABA-signaling-regulated transcription factor VvAREB2 binds directly to the ABA-responsive element (ABRE) of the MIR156b promoter and activates miR156b expression. Furthermore, two miR156b downstream targets, VvSBP8 and VvSBP13, exhibited reduced grape anthocyanin content in their overexpressors but there was a contrary result in their CRISPR-edited lines, the decrease in anthocyanin content was rescued in miR156b and SBP8/13 double overexpressors. We further demonstrated that both VvSBP8 and VvSBP13, encoding transcriptional repressors, displayed sufficient ability to interact with VvMYC1 and VvMYBA1, thereby interfering with MYB-bHLH-WD (MBW) repeat transcriptional complex formation, resulting in the repression of anthocyanin biosynthesis. Our findings demonstrate a direct functional relationship between ABA signaling and the miR156-SBP-MBW complex regulatory module in driving drought-induced anthocyanin accumulation in grape berries.

High CO2 adaptation mechanisms revealed in the miR156-regulated flowering time pathway

Elevated CO2 concentrations have been observed to accelerate flowering time in Arabidopsis through the action of a highly conserved regulatory network controlled by miR156 and miR172. However, the network's robustness to the impact of increasing CO2 concentrations on flowering time remains poorly understood. In this study, we investigate this question by conducting a comprehensive analysis of the global landscape of network dynamics, including quantifying the probabilities associated with juvenile and flowering states and assessing the speed of the transition between them. Our findings reveal that a CO2 concentration range of 400-800ppm only mildly advances flowering time, contrasting with the dramatic changes from 200 to 300ppm. Notably, the feedback regulation of miR156 by squamosal promoter binding protein-like proteins (SPLs) plays a substantial role in mitigating the effects of increasing CO2 on flowering time. Intriguingly, we consistently observe a correlation between delayed flowering time and increased variance in flowering time, and vice versa, suggesting that this might be an intrinsic adaptation mechanism embedded within the network. To gain a deeper understanding of this network's dynamics, we identified the sensitive features within the feedback loops of miR156 SPLs and miR172-APETALA2 family proteins (AP2s), with the latter proving to be the most sensitive. Strikingly, our study underscores the indispensability of all feedback regulations in maintaining both juvenile and adult states as well as the transition time between them. Together, our research provides the first physical basis in plant species, aiding in the elucidation of novel regulatory mechanisms and the robustness of the miRNAs-regulated network in response to increasing CO2, therefore influencing the control of flowering time. Moreover, this study provides a promising strategy for engineering plant flowering time to enhance their adaptation and resilience.

Versatile Roles of Microbes and Small RNAs in Rice and Planthopper Interactions

Planthopper infestation in rice causes direct and indirect damage through feeding and viral transmission. Host microbes and small RNAs (sRNAs) play essential roles in regulating biological processes, such as metabolism, development, immunity, and stress responses in eukaryotic organisms, including plants and insects. Recently, advanced metagenomic approaches have facilitated investigations on microbial diversity and its function in insects and plants, highlighting the significance of microbiota in sustaining host life and regulating their interactions with the environment. Recent research has also suggested significant roles for sRNA-regulated genes during rice-planthopper interactions. The response and behavior of the rice plant to planthopper feeding are determined by changes in the host transcriptome, which might be regulated by sRNAs. In addition, the roles of microbial symbionts and sRNAs in the host response to viral infection are complex and involve defense-related changes in the host transcriptomic profile. This review reviews the structure and potential functions of microbes and sRNAs in rice and the associated planthopper species. In addition, the involvement of the microbiota and sRNAs in the rice-planthopper-virus interactions during planthopper infestation and viral infection are discussed.

Identification of miRNAs and their target genes associated with improved maize seed vigor induced by gibberellin

High seed vigor is crucial for agricultural production owing to its potential in high quality and yield of crops and a better understanding of the molecular mechanism associated with maize seed vigor is highly necessary. To better understand the involvement and regulatory mechanism of miRNAs correlated with maize seed vigor, small RNAs and degradome sequencing of two inbred lines Yu537A and Yu82 were performed. A total of 791 mature miRNAs were obtained with different expressions, among of which 505 miRNAs were newly identified and the rest miRNAs have been reported before by comparing the miRNAs with the sequences in miRbase database. Analysis of miRNA families showed maize seeds contain fewer miRNA families and larger miRNA families compared with animals, indicating that functions of miRNAs in maize seeds were more synergistic than animals. Degradome sequencing was used to identify the targets of miRNAs and the results showed a total of 6,196 targets were obtained. Function analysis of differentially expressed miRNAs and targets showed Glycan degradation and galactose metabolism were closely correlated with improved maize seed vigor. These findings provide valuable information to understand the involvement of miRNAs with maize seed vigor and these putative genes will be valuable resources for improving the seed vigor in future maize breeding.

Biogenesis, Functions, Interactions, and Resources of Non-Coding RNAs in Plants

Plant transcriptomes encompass a large number of functional non-coding RNAs (ncRNAs), only some of which have protein-coding capacity. Since their initial discovery, ncRNAs have been classified into two broad categories based on their biogenesis and mechanisms of action, housekeeping ncRNAs and regulatory ncRNAs. With advances in RNA sequencing technology and computational methods, bioinformatics resources continue to emerge and update rapidly, including workflow for in silico ncRNA analysis, up-to-date platforms, databases, and tools dedicated to ncRNA identification and functional annotation. In this review, we aim to describe the biogenesis, biological functions, and interactions with DNA, RNA, protein, and microorganism of five major regulatory ncRNAs (miRNA, siRNA, tsRNA, circRNA, lncRNA) in plants. Then, we systematically summarize tools for analysis and prediction of plant ncRNAs, as well as databases. Furthermore, we discuss the silico analysis process of these ncRNAs and present a protocol for step-by-step computational analysis of ncRNAs. In general, this review will help researchers better understand the world of ncRNAs at multiple levels.

Small RNAs Participate in Plant-Virus Interaction and Their Application in Plant Viral Defense

Small RNAs are significant regulators of gene expression, which play multiple roles in plant development, growth, reproductive and stress response. It is generally believed that the regulation of plants’ endogenous genes by small RNAs has evolved from a cellular defense mechanism for RNA viruses and transposons. Most small RNAs have well-established roles in the defense response, such as viral response. During viral infection, plant endogenous small RNAs can direct virus resistance by regulating the gene expression in the host defense pathway, while the small RNAs derived from viruses are the core of the conserved and effective RNAi resistance mechanism. As a counter strategy, viruses evolve suppressors of the RNAi pathway to disrupt host plant silencing against viruses. Currently, several studies have been published elucidating the mechanisms by which small RNAs regulate viral defense in different crops. This paper reviews the distinct pathways of small RNAs biogenesis and the molecular mechanisms of small RNAs mediating antiviral immunity in plants, as well as summarizes the coping strategies used by viruses to override this immune response. Finally, we discuss the current development state of the new applications in virus defense based on small RNA silencing.

Plant RNA-mediated gene regulatory network

Not all transcribed RNAs are protein-coding RNAs. Many of them are non-protein-coding RNAs in diverse eukaryotes. However, some of them seem to be non-functional and are resulted from spurious transcription. A lot of non-protein-coding transcripts have a significant function in the translation process. Gene expressions depend on complex networks of diverse gene regulatory pathways. Several non-protein-coding RNAs regulate gene expression in a sequence-specific system either at the transcriptional level or post-transcriptional level. They include a significant part of the gene expression regulatory network. RNA-mediated gene regulation machinery is evolutionarily ancient. They well-evolved during the evolutionary time and are becoming much more complex than had been expected. In this review, we are trying to summarizing the current knowledge in the field of RNA-mediated gene silencing.

Current status of microRNA-mediated regulation of drought stress responses in cereals

Drought is one of the most important abiotic stress factors impeding crop productivity. With the uncovering of their role as potential regulators of gene expression, microRNAs (miRNAs) have been recognized as new targets for developing stress resistance. MicroRNAs are small noncoding RNAs whose abundance is significantly altered under stress conditions. Interestingly, plant miRNAs predominantly targets transcription factors (TFs), and some of which are also the most critical drought-responsive genes that in turn could regulate the expression of numerous loci with drought-adaptive potential. The phytohormone ABA plays important roles in regulating stomatal conductance and in initiating an adaptive response to drought stress. miRNAs are implicated in regulating ABA-(abscisic acid) and non-ABA-mediated drought resistance pathways. For instance, miR159-MYB module and miR169-NFYA module participates in an ABA-dependent pathway, whereas several other ABA-independent miRNA-target modules (miR156-SPL; miR393-TIR1; miR160-ARF10, ARF16, ARF17; miR167-ARF6 and ARF8; miR390/TAS3siRNA-ARF2, ARF3, ARF4) collectively regulate drought responses in plants. Overall, miRNA-mediated drought response manifests diverse molecular, biochemical and physiological processes. Because of their immense role in controlling gene expression, miRNA manipulation has significant potential to augment plant tolerance to drought stress. This review compiles the current understanding of drought-responsive miRNAs in major cereals. Also, potential miRNA manipulation strategies currently in use along with the challenges and future perspectives are discussed.

The Roles of Long Non-coding RNA in Osteoporosis

The Human Genome Project (HGP) announced in 2001 that it had sequenced the entire human genome, yielding nearly complete human DNA. About 98.5 percent of the human genome has been found to be non-coding sequences. Long non-coding RNA (lncRNA) is a non-coding RNA with a length between 200 and 100,000 nucleotide units. Because of shallow research on lncRNA, it was believed that it had no biological functions, but exists as a by-product of the transcription process. With the development of high-throughput sequencing technology, studies have shown that lncRNA plays important roles in many processes by participating in epigenetics, transcription, translation and protein modification. Current researches have shown that lncRNA also has an important part in the pathogenesis of osteoporosis. Osteoporosis is a common disorder of bone metabolism, also a major medical and socioeconomic challenge worldwide. It is characterized by a systemic reduction in bone mass and microstructure changes, which increases the risk of brittle fractures. It is more common in postmenopausal women and elderly men. However, the roles of lncRNA and relevant mechanisms in osteoporosis remain unclear. Based on this background, we hereby review the roles of lncRNA in osteoporosis, and how it influences the functions of osteoblasts and osteoclasts, providing reference to clinical diagnosis, treatment and prognosis of osteoporosis.

Bioinformatics Resources for Plant Abiotic Stress Responses: State of the Art and Opportunities in the Fast Evolving -Omics Era

Abiotic stresses are among the principal limiting factors for productivity in agriculture. In the current era of continuous climate changes, the understanding of the molecular aspects involved in abiotic stress response in plants is a priority. The rise of -omics approaches provides key strategies to promote effective research in the field, facilitating the investigations from reference models to an increasing number of species, tolerant and sensitive genotypes. Integrated multilevel approaches, based on molecular investigations at genomics, transcriptomics, proteomics and metabolomics levels, are now feasible, expanding the opportunities to clarify key molecular aspects involved in responses to abiotic stresses. To this aim, bioinformatics has become fundamental for data production, mining and integration, and necessary for extracting valuable information and for comparative efforts, paving the way to the modeling of the involved processes. We provide here an overview of bioinformatics resources for research on plant abiotic stresses, describing collections from -omics efforts in the field, ranging from raw data to complete databases or platforms, highlighting opportunities and still open challenges in abiotic stress research based on -omics technologies.

Integrated Analysis of Small RNA, Transcriptome and Degradome Sequencing Provides New Insights into Floral Development and Abscission in Yellow Lupine (Lupinus luteus L.)

The floral development in an important legume crop yellow lupine (Lupinus luteus L., Taper cv.) is often affected by the abscission of flowers leading to significant economic losses. Small non-coding RNAs (sncRNAs), which have a proven effect on almost all developmental processes in other plants, might be of key players in a complex net of molecular interactions regulating flower development and abscission. This study represents the first comprehensive sncRNA identification and analysis of small RNA, transcriptome and degradome sequencing data in lupine flowers to elucidate their role in the regulation of lupine generative development. As shedding in lupine primarily concerns flowers formed at the upper part of the inflorescence, we analyzed samples from extreme parts of raceme separately and conducted an additional analysis of pedicels from abscising and non-abscising flowers where abscission zone forms. A total of 394 known and 28 novel miRNAs and 316 phased siRNAs were identified. In flowers at different stages of development 59 miRNAs displayed differential expression (DE) and 46 DE miRNAs were found while comparing the upper and lower flowers. Identified tasiR-ARFs were DE in developing flowers and were strongly expressed in flower pedicels. The DEmiR-targeted genes were preferentially enriched in the functional categories related to carbohydrate metabolism and plant hormone transduction pathways. This study not only contributes to the current understanding of how lupine flowers develop or undergo abscission but also holds potential for research aimed at crop improvement.

Polymorphisms and genetic effects of PRLR, MOGAT1, MINPP1 and CHUK genes on milk fatty acid traits in Chinese Holstein

BACKGROUND

Our initial genome-wide association study (GWAS) identified 20 promising candidate genes for milk fatty acid (FA) traits in a Chinese Holstein population, including PRLR, MOGAT1, MINPP1 and CHUK genes. In this study, we performed whether they had significant genetic effects on milk FA traits in Chinese Holstein.

RESULTS

We re-sequenced the entire exons and 3000 bp of the 5' and 3' flanking regions, and identified 11 single nucleotide polymorphisms (SNPs), containing four in PRLR, two in MOGAT1, two in MINPP1, and three in CHUK. The SNP-based association analyses showed that all the 11 SNPs were significantly associated with at least one milk FA trait (P = 0.0456 ~ < 0.0001), and none of them had association with C11:0, C13:0, C15:0 and C16:0 (P > 0.05). By the linkage disequilibrium (LD) analyses, we found two, one, one, and one haplotype blocks in PRLR, MOGAT1, MINPP1, and CHUK, respectively, and each haplotype block was significantly associated with at least one milk FA trait (P = 0.0456 ~ < 0.0001). Further, g.38949011G > A in PRLR, and g.111599360A > G and g.111601747 T > A in MOGAT1 were predicted to alter the transcription factor binding sites (TFBSs). A missense mutation, g.39115344G > A, could change the PRLR protein structure. The g.20966385C > G of CHUK varied the binding sequences for microRNAs. Therefore, we deduced the five SNPs as the potential functional mutations.

CONCLUSION

In summary, we first detected the genetic effects of PRLR, MOGAT1, MINPP1 and CHUK genes on milk FA traits, and researched the potential functional mutations. These data provided the basis for further investigation on function validation of the four genes in Chinese Holstein.

Genome-wide identification and comparative analysis of drought-related microRNAs in two maize inbred lines with contrasting drought tolerance by deep sequencing

Drought has become one of the most serious abiotic stresses influencing crop production worldwide. Understanding the molecular regulatory networks underlying drought adaption and tolerance in crops is of great importance for future breeding. microRNAs (miRNAs), as important components of post-transcriptional regulation, play crucial roles in drought response and adaptation in plants. Here, we report a miRNome analysis of two maize inbred lines with contrasting levels of drought tolerance under soil drought in the field. Differential expression analysis showed 11 and 34 miRNAs were uniquely responded to drought in H082183 (drought tolerant) and Lv28 (drought sensitive), respectively, in leaves. In roots, 19 and 23 miRNAs uniquely responded to drought in H082183 and Lv28, respectively. Expression analysis of these drought-responsive miRNA-mRNA modules revealed miR164-MYB, miR164-NAC, miR159-MYB, miR156-SPL and miR160-ARF showed a negative regulatory relationship. Further analysis showed that the miR164-MYB and miR164-NAC modules in the tolerant line modulated the stress response in an ABA (abscisic acid)-dependent manner, while the miR156-SPL and miR160-ARF modules in the sensitive line participated in the inhibition of metabolism in drought-exposed leaves. Together, our results provide new insight into not only drought-tolerance-related miRNA regulation networks in maize but also key miRNAs for further characterization and improvement of maize drought tolerance.

Genome-wide identification of drought-responsive microRNAs in two sets of Malus from interspecific hybrid progenies

Drought stress can negatively impact apple fruit quality and yield. Apple microRNAs (miRNAs) participate in apple tree and fruit development, as well as in biotic stress tolerance; however, it is largely unknown whether these molecules are involved in the drought response. To identify drought-responsive miRNAs in Malus, we first examined the drought stress tolerance of ten F₁ progenies of R3 (M. × domestica) × M. sieversii. We performed Illumina sequencing on pooled total RNA from both drought-tolerant and drought-sensitive plants. The sequencing results identified a total of 206 known miRNAs and 253 candidate novel miRNAs from drought-tolerant plants and drought-sensitive plants under control or drought conditions. We identified 67 miRNAs that were differentially expressed in drought-tolerant plants compared with drought-sensitive plants under drought conditions. Under drought stress, 61 and 35 miRNAs were differentially expressed in drought-tolerant and drought-sensitive plants, respectively. We determined the expression levels of seven out of eight miRNAs by stem-loop qPCR analysis. We also predicted the target genes of all differentially expressed miRNAs and identified the expression of some genes. Gene Ontology analyses indicated that the target genes were mainly involved in stimulus response and cellular and metabolic processes. Finally, we confirmed roles of two miRNAs in apple response to mannitol. Our results reveal candidate miRNAs and their associated mRNAs that could be targeted for improving drought tolerance in Malus species, thus providing a foundation for understanding the molecular networks involved in the response of apple trees to drought stress.

In Silico Identification of Conserved MiRNAs from Physcomitrella patens ESTs and their Target Characterization

Background:

MicroRNAs (miRNAs) are groups of small non-protein-coding endogenous single stranded RNAs with approximately 18-24 nucleotides in length. High evolutionary sequence conservation of miRNAs among plant species and availability of powerful computational tools allow identification of new orthologs and paralogs.

Methods:

New conserved miRNAs in P. patens were found by EST-based homology search approaches. All candidates were screened according to a series of miRNA filtering criteria. Unigene, DFCI Gene Index (PpspGI) databases and psRNATarget algorithm were applied to identify target transcripts using P. patens putative conserved miRNA sequences.

Results:

Nineteen conserved P. patens miRNAs were identified. The sequences were homologous to known reference plant mature miRNA from 10 miRNA families. They could be folded into the typical miRNA secondary structures. RepeatMasker algorithm demonstrated that ppt-miR2919e and pptmiR1533 had simple sequence repeats in their sequences. Target sites (49 genes) were identified for 7 out of 19 miRNAs. GO and KEGG analysis of targets indicated the involvement of some in important multiple biological and metabolic processes.

Conclusion:

The majority of the registered miRNAs in databases were predicted by computational approaches while many more have remained unknown. Due to the conserved nature of miRNAs in plant species from closely to distantly related, homology search-based approaches between plants species could lead to the identification of novel miRNAs in other plant species providing baseline information for further search about the biological functions and evolution of miRNAs.

Plant Small Non-coding RNAs and Their Roles in Biotic Stresses

Non-coding RNAs (ncRNAs) have emerged as critical components of gene regulatory networks across a plethora of plant species. In particular, the 20-30 nucleotide small ncRNAs (sRNAs) play important roles in mediating both developmental processes and responses to biotic stresses. Based on variation in their biogenesis pathways, a number of different sRNA classes have been identified, and their specific functions have begun to be characterized. Here, we review the current knowledge of the biogenesis of the primary sRNA classes, microRNA (miRNA) and small nuclear RNA (snRNA), and their respective secondary classes, and discuss the roles of sRNAs in plant-pathogen interactions. sRNA mobility between species is also discussed along with potential applications of sRNA-plant-pathogen interactions in crop improvement technologies.

Wheat miRNA TaemiR408 Acts as an Essential Mediator in Plant Tolerance to Pi Deprivation and Salt Stress via Modulating Stress-Associated Physiological Processes

MicroRNAs (miRNA) families act as critical regulators for plant growth, development, and responses to abiotic stresses. In this study, we characterized TaemiR408, a miRNA family member of wheat (Triticum aestivum), for the role in mediating plant responses to Pi starvation and salt stress. TaemiR408 targets six genes that encode proteins involving biochemical metabolism, microtubule organization, and signaling transduction. 5'- and 3'-RACE analyses confirmed the mRNA cleavage of target genes mediated by this wheat miRNA. TaemiR408 showed induced expression patterns upon Pi starvation and salt stress and whose upregulated expression was gradually repressed by the normal recovery treatments. The target genes of TaemiR408 exhibited reverse expression patterns to this miRNA, whose transcripts were downregulated under Pi starvation and salt stress and the reduced expression was recovered by the followed normal condition. These results suggest the regulation of the target genes under TaemiR408 through a cleavage mechanism. Tobacco lines with TaemiR408 overexpression exhibited enhanced stress tolerance, showing improved phenotype, biomass, and photosynthesis behavior compared with wild type under both Pi starvation and salt treatments, which closely associate increased P accumulation upon Pi deprivation and elevated osmolytes under salt stress, respectively. Phosphate transporter (PT) gene NtPT2 displays upregulated transcripts in the Pi-deprived TaemiR408 overexpressors; knockdown of this PT gene reduces Pi acquisition under low-Pi stress, confirming its role in improving plant Pi taken up. Likewise, NtPYL2 and NtSAPK3, genes encoding abscisic acid (ABA) receptor and SnRK2 protein, respectively, exhibited upregulated transcripts in salt-challenged TaemiR408 overexpressors; knockdown of them caused deteriorated growth and lowered osmolytes amounts of plants upon salt treatment. Thus, TaemiR408 is crucial for plant adaptations to Pi starvation and salt stress through regulating Pi acquisition under low-Pi stress and remodel ABA signaling pathway and osmoprotects biosynthesis under salt stress.

An expression atlas of miRNAs in Arabidopsis thaliana

MicroRNAs (miRNAs) are small non-coding RNAs that regulate a variety of biological processes. MiRNA expression often exhibits spatial and temporal specificity. However, genome-wide miRNA expression patterns in different organs during development of Arabidopsis thaliana have not yet been systemically investigated. In this study, we sequenced small RNA libraries generated from 27 different organ/tissue types, which cover the entire life cycle of Arabidopsis. Analysis of the sequencing data revealed that most miRNAs are ubiquitously expressed, whereas a small set of miRNAs display highly specific expression patterns. In addition, different miRNA members within the same family have distinct spatial and temporal expression patterns. Moreover, we found that some miRNAs are produced from different arms of their hairpin precursors at different developmental stages. This work provides new insights into the regulation of miRNA biogenesis and a rich resource for future investigation of miRNA functions in Arabidopsis.

Small RNA Activity in Archeological Barley Shows Novel Germination Inhibition in Response to Environment

The recovery of ancient RNA from archeological material could enable the direct study of microevolutionary processes. Small RNAs are a rich source of information because their small size is compatible with biomolecular preservation, and their roles in gene regulation make them likely foci of evolutionary change. We present here the small RNA fraction from a sample of archeological barley generated using high-throughput sequencing that has previously been associated with localized adaptation to drought. Its microRNA profile is broadly similar to 19 globally distributed modern barley samples with the exception of three microRNAs (miRNA159, miRNA319, and miR396), all of which are known to have variable expression under stress conditions. We also found retrotransposon activity to be significantly reduced in the archeological barley compared with the controls, where one would expect the opposite under stress conditions. We suggest that the archeological barley's conflicting stress signals could be the result of long-term adaptation to its local environment.

Genome-wide identification of Hami melon miRNAs with putative roles during fruit development

MicroRNAs represent a family of small endogenous, non-coding RNAs that play critical regulatory roles in plant growth, development, and environmental stress responses. Hami melon is famous for its attractive flavor and excellent nutritional value, however, the mechanisms underlying the fruit development and ripening remains largely unknown. Here, we performed small RNA sequencing to investigate the roles of miRNAs during Hami melon fruit development. Two batches of flesh samples were collected at four fruit development stages. Small RNA sequencing yielded a total of 54,553,424 raw reads from eight libraries. 113 conserved miRNAs belonging to 30 miRNA families and nine novel miRNAs comprising nine miRNA families were identified. The expression of 42 conserved miRNAs and three Hami melon-specific miRNAs significantly changed during fruit development. Furthermore, 484 and 124 melon genes were predicted as putative targets of 29 conserved and nine Hami melon-specific miRNA families, respectively. GO enrichment analysis were performed on target genes, "transcription, DNA-dependent", "rRNA processing", "oxidation reduction", "signal transduction", "regulation of transcription, DNA-dependent", and "metabolic process" were the over-represented biological process terms. Cleavage sites of six target genes were validated using 5' RACE. Our results present a comprehensive set of identification and characterization of Hami melon fruit miRNAs and their potential targets, which provide valuable basis towards understanding the regulatory mechanisms in programmed process of normal Hami fruit development and ripening. Specific miRNAs could be selected for further research and applications in breeding practices.

Renaissance in phytomedicines: promising implications of NGS technologies

Medicinal plant research is growing significantly in faith to discover new and more biologically compatible phytomedicines. Deposition of huge genome/trancriptome sequence data assisted by NGS technologies has revealed the new possibilities for producing upgraded bioactive molecules in medicinal plants. Growing interest of investors and consumers in the herbal drugs raises the need for extensive research to open the facts and details of every inch of life canvas of medicinal plants to produce improved quality of phytomedicines. As in agriculture crops, knowledge emergence from medicinal plant's genome/transcriptome, can be used to assure their amended quality and these improved varieties are then transported to the fields for cultivation. Genome studies generate huge sequence data which can be exploited further for obtaining information regarding genes/gene clusters involved in biosynthesis as well as regulation. This can be achieved rapidly at a very large scale with NGS platforms. Identification of new RNA molecules has become possible, which can lead to the discovery of novel compounds. Sequence information can be combined with advanced phytochemical and bioinformatics tools to discover functional herbal drugs. Qualitative and quantitative analysis of small RNA species put a light on the regulatory aspect of biosynthetic pathways for phytomedicines. Inter or intra genomic as well as transcriptomic interactive processes for biosynthetic pathways can be elucidated in depth. Quality management of herbal material will also become rapid and high throughput. Enrichment of sequence information will be used to engineer the plants to get more efficient phytopharmaceuticals. The present review comprises of role of NGS technologies to boost genomic studies of pharmaceutically important plants and further, applications of sequence information aiming to produce enriched phytomedicines. Emerging knowledge from the medicinal plants genome/transcriptome can give birth to deep understanding of the processes responsible for biosynthesis of medicinally important compounds.

Physiological and transcriptional responses of two contrasting Populus clones to nitrogen stress

The aim of this study was to reveal mechanisms responsible for nitrogen (N) stress in two contrasting Populus clones. Leaves of Nanlin 1388 (N stress-insensitive clone hybrids of Populus deltoides Bart.CV. × Populus euramericana (Dode) Guineir CV) and Nanlin 895 (N stress-sensitive clone hybrids of Populus deltoides Bart.CV. × Populus euramericana (Dode) Guineir CV) were harvested and analyzed. Different responses visible in photosynthesis, N and carbon contents, physiological traits, and chlorophyll were observed. The Solexa/Illumina's digital gene expression system was used to investigate differentially expressed miRNAs and mRNAs under N stress. Target profiling, and biological network and function analyses were also performed. Randomly selected mRNAs and miRNAs were validated by quantitative reverse transcription polymerase chain reaction. In all, 110 Nanlin 1388 and 122 Nanlin 895 miRNAs were differentially expressed, among which 34 and 23 miRNAs were newly found in the two clones, respectively. Under N stress, a total of 329 and 98 mRNAs were regulated in N stress-insensitive and -sensitive clones, respectively. Notably, the miR396 family and its regulated mRNAs were altered in both clones under N stress, while miR646 was regulated only in the N stress-insensitive clone (Nanlin 1388), and miR156, miR319 and miR393 in the N stress-sensitive clone (Nanlin 895). Gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses also proved several clone-specific functions and pathways. These findings may be significant for understanding the genetic responses of Populus to N stress.

Characterization of the Small RNA Transcriptome of the Marine Coccolithophorid, Emiliania huxleyi

Small RNAs (smRNAs) control a variety of cellular processes by silencing target genes at the transcriptional or post-transcription level. While extensively studied in plants, relatively little is known about smRNAs and their targets in marine phytoplankton, such as Emiliania huxleyi (E. huxleyi). Deep sequencing was performed of smRNAs extracted at different time points as E. huxleyi cells transition from logarithmic to stationary phase growth in batch culture. Computational analyses predicted 18 E. huxleyi specific miRNAs. The 18 miRNA candidates and their precursors vary in length (18–24 nt and 71–252 nt, respectively), genome copy number (3–1,459), and the number of genes targeted (2–107). Stem-loop real time reverse transcriptase (RT) PCR was used to validate miRNA expression which varied by nearly three orders of magnitude when growth slows and cells enter stationary phase. Stem-loop RT PCR was also used to examine the expression profiles of miRNA in calcifying and non-calcifying cultures, and a small subset was found to be differentially expressed when nutrients become limiting and calcification is enhanced. In addition to miRNAs, endogenous small RNAs such as ra-siRNAs, ta-siRNAs, nat-siRNAs, and piwiRNAs were predicted along with the machinery for the biogenesis and processing of si-RNAs. This study is the first genome-wide investigation smRNAs pathways in E. huxleyi. Results provide new insights into the importance of smRNAs in regulating aspects of physiological growth and adaptation in marine phytoplankton and further challenge the notion that smRNAs evolved with multicellularity, expanding our perspective of these ancient regulatory pathways.

A Primary Sequence Analysis of the ARGONAUTE Protein Family in Plants

Small RNA (sRNA)-mediated gene silencing represents a conserved regulatory mechanism controlling a wide diversity of developmental processes through interactions of sRNAs with proteins of the ARGONAUTE (AGO) family. On the basis of a large phylogenetic analysis that includes 206 AGO genes belonging to 23 plant species, AGO genes group into four clades corresponding to the phylogenetic distribution proposed for the ten family members of Arabidopsis thaliana. A primary analysis of the corresponding protein sequences resulted in 50 sequences of amino acids (blocks) conserved across their linear length. Protein members of the AGO4/6/8/9 and AGO1/10 clades are more conserved than members of the AGO5 and AGO2/3/7 clades. In addition to blocks containing components of the PIWI, PAZ, and DUF1785 domains, members of the AGO2/3/7 and AGO4/6/8/9 clades possess other consensus block sequences that are exclusive of members within these clades, suggesting unforeseen functional specialization revealed by their primary sequence. We also show that AGO proteins of animal and plant kingdoms share linear sequences of blocks that include motifs involved in posttranslational modifications such as those regulating AGO2 in humans and the PIWI protein AUBERGINE in Drosophila. Our results open possibilities for exploring new structural and functional aspects related to the evolution of AGO proteins within the plant kingdom, and their convergence with analogous proteins in mammals and invertebrates.

Novel Meiotic miRNAs and Indications for a Role of PhasiRNAs in Meiosis

Small RNAs (sRNA) add additional layers to the regulation of gene expression, with siRNAs directing gene silencing at the DNA level by RdDM (RNA-directed DNA methylation), and micro RNAs (miRNAs) directing post-transcriptional regulation of specific target genes, mostly by mRNA cleavage. We used manually isolated male meiocytes from maize (Zea mays) to investigate sRNA and DNA methylation landscapes during zygotene, an early stage of meiosis during which steps of meiotic recombination and synapsis of paired homologous chromosomes take place. We discovered two novel miRNAs from meiocytes, zma-MIR11969 and zma-MIR11970, and identified putative target genes. Furthermore, we detected abundant phasiRNAs of 21 and 24 nt length. PhasiRNAs are phased small RNAs which occur in 21 or 24 nt intervals, at a few hundred loci, specifically in male reproductive tissues in grasses. So far, the function of phasiRNAs remained elusive. Data from isolated meiocytes now revealed elevated DNA methylation at phasiRNA loci, especially in the CHH context, suggesting a role for phasiRNAs in cis DNA methylation. In addition, we consider a role of these phasiRNAs in chromatin remodeling/dynamics during meiosis. However, this is not well supported yet and will need more additional data. Here, we only lay out the idea due to other relevant literature and our additional observation of a peculiar GC content pattern at phasiRNA loci. Chromatin remodeling is also indicated by the discovery that histone genes were enriched for sRNA of 22 nt length. Taken together, we gained clues that lead us to hypothesize sRNA-driven DNA methylation and possibly chromatin remodeling during male meiosis in the monocot maize which is in line with and extends previous knowledge.

RNA sequencing and functional analysis implicate the regulatory role of long non-coding RNAs in tomato fruit ripening

Recently, long non-coding RNAs (lncRNAs) have been shown to play critical regulatory roles in model plants, such as Arabidopsis, rice, and maize. However, the presence of lncRNAs and how they function in fleshy fruit ripening are still largely unknown because fleshy fruit ripening is not present in the above model plants. Tomato is the model system for fruit ripening studies due to its dramatic ripening process. To investigate further the role of lncRNAs in fruit ripening, it is necessary and urgent to discover and identify novel lncRNAs and understand the function of lncRNAs in tomato fruit ripening. Here it is reported that 3679 lncRNAs were discovered from wild-type tomato and ripening mutant fruit. The lncRNAs are transcribed from all tomato chromosomes, 85.1% of which came from intergenic regions. Tomato lncRNAs are shorter and have fewer exons than protein-coding genes, a situation reminiscent of lncRNAs from other model plants. It was also observed that 490 lncRNAs were significantly up-regulated in ripening mutant fruits, and 187 lncRNAs were down-regulated, indicating that lncRNAs could be involved in the regulation of fruit ripening. In line with this, silencing of two novel tomato intergenic lncRNAs, lncRNA1459 and lncRNA1840, resulted in an obvious delay of ripening of wild-type fruit. Overall, the results indicated that lncRNAs might be essential regulators of tomato fruit ripening, which sheds new light on the regulation of fruit ripening.

Identification of miRNAs and Their Target Genes Associated with Sweet Corn Seed Vigor by Combined Small RNA and Degradome Sequencing

High seed vigor is significant for agriculture. Low seed vigor of sweet corn hindered the popularization of sweet corn (Zea mays L. saccharata Sturt). To better understand the involvement and regulatory mechanism of miRNAs with seed vigor, small RNA libraries from seeds non-artificially aged and artificially aged for 2 days were generated by small RNA sequencing. A total of 27 differentially expressed miRNAs were discovered, of which 10 were further confirmed by real-time quantitative polymerase chain reaction. Furthermore, targets of miRNAs were identified by degradome sequencing. A total of 1142 targets that were potentially cleaved by 131 miRNAs were identified. Gene ontology (GO) annotations of target transcripts indicated that 26 target genes cleaved by 9 differentially expressed miRNAs might play roles in the regulation of seed vigor, such as peroxidase superfamily protein targeted by PC-5p-213179_17 playing a role in the oxidation-reduction process and response to oxidative stress. These findings provide valuable information to understand the involvement of miRNAs with seed vigor.

Long non-coding RNAs and their biological roles in plants

With the development of genomics and bioinformatics, especially the extensive applications of high-throughput sequencing technology, more transcriptional units with little or no protein-coding potential have been discovered. Such RNA molecules are called non-protein-coding RNAs (npcRNAs or ncRNAs). Among them, long npcRNAs or ncRNAs (lnpcRNAs or lncRNAs) represent diverse classes of transcripts longer than 200 nucleotides. In recent years, the lncRNAs have been considered as important regulators in many essential biological processes. In plants, although a large number of lncRNA transcripts have been predicted and identified in few species, our current knowledge of their biological functions is still limited. Here, we have summarized recent studies on their identification, characteristics, classification, bioinformatics, resources, and current exploration of their biological functions in plants.

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.

Plant systems biology: insights, advances and challenges

Plants dwelling at the base of biological food chain are of fundamental significance in providing solutions to some of the most daunting ecological and environmental problems faced by our planet. The reductionist views of molecular biology provide only a partial understanding to the phenotypic knowledge of plants. Systems biology offers a comprehensive view of plant systems, by employing a holistic approach integrating the molecular data at various hierarchical levels. In this review, we discuss the basics of systems biology including the various 'omics' approaches and their integration, the modeling aspects and the tools needed for the plant systems research. A particular emphasis is given to the recent analytical advances, updated published examples of plant systems biology studies and the future trends.

The effects of carbon dioxide and temperature on microRNA expression in Arabidopsis development

Elevated levels of CO2 and temperature can both affect plant growth and development, but the signalling pathways regulating these processes are still obscure. MicroRNAs function to silence gene expression, and environmental stresses can alter their expressions. Here we identify, using the small RNA-sequencing method, microRNAs that change significantly in expression by either doubling the atmospheric CO2 concentration or by increasing temperature 3-6 °C. Notably, nearly all CO2-influenced microRNAs are affected inversely by elevated temperature. Using the RNA-sequencing method, we determine strongly correlated expression changes between miR156/157 and miR172, and their target transcription factors under elevated CO2 concentration. Similar correlations are also found for microRNAs acting in auxin-signalling, stress responses and potential cell wall carbohydrate synthesis. Our results demonstrate that both CO2 and temperature alter microRNA expression to affect Arabidopsis growth and development, and miR156/157- and miR172-regulated transcriptional network might underlie the onset of early flowering induced by increasing CO2.

The involvement of InMIR167 in the regulation of expression of its target gene InARF8, and their participation in the vegetative and generative development of Ipomoea nil plants

The plant hormone auxin plays a critical role in regulating plant growth and development. Recent advances have been made that having improved our understanding of auxin response pathways, primarily by characterizing the genes encoding auxin response factors (ARFs) in Arabidopsis. In addition, the expression of some ARFs is regulated by microRNAs (miRNAs). In Arabidopsis thaliana, ARF6 and ARF8 are targeted by miR167, whereas ARF10, ARF16 and ARF17 are targeted by miR160. Nevertheless, little is known about any possible interactions between miRNAs and the auxin signaling pathway during plant development. In this study, we isolated the miR167 target gene InARF8 cDNA from the cotyledons of the short day plant (SDP) Ipomoea nil (named also Pharbitis nil). Additionally, the In-miR167 precursor was identified from the I. nil EST database and analyses of InARF8 mRNA, In-pre-miR167 and mature miR167 accumulation in the plant's vegetative and generative organs were performed. The identified cDNA of InARF8 contains a miR167 complementary sequence and shows significant similarity to ARF8 cDNAs of other plant species. The predicted amino acid sequence of InARF8 includes all of the characteristic domains for ARF family transcription factors (B3 DNA-binding domain, AUX/IAA-CTD and a glutamine-rich region). Quantitative RT-PCR reactions and in situ hybridization indicated that InARF8 was expressed primarily in the shoot apices, leaf primordia and hypocotyls of I. nil seedlings, as well as in flower pistils and petals. The InARF8 transcript level increased consistently during the entire period of pistil development, whereas in the stamens, the greatest transcriptional activity occurred only during the intensive elongation phase. Additionally, an expression analysis of both the precursor In-pre-miR167 molecules identified and mature miRNA was performed. We observed that, in most of the organs examined, the InARF8 expression pattern was opposite to that of MIR167, indicating that the gene's activity was regulated by mRNA cleavage. Our findings suggested that InARF8 and InMIR167 participated in the development of young tissues, especially the shoot apices and flower elements. The main function of MIR167 appears to be to regulate InARF8 organ localization.

Fast selection of miRNA candidates based on large-scale pre-computed MFE sets of randomized sequences

Background

Small RNAs are important regulators of genome function, yet their prediction in genomes is still a major computational challenge. Statistical analyses of pre-miRNA sequences indicated that their 2D structure tends to have a minimal free energy (MFE) significantly lower than MFE values of equivalently randomized sequences with the same nucleotide composition, in contrast to other classes of non-coding RNA. The computation of many MFEs is, however, too intensive to allow for genome-wide screenings.

Results

Using a local grid infrastructure, MFE distributions of random sequences were pre-calculated on a large scale. These distributions follow a normal distribution and can be used to determine the MFE distribution for any given sequence composition by interpolation. It allows on-the-fly calculation of the normal distribution for any candidate sequence composition.

Conclusion

The speedup achieved makes genome-wide screening with this characteristic of a pre-miRNA sequence practical. Although this particular property alone will not be able to distinguish miRNAs from other sequences sufficiently discriminative, the MFE-based P-value should be added to the parameters of choice to be included in the selection of potential miRNA candidates for experimental verification.

Exploring the interaction between small RNAs and R genes during Brachypodium response to Fusarium culmorum infection

The present study aims to investigate small RNA interactions with putative disease response genes in the model grass species Brachypodium distachyon. The fungal pathogen Fusarium culmorum (Fusarium herein) and phytohormone salicylic acid treatment were used to induce the disease response in Brachypodium. Initially, 121 different putative disease response genes were identified using bioinformatic and homology based approaches. Computational prediction was used to identify 33 candidate new miRNA coding sequences, of which 9 were verified by analysis of small RNA sequence libraries. Putative Brachypodium miRNA target sites were identified in the disease response genes, and a subset of which were screened for expression and possible miRNA interactions in 5 different Brachypodium lines infected with Fusarium. An NBS-LRR family gene, 1g34430, was polymorphic among the lines, forming two major genotypes, one of which has its miRNA target sites deleted, resulting in altered gene expression during infection. There were siRNAs putatively involved in regulation of this gene, indicating a role of small RNAs in the B. distachyon disease response.

Identification of stable, high copy number, medium-sized RNA degradation intermediates that accumulate in plants under non-stress conditions

It is becoming increasingly evident that the RNA degradome is a crucial component of the total cellular RNA pool. Here, we present an analysis of the medium-sized RNAs (midi RNAs) that form in Arabidopsis thaliana. Our analyses revealed that the midi RNA fraction contained mostly 20-70-nt-long fragments derived from various RNA species, including tRNA, rRNA, mRNA and snRNA. The majority of these fragments could be classified as stable RNA degradation intermediates (RNA degradants). Using two dimensional polyacrylamide gel electrophoresis, we demonstrated that high copy number RNA (hcn RNA) degradants appear in plant cells not only during stress, as it was earlier suggested. They are continuously produced also under physiological conditions. The data collected indicated that the accumulation pattern of the hcn RNA degradants is organ-specific and can be affected by various endogenous and exogenous factors. In addition, we demonstrated that selected degradants efficiently inhibit translation in vitro. Thus, the results of our studies suggest that hcn RNA degradants are likely to be involved in the regulation of gene expression in plants.

Regulation of MIR genes in response to abiotic stress in Hevea brasiliensis

Increasing demand for natural rubber (NR) calls for an increase in latex yield and also an extension of rubber plantations in marginal zones. Both harvesting and abiotic stresses lead to tapping panel dryness through the production of reactive oxygen species. Many microRNAs regulated during abiotic stress modulate growth and development. The objective of this paper was to study the regulation of microRNAs in response to different types of abiotic stress and hormone treatments in Hevea. Regulation of MIR genes differs depending on the tissue and abiotic stress applied. A negative co-regulation between HbMIR398b with its chloroplastic HbCuZnSOD target messenger is observed in response to salinity. The involvement of MIR gene regulation during latex harvesting and tapping panel dryness (TPD) occurrence is further discussed.

Identification and characterization of microRNAs in the developing maize endosperm

MicroRNAs (miRNAs) are non-coding RNAs that are approximately 20-22 nucleotides long. miRNAs have been shown to be important regulators that control a large variety of biological functions in eukaryotic cells. To investigate the roles of miRNAs in maize endosperm development, two small RNA libraries of maize endosperm at two developmental stages were sequenced. A total of 17,773,394 and 18,586,523 small RNA raw reads were obtained, respectively. Further analysis identified and characterized 95 known miRNAs belonging to 20 miRNA families. In addition, 18 novel miRNAs were identified and grouped into 11 families. Potential targets for 5 of the novel miRNA families were successfully predicted. We had also identified 12 corresponding miRNAs* of these novel miRNAs. In summary, we investigated expression patterns of miRNA in maize endosperm at key developmental stages and identified miRNAs that are likely to playing an important role in endosperm development.

Involvement of a novel intronic microRNA in cross regulation of N-methyltransferase genes involved in caffeine biosynthesis in Coffea canephora

There are numerous reports on intronic miRNAs from plants, most of which are involved in the regulation of unrelated genes. Some of the target genes are antagonistic to the host genes. Intronic miRNAs in animal systems, however, are known to have synergistic effects. This article is the first to report a similar regulatory effect of a miRNA originating from an intron in plants. NMT genes involved in caffeine biosynthesis were silenced to obtain transformants with reduced caffeine. Transcript analysis revealed the accumulation of transcripts for a related NMT gene (CaMTL1) in transformants bearing either antisense or RNAi constructs. The altered expression was assumed to relate to the silencing of the NMT genes. Bioinformatics analysis of the genes involved in biosynthesis revealed the presence of an intronic miRNA originating from the intron of the theobromine synthase gene targeting CaMTL1. The putative miRNA was cloned and sequenced. Modified 5'-RLM-RACE mapping of the cleavage site and subsequent Northern blotting experimentally demonstrated the presence and activity of such a miRNA in Coffea canephora. This novel regulatory mechanism previously unreported in plants will shed more light onto the evolution of multigene families and the role of introns in this process.

Genome-wide analysis of 24-nt siRNAs dynamic variations during rice superior and inferior grain filling

24 nt-siRNAs are the most abundant small interfering RNAs in rice grains aside from microRNAs. To investigate the roles that 24 nt-siRNAs played in the poor grain filling of rice inferior grains, dynamic variations of 24 nt-siRNAs in inferior grains were compared with those of superior grains by using small RNA deep sequencing technology. The results showed that 24 nt-siRNAs derived from multiple regions of rice genome, and the maintenance of the two strands of 24 nt-siRNA duplex was a non-random process. The amounts of 24 nt-siRNAs declined with the process of grain filling in both superior and inferior grains, but 24 nt-siRNAs in inferior grains was much higher than that of superior grains in each period we sampled. Bioinformatics prediction indicated that 24 nt-siRNAs targeted on more genes involved in most of the known KEGG rice pathways, such as the starch and sucrose biosynthesis pathway. Combined with digital gene expression profiling of target genes, 24 nt-siRNAs mapped on the antisense strands of exons were specifically investigated, but the abundance of 24 nt-siRNAs did not show negative correlations with their corresponding target genes. The results indicated that 24 nt-siRNAs were not involved in down-regulation of target genes. The potential biological meanings for this inconsistency were probably the results of methylation directed gene expression activation, or competition for small RNA stability methylation.

Transcriptome response to nitrogen starvation in rice

Nitrogen is an essential mineral nutrient required for plant growth and development. Insufficient nitrogen (N) supply triggers extensive physiological and biochemical changes in plants. In this study, we used Affymetrix GeneChip rice genome arrays to analyse the dynamics of rice transcriptome under N starvation. N starvation induced or suppressed transcription of 3518 genes, representing 10.88 percent of the genome. These changes, mostly transient, affected various cellular metabolic pathways, including stress response, primary and secondary metabolism, molecular transport, regulatory process and organismal development. 462 or 13.1 percent transcripts for N starvation expressed similarly in root and shoot. Comparative analysis between rice and Arabidopsis identified 73 orthologous groups that responded to N starvation, demonstrated the existence of conserved N stress coupling mechanism among plants. Additional analysis of transcription profiles of microRNAs revealed differential expression of miR399 and miR530 under N starvation, suggesting their potential roles in plant nutrient homeostasis.

Regulation of miR399f transcription by AtMYB2 affects phosphate starvation responses in Arabidopsis

Although a role for microRNA399 (miR399) in plant responses to phosphate (Pi) starvation has been indicated, the regulatory mechanism underlying miR399 gene expression is not clear. Here, we report that AtMYB2 functions as a direct transcriptional activator for miR399 in Arabidopsis (Arabidopsis thaliana) Pi starvation signaling. Compared with untransformed control plants, transgenic plants constitutively overexpressing AtMYB2 showed increased miR399f expression and tissue Pi contents under high Pi growth and exhibited elevated expression of a subset of Pi starvation-induced genes. Pi starvation-induced root architectural changes were more exaggerated in AtMYB2-overexpressing transgenic plants compared with the wild type. AtMYB2 directly binds to a MYB-binding site in the miR399f promoter in vitro, as well as in vivo, and stimulates miR399f promoter activity in Arabidopsis protoplasts. Transcription of AtMYB2 itself is induced in response to Pi deficiency, and the tissue expression patterns of miR399f and AtMYB2 are similar. Both genes are expressed mainly in vascular tissues of cotyledons and in roots. Our results suggest that AtMYB2 regulates plant responses to Pi starvation by regulating the expression of the miR399 gene.

Exploration of small non coding RNAs in wheat (Triticum aestivum L.)

Large numbers of noncoding RNA transcripts (ncRNAs) are being revealed in animals and plants, which can function at the transcriptional or posttranscriptional level to negatively regulate or control genes, repetitive sequences, viruses, and mobile elements. With the identification of microRNA and siRNAs in diverse organisms, increasing evidences indicate that these short npcRNAs play important roles in development, stress response and diseases by cleavage of target mRNA or interfere with translation of target genes. To explore the small RNA transcriptome in wheat (Triticum aestivum L.), a couple of small RNA libraries were constructed and sequenced by high throughput sequencing method. In this review, we focused on the discovery of wheat small RNAs including miRNA and some other non coding small RNAs, then have a view of miRNAs conservations and differences among wheat and other plant species. We also summarized the developmental and stress responsive expression of wheat miRNAs and these observations could serve as a foundation for future functional studies.

High throughput sequencing reveals novel and abiotic stress-regulated microRNAs in the inflorescences of rice

BACKGROUND

MicroRNAs (miRNAs) are small RNA molecules that play important regulatory roles in plant development and stress responses. Identification of stress-regulated miRNAs is crucial for understanding how plants respond to environmental stimuli. Abiotic stresses are one of the major factors that limit crop growth and yield. Whereas abiotic stress-regulated miRNAs have been identified in vegetative tissues in several plants, they are not well studied in reproductive tissues such as inflorescences.

RESULTS

We used Illumina deep sequencing technology to sequence four small RNA libraries that were constructed from the inflorescences of rice plants that were grown under control condition and drought, cold, or salt stress. We identified 227 miRNAs that belong to 127 families, including 70 miRNAs that are not present in the miRBase. We validated 62 miRNAs (including 10 novel miRNAs) using published small RNA expression data in DCL1, DCL3, and RDR2 RNAi lines and confirmed 210 targets from 86 miRNAs using published degradome data. By comparing the expression levels of miRNAs, we identified 18, 15, and 10 miRNAs that were regulated by drought, cold and salt stress conditions, respectively. In addition, we identified 80 candidate miRNAs that originated from transposable elements or repeats, especially miniature inverted-repeat elements (MITEs).

CONCLUSION

We discovered novel miRNAs and stress-regulated miRNAs that may play critical roles in stress response in rice inflorescences. Transposable elements or repeats, especially MITEs, are rich sources for miRNA origination.

Characterization of microRNAs expression during maize seed development

Background

MicroRNAs (miRNAs) are approximately 20-22 nt non-coding RNAs that play key roles in many biological processes in both animals and plants. Although a number of miRNAs were identified in maize, the function of miRNA in seed development was merely discussed.

Results

In this study, two small RNA libraries were sequenced, and a total reads of 9,705,761 and 9,005,563 were generated from developing seeds and growing leaves, respectively. Further analysis identified 125 known miRNAs in seeds and 127 known miRNAs in leaves. 54 novel miRNAs were identified and they were not reported in other plants. Additionally, some miRNA*s of these novel miRNAs were detected. Potential targets of all novel miRNAs were predicted based on our strict criteria. In addition to deep-sequencing, miRNA microarray study confirmed the higher expression of several miRNAs in seeds. In summary, our results indicated the distinct expression of miRNAs during seed development.

Conclusions

We had identified 125 and 127 known miRNAs from seeds and leaves in maize, and a total of 54 novel miRNAs were discovered. The different miRNA expression profile in developing seeds were revealed by both sequencing and microarray studies.

Characterization and differential expression of microRNAs elicited by sulfur deprivation in Chlamydomonas reinhardtii

Background

microRNAs (miRNAs) have been found to play an essential role in the modulation of numerous biological processes in eukaryotes. Chlamydomonas reinhardtii is an ideal model organism for the study of many metabolic processes including responses to sulfur-deprivation. We used a deep sequencing platform to extensively profile and identify changes in the miRNAs expression that occurred under sulfur-replete and sulfur-deprived conditions. The aim of our research was to characterize the differential expression of Chlamydomonas miRNAs under sulfur-deprived conditions, and subsequently, the target genes of miRNA involved in sulfur-deprivation were further predicted and analyzed.

Results

By using high-throughput sequencing, we characterized the microRNA transcriptomes under sulphur-replete and sulfur-deprived conditions in Chlamydomonas reinhardtii. We predicted a total of 310 miRNAs which included 85 known miRNAs and 225 novel miRNAs. 13 miRNAs were the specific to the sulfur-deprived conditions. 47 miRNAs showed significantly differential expressions responding to sulfur-deprivation, and most were up-regulated in the small RNA libraries with sulfur-deprivation. Using a web-based integrated system (Web MicroRNAs Designer 3) and combing the former information from a transcriptome of Chlamydomonas reinhardtii, 22 miRNAs and their targets involved in metabolism regulation with sulfur-deprivation were verified.

Conclusions

Our results indicate that sulfur-deprivation may have a significant influence on small RNA expression patterns, and the differential expressions of miRNAs and interactions between miRNA and its targets might further reveal the molecular mechanism responding to sulfur-deprivation in Chlamydomonas reinhardtii.

Identification of novel microRNAs in Hevea brasiliensis and computational prediction of their targets

BACKGROUND

Plants respond to external stimuli through fine regulation of gene expression partially ensured by small RNAs. Of these, microRNAs (miRNAs) play a crucial role. They negatively regulate gene expression by targeting the cleavage or translational inhibition of target messenger RNAs (mRNAs). In Hevea brasiliensis, environmental and harvesting stresses are known to affect natural rubber production. This study set out to identify abiotic stress-related miRNAs in Hevea using next-generation sequencing and bioinformatic analysis.

RESULTS

Deep sequencing of small RNAs was carried out on plantlets subjected to severe abiotic stress using the Solexa technique. By combining the LeARN pipeline, data from the Plant microRNA database (PMRD) and Hevea EST sequences, we identified 48 conserved miRNA families already characterized in other plant species, and 10 putatively novel miRNA families. The results showed the most abundant size for miRNAs to be 24 nucleotides, except for seven families. Several MIR genes produced both 20-22 nucleotides and 23-27 nucleotides. The two miRNA class sizes were detected for both conserved and putative novel miRNA families, suggesting their functional duality. The EST databases were scanned with conserved and novel miRNA sequences. MiRNA targets were computationally predicted and analysed. The predicted targets involved in "responses to stimuli" and to "antioxidant" and "transcription activities" are presented.

CONCLUSIONS

Deep sequencing of small RNAs combined with transcriptomic data is a powerful tool for identifying conserved and novel miRNAs when the complete genome is not yet available. Our study provided additional information for evolutionary studies and revealed potentially specific regulation of the control of redox status in Hevea.