Genome-wide analysis for discovery of rice microRNAs reveals natural antisense microRNAs (nat-miRNAs)

PEP444c encoded within the MIR444c gene regulates microRNA444c accumulation in barley

MicroRNAs are small, noncoding RNA molecules that regulate the expression of their target genes. The MIR444 gene family is present exclusively in monocotyledons, and microRNAs444 from this family have been shown to target certain MADS-box transcription factors in rice and barley. We identified three barley MIR444 (MIR444a/b/c) genes and comprehensively characterised their structure and the processing pattern of the primary transcripts (pri-miRNAs444). Pri-microRNAs444 undergo extensive alternative splicing, generating functional and nonfunctional pri-miRNA444 isoforms. We show that barley pri-miRNAs444 contain numerous open reading frames (ORFs) whose transcripts associate with ribosomes. Using specific antibodies, we provide evidence that selected ORFs encoding PEP444a within MIR444a and PEP444c within MIR444c are expressed in barley plants. Moreover, we demonstrate that CRISPR-associated endonuclease 9 (Cas9)-mediated mutagenesis of the PEP444c-encoding sequence results in a decreased level of PEP444 transcript in barley shoots and roots and a 5-fold reduced level of mature microRNA444c in roots. Our observations suggest that PEP444c encoded by the MIR444c gene is involved in microRNA444c biogenesis in barley.

Revitalizing miRNAs mediated agronomical advantageous traits improvement in rice

One of the key enigmas in conventional and modern crop improvement programmes is how to introduce beneficial traits without any penalty impairment. Rice (Oryza sativa L.), among the essential staple food crops grown and utilized worldwide, needs to improve genotypes in multifaceted ways. With the global view to feed ten billion under the climatic perturbation, only a potent functional master regulator can withstand with hope for the next green revolution and food security. miRNAs are such, miniature, fine tuners for crop improvement and provide a value addition in emerging technologies, namely large-scale genotyping, phenotyping, genome editing, marker-assisted selection, and genomic selection, to make rice production feasible. There has been surplus research output generated since the last decade on miRNAs in rice, however, recent functional knowledge is limited to reaping the benefits for conventional and modern improvements in rice to avoid ambiguity and redundancy in the generated data. Here, we present the latest functional understanding of miRNAs in rice. In addition, their biogenesis, intra- and inter-kingdom signaling and communication, implication of amiRNAs, and consequences upon integration with CRISPR-Cas9. Further, highlights refer to the application of miRNAs for rice agronomical trait improvements, broadly classified into three functional domains. The majority of functionally established miRNAs are responsible for growth and development, followed by biotic and abiotic stresses. Tabular cataloguing reveals and highlights two multifaceted modules that were extensively studied. These belong to miRNA families 156 and 396, orchestrate multifarious aspects of advantageous agronomical traits. Moreover, updated and exhaustive functional aspects of different supplemental miRNA modules that would strengthen rice improvement are also being discussed.

Weighted gene co-expression network analysis of nitrogen (N)-responsive genes and the putative role of G-quadruplexes in N use efficiency (NUE) in rice

Rice is an important target to improve crop nitrogen (N) use efficiency (NUE), and the identification and shortlisting of the candidate genes are still in progress. We analyzed data from 16 published N-responsive transcriptomes/microarrays to identify, eight datasets that contained the maximum number of 3020 common genes, referred to as N-responsive genes. These include different classes of transcription factors, transporters, miRNA targets, kinases and events of post-translational modifications. A Weighted gene co-expression network analysis (WGCNA) with all the 3020 N-responsive genes revealed 15 co-expression modules and their annotated biological roles. Protein-protein interaction network analysis of the main module revealed the hub genes and their functional annotation revealed their involvement in the ubiquitin process. Further, the occurrences of G-quadruplex sequences were examined, which are known to play important roles in epigenetic regulation but are hitherto unknown in N-response/NUE. Out of the 3020 N-responsive genes studied, 2298 contained G-quadruplex sequences. We compared these N-responsive genes containing G-quadruplex sequences with the 3601 genes we previously identified as NUE-related (for being both N-responsive and yield-associated). This analysis revealed 389 (17%) NUE-related genes containing G-quadruplex sequences. These genes may be involved in the epigenetic regulation of NUE, while the rest of the 83% (1811) genes may regulate NUE through genetic mechanisms and/or other epigenetic means besides G-quadruplexes. A few potentially important genes/processes identified as associated with NUE were experimentally validated in a pair of rice genotypes contrasting for NUE. The results from the WGCNA and G4 sequence analysis of N-responsive genes helped identify and shortlist six genes as candidates to improve NUE. Further, the hitherto unavailable segregation of genetic and epigenetic gene targets could aid in informed interventions through genetic and epigenetic means of crop improvement.

Identification of microRNAs That Provide a Low Light Stress Tolerance-Mediated Signaling Pathway during Vegetative Growth in Rice

Low light intensity affects several physiological parameters during the different growth stages in rice. Plants have various regulatory mechanisms to cope with stresses. One of them is the differential and temporal expression of genes, which is governed by post-transcriptional gene expression regulation through endogenous miRNAs. To decipher low light stress-responsive miRNAs in rice, miRNA expression profiling was carried out using next-generation sequencing of low-light-tolerant (Swarnaprabha) and -sensitive (IR8) rice genotypes through Illumina sequencing. Swarnaprabha and IR8 were subjected to 25% low light treatment for one day, three days, and five days at the active tillering stage. More than 43 million raw reads and 9 million clean reads were identified in Swarnaprabha, while more than 41 million raw reads and 8.5 million clean reads were identified in IR8 after NGS. Importantly, 513 new miRNAs in rice were identified, whose targets were mostly regulated by the genes involved in photosynthesis and metabolic pathways. Additionally, 114 known miRNAs were also identified. Five novel (osa-novmiR1, osa-novmiR2, osa-novmiR3, osa-novmiR4, and osa-novmiR5) and three known (osa-miR166c-3p, osa-miR2102-3p, and osa-miR530-3p) miRNAs were selected for their expression validation through miRNA-specific qRT-PCR. The expression analyses of most of the predicted targets of corresponding miRNAs show negative regulation. Hence, miRNAs modulated the expression of genes providing tolerance/susceptibility to low light stress. This information might be useful in the improvement of crop productivity under low light stress.

Excess nitrogen responsive HvMADS27 transcription factor controls barley root architecture by regulating abscisic acid level

Nitrogen (N) is an important element for plant growth and development. Although several studies have examined plants' response to N deficiency, studies on plants' response to excess N, which is common in fertilizer-based agrosystems, are limited. Therefore, the aim of this study was to examine the response of barley to excess N conditions, specifically the root response. Additionally, genomic mechanism of excess N response in barley was elucidated using transcriptomic technologies. The results of the study showed that barley MADS27 transcription factor was mainly expressed in the roots and its gene contained N-responsive cis-regulatory elements in the promoter region. Additionally, there was a significant decrease in HvMADS27 expression under excess N condition; however, its expression was not significantly affected under low N condition. Phenotypic analysis of the root system of HvMADS27 knockdown and overexpressing barley plants revealed that HvMADS27 regulates barley root architecture under excess N stress. Further analysis of wild-type (WT) and transgenic barley plants (hvmads27 kd and hvmads27 c-Myc OE) revealed that HvMADS27 regulates the expression of HvBG1 β-glucosidase, which in turn regulates abscisic acid (ABA) level in roots. Overall, the findings of this study showed that HvMADS27 expression is downregulated in barley roots under excess N stress, which induces HvBG1 expression, leading to the release of ABA from ABA-glucose conjugate, and consequent shortening of the roots.

Integrated analysis of small RNAs, transcriptome and degradome sequencing reveal the drought stress network in Agropyron mongolicum Keng

Agropyron mongolicum (A. mongolicum) is an excellent gramineous forage with extreme drought tolerance, which lives in arid and semiarid desert areas. However, the mechanism that underlies the response of microRNAs (miRNAs) and their targets in A. mongolicum to drought stress is not well understood. In this study, we analyzed the transcriptome, small RNAome (specifically the miRNAome) and degradome to generate a comprehensive resource that focused on identifying key regulatory miRNA-target circuits under drought stress. The most extended transcript in each collection is known as the UniGene, and a total of 41,792 UniGenes and 1,104 miRNAs were identified, and 99 differentially expressed miRNAs negatively regulated 1,474 differentially expressed target genes. Among them, eight miRNAs were unique to A. mongolicum, and there were 36 target genes. A weighted gene co-expression network analysis identified five hub genes. The miRNAs of five hub genes were screened with an integration analysis of the degradome and sRNAs, such as osa-miR444a-3p.2-MADS47, bdi-miR408-5p_1ss19TA-CCX1, tae-miR9774_L-2R-1_1ss11GT-carC, ata-miR169a-3p-PAO2, and bdi-miR528-p3_2ss15TG20CA-HOX24. The functional annotations revealed that they were involved in mediating the brassinosteroid signal pathway, transporting and exchanging sodium and potassium ions and regulating the oxidation-reduction process, hydrolase activity, plant response to water deprivation, abscisic acid (ABA) and the ABA-activated signaling pathway to regulate drought stress. Five hub genes were discovered, which could play central roles in the regulation of drought-responsive genes. These results show that the combined analysis of miRNA, the transcriptome and degradation group provides a useful platform to investigate the molecular mechanism of drought resistance in A. mongolicum and could provide new insights into the genetic engineering of Poaceae crops in the future.

Rapid sequence and functional diversification of a miRNA superfamily targeting calcium signaling components in seed plants

MicroRNA (miRNA)-directed posttranscriptional gene silencing (miR-PTGS) is an integral component of gene regulatory networks governing plant development and responses to the environment. The sequence homology between Sly-miR4376, a miRNA common to Solanaceae and reported to target autoinhibited Ca²⁺ -ATPase 10 (ACA10) messenger RNA (mRNA) in tomato, and Arabidopsis miR391 (Ath-miR391), previously annotated as a nonconserved member of the deeply conserved miR390 family, has prompted us to revisit the function of Ath-miR391, as well as its regulatory conservation. A combination of genetic, molecular, and bioinformatic analyses revealed a hidden conservation for miR-PTGS of ACA10 homologs in spermatophytes. We found that the Arabidopsis ACA10 mRNA undergoes miR391-directed cleavage in vivo. Furthermore, transgenic overexpression of miR391 recapitulated the compact inflorescence (cif) phenotypes characteristic of ACA10 loss-of-function mutants, due to miR391-directed PTGS of ACA10. Significantly, comprehensive data mining revealed robust evidence for widespread PTGS of ACA10 homologs directed by a superfamily of related miRNAs sharing a conserved sequence core. Intriguingly, the ACA-targeting miRNAs in Poaceae also direct PTGS for calmodulin-like proteins which are putative Ca²⁺ sensors. The PTGS of ACA10 homologs is therefore directed by a miRNA superfamily that is of ancient origin and has undergone rapid sequence diversification associated with functional innovation.

Nitrate-dependent regulation of miR444-OsMADS27 signalling cascade controls root development in rice

Nitrate is an important nutrient and a key signalling molecule for plant development. A number of transcription factors involved in the response to nitrate and their regulatory mechanisms have been identified. However, little is known about the transcription factors involved in nitrate sensing and their regulatory mechanisms among crop plants. In this study, we identified functions of a nitrate-responsive miR444:MADS-box transcription factor OsMADS27 module and its downstream targets mediating rice root growth and stress responses. Transgenic rice plants expressing miR444 target mimic improved rice root growth. Although miR444 has the potential to target multiple genes, we identified OsMADS27 as the major miR444 target that regulates the expression of nitrate transporters, as well as several key genes including expansins, and those associated with auxin signalling, to promote root growth. In agreement with this, overexpression of miRNA-resistant OsMADS27 improved root development and tolerance to abiotic stresses, while its silencing suppressed root growth. OsMADS27 mediated robust stress tolerance in plants through its ability to bind to the promoters of specific stress regulators, as observed in ChIP-seq analysis. Our results provide evidence of a nitrate-dependent miR444-OsMADS27 signalling cascade involved in the regulation of rice root growth, as well as its surprising role in stress responses.

Small regulatory RNAs in rice epigenetic regulation

Plant small RNAs (sRNAs) are short non-coding RNAs that are implicated in various regulatory processes involving post-transcriptional gene silencing and epigenetic gene regulation. In epigenetic regulation, sRNAs are primarily involved in RNA-directed DNA methylation (RdDM) pathways. sRNAs in the RdDM pathways play a role not only in the suppression of transposable element (TE) activity but also in gene expression regulation. Although the major components of the RdDM pathways have been well studied in Arabidopsis, recent studies have revealed that the RdDM pathways in rice have important biological functions in stress response and developmental processes. In this review, we summarize and discuss recent literature on sRNA-mediated epigenetic regulation in rice. First, we describe the RdDM mechanisms in plants. We then introduce recent discoveries on the biological roles of rice genes involved in the RdDM pathway and TE-derived sRNAs working at specific genomic loci for epigenetic control in rice.

TRANS-ACTING SIRNA3-derived short interfering RNAs confer cleavage of mRNAs in rice

Plant TRANS-ACTING SIRNA3 (TAS3)-derived short interfering RNAs (siRNAs) include tasiR-AUXIN RESPONSE FACTORs (ARFs), which are functionally conserved in targeting ARF genes, and a set of non-tasiR-ARF siRNAs, which have rarely been studied. In this study, TAS3 siRNAs were systematically characterized in rice (Oryza sativa). Small RNA sequencing results showed that an overwhelming majority of TAS3 siRNAs belong to the non-tasiR-ARF group, while tasiR-ARFs occupy a diminutive fraction. Phylogenetic analysis of TAS3 genes across dicot and monocot plants revealed that the siRNA-generating regions were highly conserved in grass species, especially in the Oryzoideae. Target genes were identified for not only tasiR-ARFs but also non-tasiR-ARF siRNAs by analyzing rice Parallel Analysis of RNA Ends datasets, and some of these siRNA-target interactions were experimentally confirmed using tas3 mutants generated by genome editing. Consistent with the de-repression of target genes, phenotypic alterations were observed for mutants in three TAS3 loci in comparison to wild-type rice. The regulatory role of ribosomes in the TAS3 siRNA-target interactions was further revealed by the fact that TAS3 siRNA-mediated target cleavage, in particular tasiR-ARFs targeting ARF2/3/14/15, occurred extensively in rice polysome samples. Altogether, our study sheds light into TAS3 genes in plants and expands our knowledge about rice TAS3 siRNA-target interactions.

Identification and characterization of heat-responsive microRNAs at the booting stage in two rice varieties, 9311 and Nagina 22

MicroRNAs (miRNAs) are small, non-coding, regulatory RNAs that play important roles in abiotic stress responses in plants, but their regulatory roles in the adaptive response to heat stress at the booting stage in two rice varieties, 9311 and Nagina 22, remain largely unknown. In this study, 464 known miRNAs and 123 potential novel miRNAs were identified. Of these miRNAs, a total of 90 differentially expressed miRNAs were obtained with 9311 libraries as the control group, of which 54 were upregulated and 36 were downregulated. To gain insight into functional significance, 2773 potential target genes of these 90 differentially expressed miRNAs were predicted. GO enrichment analysis showed that the predicted target genes of differentially expressed miRNAs included NACs, LACs, CSD, and Hsp40. KEGG pathway analysis showed that the target genes of these differentially expressed miRNAs were significantly enriched in the plant hormone signal transduction pathway. The expression levels of 10 differentially expressed miRNAs and their target genes obtained by qRT-PCR were largely consistent with the sequencing results. This study lays a foundation for the elucidation of the miRNA-mediated regulatory mechanisms in rice at elevated temperatures.

Reproductive phasiRNA loci and DICER-LIKE5, but not microRNA loci, diversified in monocotyledonous plants

In monocots other than maize (Zea mays) and rice (Oryza sativa), the repertoire and diversity of microRNAs (miRNAs) and the populations of phased, secondary, small interfering RNAs (phasiRNAs) are poorly characterized. To remedy this, we sequenced small RNAs (sRNA) from vegetative and dissected inflorescence tissue in 28 phylogenetically diverse monocots and from several early-diverging angiosperm lineages, as well as publicly available data from 10 additional monocot species. We annotated miRNAs, small interfering RNAs (siRNAs) and phasiRNAs across the monocot phylogeny, identifying miRNAs apparently lost or gained in the grasses relative to other monocot families, as well as a number of transfer RNA fragments misannotated as miRNAs. Using our miRNA database cleaned of these misannotations, we identified conservation at the 8th, 9th, 19th, and 3'-end positions that we hypothesize are signatures of selection for processing, targeting, or Argonaute sorting. We show that 21-nucleotide (nt) reproductive phasiRNAs are far more numerous in grass genomes than other monocots. Based on sequenced monocot genomes and transcriptomes, DICER-LIKE5, important to 24-nt phasiRNA biogenesis, likely originated via gene duplication before the diversification of the grasses. This curated database of phylogenetically diverse monocot miRNAs, siRNAs, and phasiRNAs represents a large collection of data that should facilitate continued exploration of sRNA diversification in flowering plants.

Comprehensive analysis of the isomiRome in the vegetative organs of the conifer Pinus pinaster under contrasting water availability

An increasing number of microRNAs (miRNAs) and miRNA-related sequences produced during miRNA biogenesis, comprising the isomiRome, have been recently highlighted in different species as critical mediators of environmental stress responses. Conifers have some of the largest known genomes but an extensive characterization of the isomiRome from any conifer species has been lacking. We provide here a comprehensive overview of the Pinus pinaster isomiRome expressed in roots, stem and needles under well-watered and drought conditions. From the 13,441 unique small RNA sequences identified, 2,980 were annotated as canonical miRNAs or miRNA* and the remaining were classified as isomiRNA or miRNA-like sequences. A survey of their expression patterns highlighted roots as the most responsive organ under drought, where specific sequences of which a 24-nt novel miRNA stood out, were strongly down-regulated. Given the putative roles of the miRNA-targeted transcripts validated specifically in root tissues, some of the miRNAs, conserved and novel, are shortlisted as potential regulators of drought response. These results provide a valuable resource for comparative studies between gymnosperms and angiosperms. Furthermore, it evidences high transferability of the isomiRome between pine species being a useful basis for further molecular regulation and physiological studies, and especially those focused on adaptation to drought conditions.

Development and use of miRNA-derived SSR markers for the study of genetic diversity, population structure, and characterization of genotypes for breeding heat tolerant wheat varieties

Heat stress is an important abiotic factor that limits wheat production globally, including south-east Asia. The importance of micro (mi) RNAs in gene expression under various biotic and abiotic stresses is well documented. Molecular markers, specifically simple sequence repeats (SSRs), play an important role in the wheat improvement breeding programs. Given the role of miRNAs in heat stress-induced transcriptional regulation and acclimatization, the development of miRNA-derived SSRs would prove useful in studying the allelic diversity at the heat-responsive miRNA-genes in wheat. In the present study, efforts have been made to identify SSRs from 96 wheat heat-responsive miRNA-genes and their characterization using a panel of wheat genotypes with contrasting reactions (tolerance/susceptible) to heat stress. A set of 13 miRNA-derived SSR markers were successfully developed as an outcome. These miRNA-SSRs are located on 11 different common wheat chromosomes (2A, 3A, 3B, 3D, 4D, 5A, 5B, 5D, 6A, 6D, and 7A). Among 13 miRNA-SSRs, seven were polymorphic on a set of 37 selected wheat genotypes. Within these polymorphic SSRs, three makers, namely HT-169j, HT-160a, and HT-160b, were found promising as they could discriminate heat-tolerant and heat-susceptible genotypes. This is the first report of miRNA-SSR development in wheat and their deployment in genetic diversity and population structure studies and characterization of trait-specific germplasm. The study suggests that this new class of molecular makers has great potential in the marker-assisted breeding (MAB) programs targeted at improving heat tolerance and other adaptability or developmental traits in wheat and other crops.

Northern Blotting Technique for Detection and Expression Analysis of mRNAs and Small RNAs

Northern analysis is a conventional but gold standard method for detection and quantification of gene expression changes. It not only detects the presence of a transcript but also indicates size and relative comparison of transcript abundance on a single membrane. In recent years it has been aptly adapted to validate and study the size and expression of small noncoding RNAs. Here, we describe protocols employed in our laboratory for conventional northern analysis with total RNA/mRNA to study gene expression and validation of small noncoding RNAs using low molecular weight fraction of RNAs. A brief account on the recent advancements for improving the sensitivity and efficiency of northern blot detection is also included in this chapter.

RNAi Mediated Hypoxia Stress Tolerance in Plants

Small RNAs regulate various biological process involved in genome stability, development, and adaptive responses to biotic or abiotic stresses. Small RNAs include microRNAs (miRNAs) and small interfering RNAs (siRNAs). MicroRNAs (miRNAs) are regulators of gene expression that affect the transcriptional and post-transcriptional regulation in plants and animals through RNA interference (RNAi). miRNAs are endogenous small RNAs that originate from the processing of non-coding primary miRNA transcripts folding into hairpin-like structures. The mature miRNAs are incorporated into the RNA-induced silencing complex (RISC) and drive the Argonaute (AGO) proteins towards their mRNA targets. siRNAs are generated from a double-stranded RNA (dsRNA) of cellular or exogenous origin. siRNAs are also involved in the adaptive response to biotic or abiotic stresses. The response of plants to hypoxia includes a genome-wide transcription reprogramming. However, little is known about the involvement of RNA signaling in gene regulation under low oxygen availability. Interestingly, miRNAs have been shown to play a role in the responses to hypoxia in animals, and recent evidence suggests that hypoxia modulates the expression of various miRNAs in plant systems. In this review, we describe recent discoveries on the impact of RNAi on plant responses to hypoxic stress in plants.

Spatiotemporal Resolved Leaf Angle Establishment Improves Rice Grain Yield via Controlling Population Density

Leaf angle is mainly determined by the lamina joint (LJ) and contributes to ideal crop architecture for high yield. Here, we dissected five successive stages with distinct cytological features of LJs spanning organogenesis to leaf angle formation and obtained the underlying stage-specific mRNAs and small RNAs, which well explained the cytological dynamics during LJ organogenesis and leaf angle plasticity. Combining the gene coexpression correlation with high-throughput promoter analysis, we identified a set of transcription factors (TFs) determining the stage- and/or cytological structure-specific profiles. The functional studies of these TFs demonstrated that cytological dynamics determined leaf angle and that the knockout rice of these TFs with erect leaves significantly enhanced yield by maintaining the proper tiller number under dense planting. This work revealed the high-resolution mechanisms of how the cytological dynamics of LJ determined leaf erectness and served as a valuable resource to remodel rice architecture for high yield by controlling population density.

Genome-Wide Identification of RNA Silencing-Related Genes and Their Expressional Analysis in Response to Heat Stress in Barley (Hordeum vulgare L.)

Barley (Hordeum vulgare L.) is an economically important crop cultivated in temperate climates all over the world. Adverse environmental factors negatively affect its survival and productivity. RNA silencing is a conserved pathway involved in the regulation of growth, development and stress responses. The key components of RNA silencing are the Dicer-like proteins (DCLs), Argonautes (AGOs) and RNA-dependent RNA polymerases (RDRs). Despite its economic importance, there is no available comprehensive report on barley RNA silencing machinery and its regulation. In this study, we in silico identified five DCL (HvDCL), eleven AGO (HvAGO) and seven RDR (HvRDR) genes in the barley genome. Genomic localization, phylogenetic analysis, domain organization and functional/catalytic motif identification were also performed. To understand the regulation of RNA silencing, we experimentally analysed the transcriptional changes in response to moderate, persistent or gradient heat stress treatments: transcriptional accumulation of siRNA- but not miRNA-based silencing factor was consistently detected. These results suggest that RNA silencing is dynamically regulated and may be involved in the coordination of development and environmental adaptation in barley. In summary, our work provides information about barley RNA silencing components and will be a ground for the selection of candidate factors and in-depth functional/mechanistic analyses.

Promotion of BR Biosynthesis by miR444 Is Required for Ammonium-Triggered Inhibition of Root Growth

Rice (Oryza sativa), the staple food for almost half of the world's population, prefers ammonium (NH₄ ⁺) as the major nitrogen resource, and while NH₄ ⁺ has profound effects on rice growth and yields, the underlying regulatory mechanisms remain largely unknown. Brassinosteroids (BRs) are a class of steroidal hormones playing key roles in plant growth and development. In this study, we show that NH₄ ⁺ promotes BR biosynthesis through miR444 to regulate rice root growth. miR444 targeted five homologous MADS-box transcription repressors potentially forming homologous or heterogeneous complexes in rice. miR444 positively regulated BR biosynthesis through its MADS-box targets, which directly repress the transcription of BR-deficient dwarf 1 (OsBRD1), a key BR biosynthetic gene. NH₄ ⁺ induced the miR444-OsBRD1 signaling cascade in roots, thereby increasing the amount of BRs, whose biosynthesis and signaling were required for NH₄ ⁺ -dependent root elongation inhibition. Consistently, miR444-overexpressing rice roots were hypersensitive to NH₄ ⁺ depending on BR biosynthesis, and overexpression of miR444's target, OsMADS57, resulted in rice hyposensitivity to NH₄ ⁺ in root elongation, which was associated with a reduction of BR content. In summary, our findings reveal a cross talk mechanism between NH₄ ⁺ and BR in which NH₄ ⁺ activates miR444-OsBRD1, an undescribed BR biosynthesis-promoting signaling cascade, to increase BR content, inhibiting root elongation in rice.

Analysis of miRNAs in Two Wheat Cultivars Infected With Puccinia striiformis f. sp. tritici

MicroRNAs are small RNAs that regulate gene expression in eukaryotes. In this study, we analyzed the small RNA profiles of two cultivars that exhibit different reactions to stripe rust infection: one susceptible, the other partially resistant. Using small RNA libraries prepared from the two wheat cultivars infected with stripe rust fungus (Puccinia striiformis f. sp. tritici), we identified 182 previously known miRNAs, 91 variants of known miRNAs, and 163 candidate novel wheat miRNAs. Known miRNA loci were usually copied in all three wheat sub-genomes, whereas novel miRNA loci were often specific to a single sub-genome. DESeq2 analysis of differentially expressed microRNAs revealed 23 miRNAs that exhibit cultivar-specific differences. TA078/miR399b showed cultivar-specific differential regulation in response to infection. Using different target prediction algorithms, 145 miRNAs were predicted to target wheat genes, while 69 miRNAs were predicted to target fungal genes. We also confirmed reciprocal expression of TA078/miR399b and tae-miR9664 and their target genes in different treatments, providing evidence for miRNA-mediated regulation during infection. Both known and novel miRNAs were predicted to target fungal genes, suggesting trans-kingdom regulation of gene expression. Overall, this study contributes to the current repository of wheat miRNAs and provides novel information on the yet-uncharacterized roles for miRNAs in the wheat-stripe rust pathosystem.

Independent origin of MIRNA genes controlling homologous target genes by partial inverted duplication of antisense-transcribed sequences

Some microRNAs (miRNAs) are key regulators of developmental processes, mainly by controlling the accumulation of transcripts encoding transcription factors that are important for morphogenesis. MADS-box genes encode a family of transcription factors which control diverse developmental processes in flowering plants. Here we study the convergent evolution of two MIRNA (MIR) gene families, named MIR444 and MIR824, targeting members of the same clade of MIKC^C -group MADS-box genes. We show that these two MIR genes most likely originated independently in monocots (MIR444) and in Brassicales (eudicots, MIR824). We provide evidence that, in both cases, the future target gene was transcribed in antisense prior to the evolution of the MIR genes. Both MIR genes then likely originated by a partial inverted duplication of their target genes, resulting in natural antisense organization of the newly evolved MIR gene and its target gene at birth. We thus propose a model for the origin of MIR genes, MEPIDAS (MicroRNA Evolution by Partial Inverted Duplication of Antisense-transcribed Sequences). MEPIDAS is a refinement of the inverted duplication hypothesis. According to MEPIDAS, a MIR gene evolves at a genomic locus at which the future target gene is also transcribed in the antisense direction. A partial inverted duplication at this locus causes the antisense transcript to fold into a stem-loop structure that is recognized by the miRNA biogenesis machinery to produce a miRNA that regulates the gene at this locus. Our analyses exemplify how to elucidate the origin of conserved miRNAs by comparative genomics and will guide future studies. OPEN RESEARCH BADGE: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/genbank/.

SMA1, a homolog of the splicing factor Prp28, has a multifaceted role in miRNA biogenesis in Arabidopsis

MicroRNAs (miRNAs) are a class of small non-coding RNAs that repress gene expression. In plants, the RNase III enzyme Dicer-like (DCL1) processes primary miRNAs (pri-miRNAs) into miRNAs. Here, we show that SMALL1 (SMA1), a homolog of the DEAD-box pre-mRNA splicing factor Prp28, plays essential roles in miRNA biogenesis in Arabidopsis. A hypomorphic sma1-1 mutation causes growth defects and reduces miRNA accumulation correlated with increased target transcript levels. SMA1 interacts with the DCL1 complex and positively influences pri-miRNA processing. Moreover, SMA1 binds the promoter region of genes encoding pri-miRNAs (MIRs) and is required for MIR transcription. Furthermore, SMA1 also enhances the abundance of the DCL1 protein levels through promoting the splicing of the DCL1 pre-mRNAs. Collectively, our data provide new insights into the function of SMA1/Prp28 in regulating miRNA abundance in plants.

MicroRNA Gene Regulation in Extremely Young and Parallel Adaptive Radiations of Crater Lake Cichlid Fish

Cichlid fishes provide textbook examples of explosive phenotypic diversification and sympatric speciation, thereby making them ideal systems for studying the molecular mechanisms underlying rapid lineage divergence. Despite the fact that gene regulation provides a critical link between diversification in gene function and speciation, many genomic regulatory mechanisms such as microRNAs (miRNAs) have received little attention in these rapidly diversifying groups. Therefore, we investigated the posttranscriptional regulatory role of miRNAs in the repeated sympatric divergence of Midas cichlids (Amphilophus spp.) from Nicaraguan crater lakes. Using miRNA and mRNA sequencing of embryos from five Midas species, we first identified miRNA binding sites in mRNAs and highlighted the presences of a surprising number of novel miRNAs in these adaptively radiating species. Then, through analyses of expression levels, we identified putative miRNA/gene target pairs with negatively correlated expression level that were consistent with the role of miRNA in downregulating mRNA. Furthermore, we determined that several miRNA/gene pairs show convergent expression patterns associated with the repeated benthic/limnetic sympatric species divergence implicating these miRNAs as potential molecular mechanisms underlying replicated sympatric divergence. Finally, as these candidate miRNA/gene pairs may play a central role in phenotypic diversification in these cichlids, we characterized the expression domains of selected miRNAs and their target genes via in situ hybridization, providing further evidence that miRNA regulation likely plays a role in the Midas cichlid adaptive radiation. These results provide support for the hypothesis that extremely quickly evolving miRNA regulation can contribute to rapid evolutionary divergence even in the presence of gene flow.

Small RNA discovery in the interaction between barley and the powdery mildew pathogen

Background

Plants encounter pathogenic and non-pathogenic microorganisms on a nearly constant basis. Small RNAs such as siRNAs and miRNAs/milRNAs influence pathogen virulence and host defense responses. We exploited the biotrophic interaction between the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh), and its diploid host plant, barley (Hordeum vulgare) to explore fungal and plant sRNAs expressed during Bgh infection of barley leaf epidermal cells.

Results

RNA was isolated from four fast-neutron immune-signaling mutants and their progenitor over a time course representing key stages of Bgh infection, including appressorium formation, penetration of epidermal cells, and development of haustorial feeding structures. The Cereal Introduction (CI) 16151 progenitor carries the resistance allele Mla6, while Bgh isolate 5874 harbors the AVR_a6 avirulence effector, resulting in an incompatible interaction. Parallel Analysis of RNA Ends (PARE) was used to verify sRNAs with likely transcript targets in both barley and Bgh. Bgh sRNAs are predicted to regulate effectors, metabolic genes, and translation-related genes. Barley sRNAs are predicted to influence the accumulation of transcripts that encode auxin response factors, NAC transcription factors, homeodomain transcription factors, and several splicing factors. We also identified phasing small interfering RNAs (phasiRNAs) in barley that overlap transcripts that encode receptor-like kinases (RLKs) and nucleotide-binding, leucine-rich domain proteins (NLRs).

Conclusions

These data suggest that Bgh sRNAs regulate gene expression in metabolism, translation-related, and pathogen effectors. PARE-validated targets of predicted Bgh milRNAs include both EKA (effectors homologous to AVR_k1 and AVR_a10) and CSEP (candidate secreted effector protein) families. We also identified barley phasiRNAs and miRNAs in response to Bgh infection. These include phasiRNA loci that overlap with a significant proportion of receptor-like kinases, suggesting an additional sRNA control mechanism may be active in barley leaves as opposed to predominant R-gene phasiRNA overlap in many eudicots. In addition, we identified conserved miRNAs, novel miRNA candidates, and barley genome mapped sRNAs that have PARE validated transcript targets in barley. The miRNA target transcripts are enriched in transcription factors, signaling-related proteins, and photosynthesis-related proteins. Together these results suggest both barley and Bgh control metabolism and infection-related responses via the specific accumulation and targeting of genes via sRNAs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5947-z) contains supplementary material, which is available to authorized users.

Transcriptomic analyses of rice (Oryza sativa) genes and non-coding RNAs under nitrogen starvation using multiple omics technologies

BACKGROUND

Nitrogen (N) is a key macronutrient essential for plant growth, and its availability has a strong influence on crop development. The application of synthetic N fertilizers on crops has increased substantially in recent decades; however, the applied N is not fully utilized due to the low N use efficiency of crops. To overcome this limitation, it is important to understand the genome-wide responses and functions of key genes and potential regulatory factors in N metabolism.

RESULTS

Here, we characterized changes in the rice (Oryza sativa) transcriptome, including genes, newly identified putative long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) and their target mRNAs in response to N starvation using four different transcriptome approaches. Analysis of rice genes involved in N metabolism and/or transport using strand-specific RNA-Seq identified 2588 novel putative lncRNA encoding loci. Analysis of previously published RNA-Seq datasets revealed a group of N starvation-responsive lncRNAs showing differential expression under other abiotic stress conditions. Poly A-primed sequencing (2P-Seq) revealed alternatively polyadenylated isoforms of N starvation-responsive lncRNAs and provided precise 3' end information on the transcript models of these lncRNAs. Analysis of small RNA-Seq data identified N starvation-responsive miRNAs and down-regulation of miR169 family members, causing de-repression of NF-YA, as confirmed by strand-specific RNA-Seq and qRT-PCR. Moreover, we profiled the N starvation-responsive down-regulation of root-specific miRNA, osa-miR444a.4-3p, and Degradome sequencing confirmed MADS25 as a novel target gene.

CONCLUSIONS

In this study, we used a combination of multiple RNA-Seq analyses to extensively profile the expression of genes, newly identified lncRNAs, and microRNAs in N-starved rice roots and shoots. Data generated in this study provide an in-depth understanding of the regulatory pathways modulated by N starvation-responsive miRNAs. The results of comprehensive, large-scale data analysis provide valuable information on multiple aspects of the rice transcriptome, which may be useful in understanding the responses of rice plants to changes in the N supply status of soil.

Coupling of microRNA-directed phased small interfering RNA generation from long noncoding genes with alternative splicing and alternative polyadenylation in small RNA-mediated gene silencing

MicroRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs) play vital regulatory roles in plant growth and development. Little is known about these small RNAs in litchi (Litchi chinensis), an economically important fruit crop widely cultivated in Southeast Asia. We profiled the litchi small RNA population with various deep-sequencing techniques and in-depth bioinformatic analyses. The genome-wide identification of miRNAs, their target genes, and phasiRNA-generating (PHAS) genes/loci showed that the function of miR482/2118 has expanded, relative to its canonical function. We also discovered that, for 29 PHAS loci, miRNA-mediated phasiRNA production was coupled with alternative splicing (AS) and alternative polyadenylation (APA). Most of these loci encoded long noncoding RNAs. An miR482/2118 targeted locus gave rise to four main transcript isoforms through AS/APA, and diverse phasiRNAs generated from these isoforms appeared to target long terminal repeat (LTR) retrotransposons and other unrelated genes. This coupling enables phasiRNA production from different exons of noncoding PHAS genes and yields diverse phasiRNA populations, both broadening and altering the range of downstream phasiRNA-regulated genes. Our results reveal the diversity of miRNA and phasiRNA in litchi, and demonstrate AS/APA as a new layer of regulation in small RNA-mediated gene silencing.

Genome-wide analysis of rice cis-natural antisense transcription under cadmium exposure using strand-specific RNA-Seq

BACKGROUND

The elucidation of novel transcripts and their expression in response to various stress conditions is necessary to understand the transcriptional network of plants as an adaptation to biotic and abiotic stresses. We performed strand-specific RNA-Seq (ssRNA-Seq) on rice exposed to cadmium (Cd) for 24 h and investigated the expression of cis-natural antisense transcripts (cis-NATs), a class of endogenous coding or non-protein-coding RNAs with sequence complementarity to the opposite strands of RAP transcripts.

RESULTS

Many RAP transcripts possessed cis-NATs and these cis-NATs were responsive to some extent. Cis-NATs were upregulated from 26, 266 and 409 RAP gene loci, while 2054, 2501 and 2825 RAP transcripts were upregulated from 38,123 RAP loci under high Cd exposure in roots at 1, 12 and 24 h, respectively. In addition, most of the upregulated cis-NATs showed little upregulation under ABA or cold treatment. A number of cis-NATs were upregulated from less than 35 RAP gene loci in different tissue and time-point combinations under low Cd exposure, suggesting that cis-NATs respond to environmental stress. Furthermore, 409 RAP transcripts with upregulated cis-NATs were classified into three groups based on the expression of the RAP transcripts from the opposite DNA strand, including 138 upregulated, 128 invariable, and 143 downregulated transcripts, although the responses of cis-NATs and RAP transcripts were not always correlated.

CONCLUSIONS

We have shown that the cis-NATs identified by ssRNA-Seq analysis are novel genes and that some of them are stress-specific and show different responses depending on the degree of stress and tissue. These results improve our understanding of the complete molecular mechanism of plant adaptation to Cd exposure.

Response of microRNAs to cold treatment in the young spikes of common wheat

Background

MicroRNAs (miRNAs) are a class of small non-coding RNAs that play important roles in biotic and abiotic stresses by regulating their target genes. For common wheat, spring frost damage frequently occurs, especially when low temperature coincides with plants at early floral organ differentiation, which may result in significant yield loss. Up to date, the role of miRNAs in wheat response to frost stress is not well understood.

Results

We report here the sequencing of small RNA transcriptomes from the young spikes that were treated with cold stress and the comparative analysis with those of the control. A total of 192 conserved miRNAs from 105 families and nine novel miRNAs were identified. Among them, 34 conserved and five novel miRNAs were differentially expressed between the cold-stressed samples and the controls. The expression patterns of 18 miRNAs were further validated by quantitative real time polymerase chain reaction (qRT-PCR). Moreover, nearly half of the miRNAs were cross inducible by biotic and abiotic stresses when compared with previously published work. Target genes were predicted and validated by degradome sequencing. Gene Ontology (GO) enrichment analysis showed that the target genes of differentially expressed miRNAs were enriched for response to the stimulus, regulation of transcription, and ion transport functions. Since many targets of differentially expressed miRNAs were transcription factors that are associated with floral development such as ARF, SPB (Squamosa Promoter Binding like protein), MADS-box (MCM1, AG, DEFA and SRF), MYB, SPX (SYG1, Pho81 and XPR1), TCP (TEOSINTE BRANCHED, Cycloidea and PCF), and PPR (PentatricoPeptide Repeat) genes, cold-altered miRNA expression may cause abnormal reproductive organ development.

Conclusion

Analysis of small RNA transcriptomes and their target genes provide new insight into miRNA regulation in developing wheat inflorescences under cold stress. MiRNAs provide another layer of gene regulation in cold stress response that can be genetically manipulated to reduce yield loss in wheat.

Electronic supplementary material

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

Evolutionary plasticity of restorer-of-fertility-like proteins in rice

Hybrid seed production in rice relies on cytoplasmic male sterility (CMS) induced by specific mitochondrial proteins, whose deleterious effects are suppressed by nuclear Restorer of Fertility (RF) genes. The majority of RF proteins belong to a specific clade of the RNA-binding pentatricopeptide repeat protein family. We have characterised 'restorer-of-fertility-like' (RFL) sequences from 13 Oryza genomes and the Brachypodium distachyon genome. The majority of the RFL sequences are found in genomic clusters located at two or three chromosomal loci with only a minor proportion being present as isolated genes. The RFL genomic cluster located on Oryza chromosome 10, the location of almost all known active rice RF genes, shows extreme variation in structure and gene content between species. We show evidence for homologous recombination events as an efficient mechanism for generating the huge repertoire of RNA sequence recognition motifs within RFL proteins and a major driver of RFL sequence evolution. The RFL sequences identified here will improve our understanding of the molecular basis of CMS and fertility restoration in plants and will accelerate the development of new breeding strategies.

Epigenetic regulation and epigenomic landscape in rice

Epigenetic regulation has been implicated in the control of complex agronomic traits in rice (Oryza sativa), a staple food crop and model monocot plant. Recent advances in high-throughput sequencing and the moderately complex genome of rice have made it possible to study epigenetic regulation in rice on a genome-wide scale. This review discusses recent advances in our understanding of epigenetic regulation in rice, with an emphasis on the roles of key epigenetic regulators, the epigenomic landscape, epigenetic variation, transposon repression, and plant development.

Generation of different sizes and classes of small RNAs in barley is locus, chromosome and/or cultivar-dependent

Background

Various small RNA (sRNA) sizes and varieties have been identified, but their relationship as well as relationship with their origins and allocations have not been well understood or investigated.

Results

By comparing sRNAs generated from two barley cultivars, Golden Promise (GP) and Pallas, we identified that the generation of different sizes and types of sRNAs in barley was locus-, chromosome- and/or cultivar-dependent. 20-nt sRNAs mainly comprising miRNAs and chloroplast-derived sRNAs were significantly over-expressed in Pallas vs. GP on chromosomes 3H and 6H. MiRNAs-enriched 21-nt sRNAs were significantly over-expressed in Pallas vs. GP only on chromosome 4H. On chromosome 5H this size of sRNAs was significantly under-expressed in Pallas, so were 22-nt sRNAs mainly comprising miRNAs and repeat-derived sRNAs. 24-nt sRNAs mostly derived from repeats were evenly distributed in all chromosomes and expressed similarly between GP and Pallas. Unlike other sizes of sRNAs, 24-nt sRNAs were little conserved in other plant species. Abundant sRNAs were mostly generated from 3’ terminal regions of chromosome 1H and 5’ terminal regions of chromosome 5H. Over-expressed miRNAs in GP vs. Pallas primarily function in stress responses and iron-binding.

Conclusions

Our study indicates that 23−24-nt sRNAs may be linked to repressive chromatin modifications and function in genome stability while 20−21-nt sRNAs may be important for the cultivar specificity. This study provides a novel insight into the mechanism of sRNA expression and function in barley.

Electronic supplementary material

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

Comprehensive analysis of microRNA-Seq and target mRNAs of rice sheath blight pathogen provides new insights into pathogenic regulatory mechanisms

MicroRNAs (miRNAs) are ∼22 nucleotide non-coding RNAs that regulate gene expression by targeting mRNAs for degradation or inhibiting protein translation. To investigate whether miRNAs regulate the pathogenesis in necrotrophic fungus Rhizoctonia solani AG1 IA, which causes significant yield loss in main economically important crops, and to determine the regulatory mechanism occurring during pathogenesis, we constructed hyphal small RNA libraries from six different infection periods of the rice leaf. Through sequencing and analysis, 177 miRNA-like small RNAs (milRNAs) were identified, including 15 candidate pathogenic novel milRNAs predicted by functional annotations of their target mRNAs and expression patterns of milRNAs and mRNAs during infection. Reverse transcription-quantitative polymerase chain reaction results for randomly selected milRNAs demonstrated that our novel comprehensive predictions had a high level of accuracy. In our predicted pathogenic protein-protein interaction network of R. solani, we added the related regulatory milRNAs of these core coding genes into the network, and could understand the relationships among these regulatory factors more clearly at the systems level. Furthermore, the putative pathogenic Rhi-milR-16, which negatively regulates target gene expression, was experimentally validated to have regulatory functions by a dual-luciferase reporter assay. Additionally, 23 candidate rice miRNAs that may involve in plant immunity against R. solani were discovered. This first study on novel pathogenic milRNAs of R. solani AG1 IA and the recognition of target genes involved in pathogenicity, as well as rice miRNAs, participated in defence against R. solani could provide new insights into revealing the pathogenic mechanisms of the severe rice sheath blight disease.

Genome-wide identification and characterization of novel lncRNAs in Populus under nitrogen deficiency

Long non-coding RNAs (lncRNAs) have been identified as important regulatory factors of gene expression in eukaryotic species, such as Homo sapiens, Arabidopsis thaliana, and Oryza sativa. However, the systematic identification of potential lncRNAs in trees is comparatively rare. In particular, the characteristics, expression, and regulatory roles of lncRNAs in trees under nutrient stress remain largely unknown. A genome-wide strategy was used in this investigation to identify and characterize novel and low-nitrogen (N)-responsive lncRNAs in Populus tomentosa; 388 unique lncRNA candidates belonging to 380 gene loci were detected and only seven lncRNAs were found to belong to seven conserved non-coding RNA families indicating the majority of P. tomentosa lncRNAs are species-specific. In total, 126 lncRNAs were significantly altered under low-N stress; 8 were repressed, and 118 were induced. Furthermore, 9 and 5 lncRNAs were detected as precursors of 11 known and 14 novel Populus miRNAs, respectively, whereas 4 lncRNAs were targeted by 29 miRNAs belonging to 5 families, including 22 conserved and 7 non-conserved miRNAs. In addition, 15 antisense lncRNAs were identified to be generated from opposite strands of 14 corresponding protein-coding genes. In total, 111 protein-coding genes with regions complementary to 38 lncRNAs were also predicted with some lncRNAs corresponding to multiple genes and vice versa, and their functions were annotated, which further demonstrated the complex regulatory relationship between lncRNAs and protein-coding genes in plants. Moreover, an interaction network among lncRNAs, miRNAs, and mRNAs was investigated. These findings enrich our understanding of lncRNAs in Populus, expand the methods of miRNA identification. Our results present the first global characterization of lncRNAs and their potential target genes in response to nitrogen stress in trees, which provides more information on low-nutrition adaptation mechanisms in woody plants.

A Signaling Cascade from miR444 to RDR1 in Rice Antiviral RNA Silencing Pathway

Plant RNA-DEPENDENT RNA POLYMERASE1 (RDR1) is a key component of the antiviral RNA-silencing pathway, contributing to the biogenesis of virus-derived small interfering RNAs. This enzyme also is responsible for producing virus-activated endogenous small interfering RNAs to stimulate the broad-spectrum antiviral activity through silencing host genes. The expression of RDR1 orthologs in various plants is usually induced by virus infection. However, the molecular mechanisms of activation of RDR1 expression in response to virus infection remain unknown. Here, we show that a monocot-specific microRNA, miR444, is a key factor in relaying the antiviral signaling from virus infection to OsRDR1 expression. The expression of miR444 is enhanced by infection with Rice stripe virus (RSV), and overexpression of miR444 improves rice (Oryza sativa) resistance against RSV infection accompanied by the up-regulation of OsRDR1 expression. We further show that three miR444 targets, the MIKC(C)-type MADS box proteins OsMADS23, OsMADS27a, and OsMADS57, form homodimers and heterodimers between them to repress the expression of OsRDR1 by directly binding to the CArG motifs of its promoter. Consequently, an increased level of miR444 diminishes the repressive roles of OsMADS23, OsMADS27a, and OsMADS57 on OsRDR1 transcription, thus activating the OsRDR1-dependent antiviral RNA-silencing pathway. We also show that overexpression of miR444-resistant OsMADS57 reduced OsRDR1 expression and rice resistance against RSV infection, and knockout of OsRDR1 reduced rice resistance against RSV infection. In conclusion, our results reveal a molecular cascade in the rice antiviral pathway in which miR444 and its MADS box targets directly control OsRDR1 transcription.

NATpipe: an integrative pipeline for systematical discovery of natural antisense transcripts (NATs) and phase-distributed nat-siRNAs from de novo assembled transcriptomes

Nat-siRNAs (small interfering RNAs originated from natural antisense transcripts) are a class of functional small RNA (sRNA) species discovered in both plants and animals. These siRNAs are highly enriched within the annealed regions of the NAT (natural antisense transcript) pairs. To date, great research efforts have been taken for systematical identification of the NATs in various organisms. However, developing a freely available and easy-to-use program for NAT prediction is strongly demanded by researchers. Here, we proposed an integrative pipeline named NATpipe for systematical discovery of NATs from de novo assembled transcriptomes. By utilizing sRNA sequencing data, the pipeline also allowed users to search for phase-distributed nat-siRNAs within the perfectly annealed regions of the NAT pairs. Additionally, more reliable nat-siRNA loci could be identified based on degradome sequencing data. A case study on the non-model plant Dendrobium officinale was performed to illustrate the utility of NATpipe. Finally, we hope that NATpipe would be a useful tool for NAT prediction, nat-siRNA discovery, and related functional studies. NATpipe is available at www.bioinfolab.cn/NATpipe/NATpipe.zip.

Identification of miRNAs and Their Targets in the Liverwort Marchantia polymorpha by Integrating RNA-Seq and Degradome Analyses

Bryophytes (liverworts, hornworts and mosses) comprise the three earliest diverging lineages of land plants (embryophytes). Marchantia polymorpha, a complex thalloid Marchantiopsida liverwort that has been developed into a model genetic system, occupies a key phylogenetic position. Therefore, M. polymorpha is useful in studies aiming to elucidate the evolution of gene regulation mechanisms in plants. In this study, we used computational, transcriptomic, small RNA and degradome analyses to characterize microRNA (miRNA)-mediated pathways of gene regulation in M. polymorpha. The data have been integrated into the open access ContigViews-miRNA platform for further reference. In addition to core components of the miRNA pathway, 129 unique miRNA sequences, 11 of which could be classified into seven miRNA families that are conserved in embryophytes (miR166a, miR390, miR529c, miR171-3p, miR408a, miR160 and miR319a), were identified. A combination of computational and degradome analyses allowed us to identify and experimentally validate 249 targets. In some cases, the target genes are orthologous to those of other embryophytes, but in other cases, the conserved miRNAs target either paralogs or members of different gene families. In addition, the newly discovered Mpo-miR11707.1 and Mpo-miR11707.2 are generated from a common precursor and target MpARGONAUTE1 (LW1759). Two other newly discovered miRNAs, Mpo-miR11687.1 and Mpo-miR11681.1, target the MADS-box transcription factors MpMADS1 and MpMADS2, respectively. Interestingly, one of the pentatricopeptide repeat (PPR) gene family members, MpPPR_66 (LW9825), the protein products of which are generally involved in various steps of RNA metabolism, has a long stem-loop transcript that can generate Mpo-miR11692.1 to autoregulate MpPPR_66 (LW9825) mRNA. This study provides a foundation for further investigations of the RNA-mediated silencing mechanism in M. polymorpha as well as of the evolution of this gene silencing pathway in embryophytes.

Systematic study of cis-antisense miRNAs in animal species reveals miR-3661 to target PPP2CA in human cells

MicroRNAs (miRNAs) suppress targeting gene expression through blocking translation or triggering mRNA degradation and, in general, act in trans, through a partially complementary interaction with the 3' untranslated region (3' UTR) or coding regions of a target gene. Although it has been reported previously that some miRNAs suppress their target genes on the opposite strand with a fully complementary sequence (i.e., natural antisense miRNAs that act in cis), there is no report to systematically study such cis-antisense miRNAs in different animal species. Here we report that cis-antisense miRNAs do exist in different animal species: 48 in Caenorhabditis elegans, 17 in Drosophila, 36 in Mus musculus, and 52 in Homo sapiens using a systematical bioinformatics approach. We show that most of these cis-antisense miRNAs can efficiently reduce the expression levels of their target genes in human cells. We further investigate hsa-miR-3661, one of the predicted cis-antisense miRNAs, in detail and demonstrate that this miRNA directly targets the coding sequence of PPP2CA located on the opposite DNA strand and inhibits the PPP2CA expression. Taken together, these results indicate that cis-antisense miRNAs are conservative and functional in animal species including humans.

Identification of Novel and Conserved miRNAs from Extreme Halophyte, Oryza coarctata, a Wild Relative of Rice

Oryza coarctata, a halophyte and wild relative of rice, is grown normally in saline water. MicroRNAs (miRNAs) are non-coding RNAs that play pivotal roles in every domain of life including stress response. There are very few reports on the discovery of salt-responsive miRNAs from halophytes. In this study, two small RNA libraries, one each from the control and salt-treated (450 mM NaCl for 24 h) leaves of O. coarctata were sequenced, which yielded 338 known and 95 novel miRNAs. Additionally, we used publicly available transcriptomics data of O. coarctata which led to the discovery of additional 48 conserved miRNAs along with their pre-miRNA sequences through in silico analysis. In total, 36 known and 7 novel miRNAs were up-regulated whereas, 12 known and 7 novel miRNAs were down-regulated under salinity stress. Further, 233 and 154 target genes were predicted for 48 known and 14 novel differentially regulated miRNAs respectively. These targets with the help of gene ontology analysis were found to be involved in several important biological processes that could be involved in salinity tolerance. Relative expression trends of majority of the miRNAs as detected by real time-PCR as well as predicted by Illumina sequencing were found to be coherent. Additionally, expression of most of the target genes was negatively correlated with their corresponding miRNAs. Thus, the present study provides an account of miRNA-target networking that is involved in salinity adaption of O. coarctata.

The mirror RNA expression pattern in human tissues

It has been realized in recent years that non-coding RNAs are playing important roles in genome functions and human diseases. Here we developed a new technology and observed a new pattern of gene expression. We observed that over 72% of RNAs in human genome are expressed in forward-reverse pairs, just like mirror images of each other between forward expression and reverse expression; the overview showed that it cannot be simply described as transcript overlapping, so we designated it as mirror expression. Furthermore, we found that the mirror expression is gene-specific and tissue-specific, and less common in the proximal promoter regions. The size of the shadows varies between different genes, different tissues and different classes. The shadow expression is most significant in the Alu element, it was also observed among L1, Simple Repeats and LTR elements, but rare in other repeats such as low-complexity, LINE/L2, DNA and MIRs. Although there is no evidence yet about the relationship of this mirror pattern and double-strand RNA (dsRNA), this new striking pattern provides a new clue and a new direction to unveil the role of RNAs in the genome functions and diseases.

Temporal small RNA transcriptome profiling unraveled partitioned miRNA expression in developing maize endosperms between reciprocal crosses

In angiosperms, the endosperm nurtures the embryo and provides nutrients for seed germination. To identify the expression pattern of small interfering RNA in the developing maize endosperm, we have performed high-throughput small RNA transcriptome sequencing of kernels at 0, 3, and 5 days after pollination (DAP) and endosperms at 7, 10, and 15 DAP using B73 and Mo17 reciprocal crosses in previous study. Here, we further explored these small RNA-seq data to investigate the potential roles of miRNAs in regulating the gene expression process. In total, 57 conserved miRNAs and 18 novel miRNAs were observed highly expressed in maize endosperm. Temporal expression profiling indicated that these miRNAs exhibited dynamic and partitioned expression patterns at different developmental stages between maize reciprocal crosses, and quantitative RT-PCR results further confirmed our observation. In addition, we found a subset of distinct tandem miRNAs are generated from a single stem-loop structure in maize that might be conserved in monocots. Furthermore, a SNP variation of Zma-miR408-5p at 11th base position was characterized between B73 and Mo17 which might lead to completely different functions in repressing targets. More interestingly, Zma-miR408-5p exhibited B73-biased expression pattern in the B73 and Mo17 reciprocal hybrid endosperms at 7, 10, and 15 DAP according to the reads abundance with SNPs and CAPS experiment. Together, this study suggests that miRNA plays a crucial role in regulating endosperm development, and exhibited distinct expression patterns in developing endosperm between maize reciprocal crosses.

MicroRNA nomenclature and the need for a revised naming prescription

A central environment and interface for microRNA (miRNA) registry and repository and a general standardized framework for their systematic annotation was established over a decade ago. However, the numbers of experimentally and computationally identified miRNAs are swiftly accumulating, and new aspects of miRNA-mediated gene regulation are being revealed. Currently, it is of great significance that the annotation framework should be redefined to include newly discovered miRNA species such as the variants of mature miRNAs (isomiRNAs), and organellar miRNAs: cipomiRNAs and mitomiRNAs. It is also of great importance that key terminology referring to the novelty, evolutionary history or biogenesis of miRNAs, as well as the confidence of their identification are standardized in the literature and disseminated in a central miRNA registry. Here, we review the status of miRNA nomenclature, curation and critical points of need for a revision of miRNA nomenclature and terminology.

Identification and characterization of microRNAs at different flowering developmental stages in moso bamboo (Phyllostachys edulis) by high-throughput sequencing

Researching moso bamboo flowering has been difficult because of its unknown flowering interval and the rarity of florescent samples. To identify microRNAs (miRNAs) and study their expression patterns during the flower developmental process of moso bamboo, small RNAs from non-flowering leaves and four flower developmental periods were sequenced using Illumina technology. In total, 409 known miRNAs and 492 differentially expressed novel miRNAs were identified in moso bamboo. Of the known miRNAs that were differentially expressed between non-flowering and flowering samples, 64 were predicted to have a total of 308 targets. Among the miRNAs, seven known and five novel miRNAs were selected, as were four of their target genes, and their expression profiles were validated using qRT-PCR. The results indicated that the miRNA expression levels were negatively correlated with those of their targets. The research comprehensively revealed that the differentially expressed miRNAs and their targets participated in diverse biological pathways and played significant regulatory roles in moso bamboo flowering. The data provide a significant resource for understanding the molecular mechanisms in moso bamboo flowering and senescence, and serve as the primary foundation for further studies on metabolic regulatory networks that involve miRNAs.

Identification and characterisation of tobacco microRNA transcriptome using high-throughput sequencing

MicroRNAs (miRNAs) are post-transcriptional regulators that are involved in numerous biological processes in plants. In this study, we investigate miRNAs in Honghua Dajinyuan, an agronomically important species of tobacco in China. Here, we report a comprehensive analysis of miRNA expression profiles in the leaf, stem and root using a high-throughput sequencing approach. A total of 165 miRNAs, representing 55 conserved families, and 50 novel miRNAs, representing 19 families, were identified in three libraries. In addition, 12 miRNAs were randomly selected from a differentially expressed conserved miRNA family in three libraries with expression alterations and subjected to qRT-PCR validation. Of these, the expression level of nta-miR167d is highly enriched in the leaf tissue. In addition, the expression level of nta-miR319a is prominently enriched in the stem, while nta-miR160c is highly enriched in the root. Moreover, the target prediction showed that most of the targets coded for transcription factors that are involved in cellular and metabolic processes. GO analysis showed that most of the targets were involved in organelle function, served binding functions, and take part in cellular and metabolic processes. This study helps shed new light on understanding the role of miRNAs in different parts of the tobacco plant and adds a significant number of novel miRNAs to the tobacco miRNA transcriptome.

Lipopolysaccharide perception leads to dynamic alterations in the microtranscriptome of Arabidopsis thaliana cells and leaf tissues

BACKGROUND

MicroRNAs (miRNAs) are non-coding RNA molecules which have recently emerged as important gene regulators in plants and their gene expression analysis is becoming increasingly important. miRNAs regulate gene expression at the post-transcriptional level by translational repression or target degradation of specific mRNAs and gene silencing. In order to profile the microtranscriptome of Arabidopsis thaliana leaf and callus tissues in response to bacterial lipopolysaccharide (LPS), small RNA libraries were constructed at 0 and 3 h post induction with LPS and sequenced by Illumina sequencing technology.

RESULTS

Differential regulation of subset of miRNAs in response to LPS treament was observed. Small RNA reads were mapped to the miRNA database and 358 miRNAs belonging to 49 miRNA families in the callus tissues and 272 miRNAs belonging to 40 miRNA families in the leaf tissues were identified. Moreover, target genes for all the identified miRNAs families in the leaf tissues and 44 of the 49 miRNAs families in the callus tissues were predicted. The sequencing analysis showed that in both callus and leaf tissues, various stress regulated-miRNAs were differentially expressed and real time PCR validated the expression profile of miR156, miR158, miR159, miR169, miR393, miR398, miR399 and miR408 along with their target genes.

CONCLUSION

A. thaliana callus and leaf callus tissues respond to LPS as a microbe-associated molecular pattern molecule through dynamic changes to the microtranscriptome associated with differential transcriptional regulation in support of immunity and basal resistance.

Stress-induced endogenous siRNAs targeting regulatory intron sequences in Brachypodium

Exposure to abiotic stresses triggers global changes in the expression of thousands of eukaryotic genes at the transcriptional and post-transcriptional levels. Small RNA (smRNA) pathways and splicing both function as crucial mechanisms regulating stress-responsive gene expression. However, examples of smRNAs regulating gene expression remain largely limited to effects on mRNA stability, translation, and epigenetic regulation. Also, our understanding of the networks controlling plant gene expression in response to environmental changes, and examples of these regulatory pathways intersecting, remains limited. Here, to investigate the role of smRNAs in stress responses we examined smRNA transcriptomes of Brachypodium distachyon plants subjected to various abiotic stresses. We found that exposure to different abiotic stresses specifically induced a group of novel, endogenous small interfering RNAs (stress-induced, UTR-derived siRNAs, or sutr-siRNAs) that originate from the 3' UTRs of a subset of coding genes. Our bioinformatics analyses predicted that sutr-siRNAs have potential regulatory functions and that over 90% of sutr-siRNAs target intronic regions of many mRNAs in trans. Importantly, a subgroup of these sutr-siRNAs target the important intron regulatory regions, such as branch point sequences, that could affect splicing. Our study indicates that in Brachypodium, sutr-siRNAs may affect splicing by masking or changing accessibility of specific cis-elements through base-pairing interactions to mediate gene expression in response to stresses. We hypothesize that this mode of regulation of gene expression may also serve as a general mechanism for regulation of gene expression in plants and potentially in other eukaryotes.

NPK macronutrients and microRNA homeostasis

Macronutrients are essential elements for plant growth and development. In natural, non-cultivated systems, the availability of macronutrients is not a limiting factor of growth, due to fast recycling mechanisms. However, their availability might be an issue in modern agricultural practices, since soil has been frequently over exploited. From a crop management perspective, the nitrogen (N), phosphorus (P), and potassium (K) are three important limiting factors and therefore frequently added as fertilizers. NPK are among the nutrients that have been reported to alter post-embryonic root developmental processes and consequently, impairs crop yield. To cope with nutrients scarcity, plants have evolved several mechanisms involved in metabolic, physiological, and developmental adaptations. In this scenario, microRNAs (miRNAs) have emerged as additional key regulators of nutrients uptake and assimilation. Some studies have demonstrated the intrinsic relation between miRNAs and their targets, and how they can modulate plants to deal with the NPK availability. In this review, we focus on miRNAs and their regulation of targets involved in NPK metabolism. In general, NPK starvation is related with miRNAs that are involved in root-architectural changes and uptake activity modulation. We further show that several miRNAs were discovered to be involved in plant-microbe symbiosis during N and P uptake, and in this way we present a global view of some studies that were conducted in the last years. The integration of current knowledge about miRNA-NPK signaling may help future studies to focus in good candidates genes for the development of important tools for plant nutritional breeding.