Transcriptome-wide identification of microRNA targets in rice

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

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

The Mh-miR393a-TIR1 module regulates Alternaria alternata resistance of Malus hupehensis mainly by modulating the auxin signaling

miRNAs govern gene expression and regulate plant defense. Alternaria alternata is a destructive fungal pathogen that damages apple. The wild apple germplasm Malus hupehensis is highly resistant to leaf spot disease caused by this fungus. Herein, we elucidated the regulatory and functional role of miR393a in apple resistance against A. alternata by targeting Transport Inhibitor Response 1. Mature miR393 accumulation in infected M. hupehensis increased owing to the transcriptional activation of MIR393a, determined to be a positive regulator of A. alternata resistance to either 'Orin' calli or 'Gala' leaves. 5' RLM-RACE and co-transformation assays showed that the target of miR393a was MhTIR1, a gene encoding a putative F-box auxin receptor that compromised apple immunity. RNA-seq analysis of transgenic calli revealed that MhTIR1 upregulated auxin signaling gene transcript levels and influenced phytohormone pathways and plant-pathogen interactions. miR393a compromised the sensitivity of several auxin-signaling genes to A. alternata infection, whereas MhTIR1 had the opposite effect. Using exogenous indole-3-acetic acid or the auxin synthesis inhibitor L-AOPP, we clarified that auxin enhances apple susceptibility to this pathogen. miR393a promotes SA biosynthesis and impedes pathogen-triggered ROS bursts by repressing TIR1-mediated auxin signaling. We uncovered the mechanism underlying the miR393a-TIR1 module, which interferes with apple defense against A. alternata by modulating the auxin signaling pathway.

Switching action modes of miR408-5p mediates auxin signaling in rice

MicroRNAs (miRNAs) play fundamental roles in many developmental and physiological processes in eukaryotes. MiRNAs in plants generally regulate their targets via either mRNA cleavage or translation repression; however, which approach plays a major role and whether these two function modes can shift remains elusive. Here, we identify a miRNA, miR408-5p that regulates AUXIN/INDOLE ACETIC ACID 30 (IAA30), a critical repressor in the auxin pathway via switching action modes in rice. We find that miR408-5p usually inhibits IAA30 protein translation, but in a high auxin environment, it promotes the decay of IAA30 mRNA when it is overproduced. We further demonstrate that IDEAL PLANT ARCHITECTURE1 (IPA1), an SPL transcription factor regulated by miR156, mediates leaf inclination through association with miR408-5p precursor promoter. We finally show that the miR156-IPA1-miR408-5p-IAA30 module could be controlled by miR393, which silences auxin receptors. Together, our results define an alternative auxin transduction signaling pathway in rice that involves the switching of function modes by miR408-5p, which contributes to a better understanding of the action machinery as well as the cooperative network of miRNAs in plants.

Control of Plant Height and Lateral Root Development via Stu-miR156 Regulation of SPL9 Transcription Factor in Potato

MicroRNAs (miRNAs) are a class of endogenous, non-coding small-molecule RNAs that usually regulate the expression of target genes at the post-transcriptional level. miR156 is one of a class of evolutionarily highly conserved miRNA families. SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factor is one of the target genes that is regulated by miR156. SPL transcription factors are involved in regulating plant growth and development, hormone response, stress response, and photosynthesis. In the present study, transgenic potato plants with overexpressed miR156 were obtained via the Agrobacterium-mediated transformation method. The results showed that the expression levels of the target gene, StSPL9, were all downregulated in the transgenic plants with overexpressed Stu-miR156. Compared with those of the control plants, the plant height and root length of the transgenic plants were significantly decreased, while the number of lateral roots was significantly increased. These results revealed that the miR156/SPLs module was involved in regulating potato plant height and root growth.

Function identification of miR159a, a positive regulator during poplar resistance to drought stress

Drought seriously affects the growth and development of plants. MiR159 is a highly conserved and abundant microRNA family that plays a crucial role in plant growth and stress responses. However, studies of its function in woody plants are still lacking. Here, the expression of miR159a was significantly upregulated after drought treatment in poplar, and the overexpression of miR159a (OX159a) significantly reduced the open area of the stomata and improved water-use efficiency in poplar. After drought treatment, OX159a lines had better scavenging ability of reactive oxygen species and damage of the membrane system was less than that in wild-type lines. MYB was the target gene of miR159a, as verified by psRNATarget prediction, RT-qPCR, degradome sequencing, and 5' rapid amplification of cDNA ends (5' RACE). Additionally, miR159a-short tandem target mimic suppression (STTM) poplar lines showed increased sensitivity to drought stress. Transcriptomic analysis comparing OX159a lines with wild-type lines revealed upregulation of a series of genes related to response to water deprivation and metabolite synthesis. Moreover, drought-responsive miR172d and miR398 were significantly upregulated and downregulated respectively in OX159a lines. This investigation demonstrated that miR159a played a key role in the tolerance of poplar to drought by reducing stomata open area, increasing the number and total area of xylem vessels, and enhancing water-use efficiency, and provided new insights into the role of plant miR159a and crucial candidate genes for the molecular breeding of trees with tolerance to drought stress.

Phased secondary small interfering RNAs in Camellia sinensis var. assamica

Phased secondary small interfering RNAs (phasiRNAs) in plants play important roles in regulating genome stability, plant development and stress adaption. Camellia sinensis var. assamica has immense economic, medicinal and cultural significance. However, there are still no studies of phasiRNAs and their putative functions in this valuable plant. We identified 476 and 43 PHAS loci which generated 4290 twenty one nucleotide (nt) and 264 twenty four nt phasiRNAs, respectively. Moreover, the analysis of degradome revealed more than 35000 potential targets for these phasiRNAs. We identified several conserved 21 nt phasiRNA generation pathways in tea plant, including miR390 → TAS3, miR482/miR2118 → NB-LRR, miR393 → F-box, miR828 → MYB/TAS4, and miR7122 → PPR in this study. Furthermore, we found that some transposase and plant mobile domain genes could generate phasiRNAs. Our results show that phasiRNAs target genes in the same family in cis- or trans-manners, and different members of the same gene family may generate the same phasiRNAs. The phasiRNAs, generated by transposase and plant mobile domain genes, and their targets, suggest that phasiRNAs may be involved in the inhibition of transposable elements in tea plant. To summarize, these results provide a comprehensive view of phasiRNAs in Camellia sinensis var. assamica.

zma-miR159 targets ZmMYB74 and ZmMYB138 transcription factors to regulate grain size and weight in maize

Endosperm cell number is critical in determining grain size in maize (Zea mays). Here, zma-miR159 overexpression led to enlarged grains in independent transgenic lines, suggesting that zma-miR159 contributes positively to maize grain size. Targeting of ZmMYB74 and ZmMYB138 transcription factor genes by zma-miR159 was validated using 5' RACE and dual-luciferase assay. Lines in which ZmMYB74 was knocked out using clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) presented a similar enlarged grain phenotype as those with zma-miR159 overexpression. Downstream genes regulating cell division were identified through DNA affinity purification sequencing using ZmMYB74 and ZmMYB138. Our results suggest that zma-miR159-ZmMYB modules act as an endosperm development hub, participating in the division and proliferation of endosperm cells.

Important Factors Controlling Gibberellin Homeostasis in Plant Height Regulation

Plant height is an important agronomic trait that is closely associated with crop yield and quality. Gibberellins (GAs), a class of highly efficient plant growth regulators, play key roles in regulating plant height. Increasing reports indicate that transcriptional regulation is a major point of regulation of the GA pathways. Although substantial knowledge has been gained regarding GA biosynthetic and signaling pathways, important factors contributing to the regulatory mechanisms homeostatically controlling GA levels remain to be elucidated. Here, we provide an overview of current knowledge regarding the regulatory network involving transcription factors, noncoding RNAs, and histone modifications involved in GA pathways. We also discuss the mechanisms of interaction between GAs and other hormones in plant height development. Finally, future directions for applying knowledge of the GA hormone in crop breeding are described.

Detoxifying the heavy metals: a multipronged study of tolerance strategies against heavy metals toxicity in plants

Heavy metal concentrations exceeding permissible limits threaten human life, plant life, and all other life forms. Different natural and anthropogenic activities emit toxic heavy metals in the soil, air, and water. Plants consume toxic heavy metals from their roots and foliar part inside the plant. Heavy metals may interfere with various aspects of the plants, such as biochemistry, bio-molecules, and physiological processes, which usually translate into morphological and anatomical changes. They use various strategies to deal with the toxic effects of heavy metal contamination. Some of these strategies include restricting heavy metals to the cell wall, vascular sequestration, and synthesis of various biochemical compounds, such as phyto-chelators and organic acids, to bind the free moving heavy metal ions so that the toxic effects are minimized. This review focuses on several aspects of genetics, molecular, and cell signaling levels, which integrate to produce a coordinated response to heavy metal toxicity and interpret the exact strategies behind the tolerance of heavy metals stress. It is suggested that various aspects of some model plant species must be thoroughly studied to comprehend the approaches of heavy metal tolerance to put that knowledge into practical use.

MicroRNAs from Holarrhena pubescens stems: Identification by small RNA Sequencing and their Potential Contribution to Human Gene Targets

Holarrhena pubescens is an effective medicinal plant from the Apocynaceae family, widely distributed over the Indian subcontinent and extensively used by Ayurveda and ethno-medicine systems without apparent side effects. We postulated that miRNAs, endogenous non-coding small RNAs that regulate gene expression at the post-transcriptional level, may, after ingestion into the human body, contribute to the medicinal properties of plants of this species by inducing regulated human gene expression to modulate. However, knowledge is scarce about miRNA in Holarrhena. In addition, to test the hypothesis on the potential pharmacological properties of miRNA, we performed a high-throughput sequencing analysis using the Next Generation Sequencing Illumina platform; 42,755,236 raw reads have been generated from H. pubescens stems from a library of small RNA isolated, identifying 687 known and 50 new miRNAs led. The novel H. pubescens miRNAs were predicted to regulate specific human genes, and subsequent annotations of gene functions suggested a possible role in various biological processes and signaling pathways, such as Wnt, MAPK, PI3K-Akt, and AMPK signaling pathways and endocytosis. The association of these putative targets with many diseases, including cancer, congenital malformations, nervous system disorders, and cystic fibrosis, has been demonstrated. The top hub proteins STAT3, MDM2, GSK3B, NANOG, IGF1, PRKCA, SNAP25, SRSF1, HTT, and SNCA show their interaction with human diseases, including cancer and cystic fibrosis. To our knowledge, this is the first report of uncovering H. pubescens miRNAs based on high-throughput sequencing and bioinformatics analysis. This study has provided new insight into a potential cross-species control of human gene expression. The potential for miRNA transfer should be evaluated as one possible mechanism of action to account for the beneficial properties of this valuable species.

Integration of transcriptome, miRNA and degradome sequencing reveals the early browning response in fresh-cut apple

Surface browning negatively impacts the shelf-life of fresh-cut apple. Herein, we found that the browning of fresh-cut apple aggravated rapidly after 24 h post-cutting, then the transcriptomic and miRNA expression profiles of fresh-cut apple immediately after cutting (T0) and 24 h post-cutting (T24) were analyzed to explore the molecular mechanism of early browning response. A total of 3156 differentially expressed mRNAs (DEGs) and 23 differentially expressed miRNAs (DEmiRNAs) were identified in T24 vs T0. Most DEGs related to respiratory, energy, antioxidant, lipid and secondary metabolism were activated in the early stage of browning. There were 63 target genes of 10 DEmiRNAs validated by degradome sequencing and among them, mdm-miR156aa_L + 1_1 targets 12-oxophytodienoate reductase, ptc-miR6478_R-1 targets patatin-like protein, mdm-miR156aa_L + 1_1 and mdm-miR156aa_L + 1_2 co-target SPLs might participate in the early browning response through regulating antioxidant, lipid and secondary metabolism. Our results will be beneficial for the technological innovation of browning amelioration for fresh-cut apple.

miR319 and its target TCP4 involved in plant architecture regulation in Brassica napus

Plant architecture is a collection of genetically controlled crop productivity and adaptation. MicroRNAs (miRNAs) have been proved to function in various biological processes, but little is known about how miRNA regulates plant architecture in rapeseed (Brassica napus L.). In this study, four small RNA libraries and two degradome libraries from shoot apex of normal and rod-like plants were sequenced. A total of 639 miRNA precursors and 16 differentially expressed miRNAs were identified in this study. In addition, 322 targets were identified through degradome sequencing. Among them, 14 targets were further validated via RNA ligase-mediated 5' rapid amplification of cDNA ends. Transgenic approach showed that increased TCP4 activity in Arabidopsis resulted in premature onset of maturation and reduced plant size along with early flowering and shortened flowering time. miR319-OE lines in Brassica napus exhibited serrated leaves and abnormal development of shoot apical meristem (SAM), which led to the deformed growth of stem and reduced plant height. In conclusion, our study lays the foundation for elucidating miRNA regulate plant architecture and provides new insight into the miR319/TCP4 module regulates plant architecture in rapeseed.

The identification of conserved sequence features of co-translationally decayed mRNAs and upstream open reading frames in angiosperm transcriptomes

RNA turnover is essential in maintaining messenger RNA (mRNA) homeostasis during various developmental stages and stress responses. Co-translational mRNA decay (CTRD), a process in which mRNAs are degraded while still associated with translating ribosomes, has recently been discovered to function in yeast and three angiosperm transcriptomes. However, it is still unclear how prevalent CTRD across the plant lineage. Moreover, the sequence features of co-translationally decayed mRNAs have not been well-studied. Here, utilizing a collection of publicly available degradome sequencing datasets for another seven angiosperm transcriptomes, we have confirmed that CTRD is functioning in at least 10 angiosperms and likely throughout the plant lineage. Additionally, we have identified sequence features shared by the co-translationally decayed mRNAs in these species, implying a possible conserved triggering mechanism for this pathway. Given that degradome sequencing datasets can also be used to identify actively translating upstream open reading frames (uORFs), which are quite understudied in plants, we have identified numerous actively translating uORFs in the same 10 angiosperms. These findings reveal that actively translating uORFs are prevalent in plant transcriptomes, some of which are conserved across this lineage. We have also observed conserved sequence features in the regions flanking these uORFs' stop codons that might contribute to ribosome stalling at these sequences. Finally, we discovered that there were very few overlaps between the mRNAs harboring actively translating uORFs and those sorted into the co-translational decay pathway in the majority of the studied angiosperms, suggesting that these two processes might be nearly mutually exclusive in those species. In total, our findings provide the identification of CTRD and actively translating uORFs across a broad collection of plants and provide novel insights into the important sequence features associated with these collections of mRNAs and regulatory elements, respectively.

Research progress on plant noncoding RNAs in response to low-temperature stress

Low temperature (LT) is an important factor limiting plant growth and distribution. Plants have evolved sophisticated adaptive mechanisms to cope with hypothermia. RNA silencing is the orchestrator of these cellular responses. RNA silencing, which modifies gene expression through noncoding RNAs (ncRNAs), is a strategy used by plants to combat environmental stress. ncRNAs, which have very little protein-coding capacity, work by binding reverse complementary endogenous transcripts. In plants, ncRNAs include small non-coding RNAs (sncRNAs), medium-sized non-coding RNAs (mncRNAs), and long non-coding RNAs (lncRNAs). Apart from describing the biogenesis of different ncRNAs (miRNAs, siRNAs, and lncRNAs), we thoroughly discuss the functions of these ncRNAs during cold acclimation. Two major classes of sncRNAs, microRNAs and siRNAs, play essential regulatory roles in cold response processes through the posttranscriptional gene silencing (PTGS) pathway or transcriptional gene silencing (TGS) pathway. Microarray or transcriptome sequencing analysis can reveal a large number of cold-responsive miRNAs in plants. In this review, the cold-response patterns of miRNAs verified by Northern blotting or quantitative PCR in Arabidopsis thaliana, rice, and many other important crops are discussed. The detailed molecular mechanisms of several miRNAs in Arabidopsis (miR397, miR408, miR402, and miR394) and rice (Osa-miR156, Osa-miR319, and Osa-miR528) that regulate plant cold resistance are elucidated. In addition, the regulatory mechanism of the lncRNA SVALKA in the cold signaling pathway is explained in detail. Finally, we present the challenges for understanding the roles of small ncRNAs in cold signal transduction.

Comprehensive, integrative genomic analysis of microRNA expression profiles in different tissues of two wheat cultivars with different traits

Wheat is one of the most important food sources on Earth. MicroRNAs (miRNA) play important roles in wheat productivity. To identify wheat miRNAs, we constructed and sequenced sRNA libraries from leaves and roots of two wheat cultivars (RAC875 and Kukri) with many different traits. Given that available miRNA wheat complement in the plant-specific database PmiREN ( https://pmiren.com ) does not include root tissues and root-associated miRNAs might thus be missing, we performed first the prediction of novel miRNAs using the sRNAbench tool. We found a total of 150 putatively novel miRNA genes with expression of both arms from 289 unique mature sequences and nearly 30% of all miRNA reads in roots corresponded to novel miRNAs. In contrast, this figure in leaves dropped to under 3%, confirming the undersampling of roots in the complement of known miRNAs. By using 120 publicly available wheat datasets, 598 Zea mays small RNA libraries, 64 plant species genomes, wheat degradome library, and functional enrichment analysis, a subset of novel miRNAs were confirmed as bona-fide miRNAs. Of the total 605 miRNAs identified in this study inclusive of 316 known miRNAs, 528 miRNAs were shared by both cultivars, 429 miRNAs were shared by both root tissues and 329 miRNAs were shared by both leaf tissues. In addition, 32 miRNAs were specific to Kukri while 45 miRNAs were specific to RAC875. These miRNAs had diverse functions, such as regulation of gene transcription, protein translation, energy metabolism, and cell cycle progression. Our data provide a genome-wide miRNA expression profile in these two wheat cultivars and help functional studies of wheat genomics.

Integrated High-Throughput Sequencing, Microarray Hybridization and Degradome Analysis Uncovers MicroRNA-Mediated Resistance Responses of Maize to Pathogen Curvularia lunata

Curvularia lunata (Wakker) Boed, the causal agent of leaf spot in maize, is prone to mutation, making it difficult to control. RNAi technology has proven to be an important tool of genetic engineering and functional genomics aimed for crop improvement. MicroRNAs (miRNAs), which act as post-transcriptional regulators, often cause translational repression and gene silencing. In this article, four small RNA (sRNA) libraries were generated from two maize genotypes inoculated by C. lunata; among these, ltR1 and ltR2 were from the susceptible variety Huangzao 4 (HZ), ltR3 and ltR4, from the resistant variety Luyuan (LY), and 2286, 2145, 1556 and 2504 reads were annotated as miRNA in these four sRNA libraries, respectively. Through the combined analysis of high-throughput sequencing, microarray hybridization and degradome, 48 miRNAs were identified as being related to maize resistance to C. lunata. Among these, PC-732 and PC-169, two new maize miRNAs discovered, were predicted to cleave mRNAs of metacaspase 1 (AMC1) and thioredoxin family protein (Trx), respectively, possibly playing crucial roles in the resistance of maize to C. lunata. To further confirm the role of PC-732 in the interaction of maize and C. lunata, the miRNA was silenced through STTM (short tandem target mimic) technology, and we found that knocking down PC-732 decreased the susceptibility of maize to C. lunata. Precisely speaking, the target gene of PC-732 might inhibit the expression of disease resistance-related genes during the interaction between maize and C. lunata. Overall, the findings of this study indicated the existence of miRNAs involved in the resistance of maize to C. lunata and will contribute to rapidly clarify the resistant mechanism of maize to C. lunata.

Molecular manipulations of miR398 increase rice grain yield under different conditions

Rice miR398 targets two stress-tolerant genes, CSD1-2 (Cu/Zn Superoxide Dismutases1-2) and CCS (copper chaperone of CSD), which usually boost plants' tolerance by inhibiting growth. So, how to accurately regulate the activities of miR398 targets and thus make rice better able to adapt to different conditions has great significances in producing rice yields under the current circumstances of shrinking arable lands resulting from global urbanization and increasing salty soil caused by irrigation. Through controlling the expressions of miR398 in different levels, we found down-regulated expression of miR398 targets can promote growth under good growth conditions while up-regulated expressions of the targets can help rice tolerate salt. In this study, we over-expressed miR398 highly, moderately, and lowly, then three concomitantly inverse levels of its targets' expression were obtained. Under normal growth conditions, the transgenic lines with low and moderate levels of over-expressions of miR398 could increase grain yields 14.5% and 7.3%, respectively, although no transgenic lines could survive well under salty conditions simulating real saline-alkali soil. Using short tandem target mimic (STTM) technology to silence miR398 highly, moderately, and lowly respectively, also three inverse levels of its targets' expression were obtained. All three transgenic lines exhibited good agronomic performances under salt stress in inverse to their degrees of STTM, but their growth was inhibited differently under normal conditions. Altogether, we suggest that flexibly manipulating the expression of miR398 is an ideal strategy to help rice survive better and achieve optimized yields under specific conditions.

Identification and characterization of miRNAome and target genes in Pseudostellaria heterophylla

miRNAs play a crucial role in the development and growth of plants by inhibiting the function of targeted genes at the post-transcription level. However, no miRNAs in Pseudostellaria heterophylla have been reported and their function in the morphogenesis of organs is still unclear. In this study, a total of 159 conserved miRNAs (belonging to 64 families) and 303 level miRNAs were identified from P. heterophylla. Some of them showed specifically up or down-regulated expression in different tissues and numbers of unigenes involved in Plant-pathogen interaction and MAPK signaling pathway-plant were targeted. The significant negative correlation of expression profiles between 30 miRNAs and their target genes (37 unigenes) was observed, respectively. Further, a large number of genes involved with signal transduction of auxin, zeatin, abscisic acid and, jasmonic acid were targeted. Predicated targets of two miRNAs were validated by 5'RLM-RACE, respectively. A large number of mRNAs from four pathogens were targeted by miRNAs from P. heterophylla and some of them were targeted by miR414. In summary, we reported a population of miRNAs from four different vegetative tissues of P. heterophylla by high throughput sequencing, which was analyzed by combining with the constructed transcriptome. These results may help to explain the function of miRNAs in the morphogenesis of organs and defense of pathogens, and may provide theoretical basis for breeding and genetic improvement of P. heterophylla.

Parallel analysis of RNA ends reveals global microRNA-mediated target RNA cleavage in maize

MicroRNAs (miRNAs) are endogenous 20-24-nucleotide non-coding RNAs that play important regulatory roles in many biological processes in eukaryotes. miRNAs modulate the expression of target genes at the post-transcriptional level by transcript cleavage or translational inhibition. The identification of miRNA target genes has been extensively investigated in Arabidopsis and rice, but an in-depth global analysis of miRNA-mediated target regulation is still lacking in maize. Here, we report a transcriptome-wide identification of miRNA targets by analyzing parallel analysis of RNA ends (PARE) datasets derived from nine different tissues at five developmental stages of the maize (Zea mays L.) B73 cultivar. In total, 246 targets corresponding to 60 miRNAs from 25 families were identified, including transcription factors and other genes. In addition, PARE analysis revealed that miRNAs guide specific target transcript cleavage in a tissue-preferential manner. Primary transcripts of MIR159c and MIR169e were found to be cleaved by mature miR159 and miR169, respectively, indicating a negative-feedback regulatory mechanism in miRNA biogenesis. Moreover, several miRNA-target gene pairs involved in seed germination were identified and experimentally validated. Our PARE analyses generated a wide and detailed miRNA-target interaction atlas, which provides a valuable resource for investigating the roles of miRNAs and their targets in maize.

MicroRNA398: A Master Regulator of Plant Development and Stress Responses

MicroRNAs (miRNAs) play crucial roles in plant development and stress responses, and a growing number of studies suggest that miRNAs are promising targets for crop improvement because they participate in the regulation of diverse, important agronomic traits. MicroRNA398 (miR398) is a conserved miRNA in plants and has been shown to control multiple stress responses and plant growth in a variety of species. There are many studies on the stress response and developmental regulation of miR398. To systematically understand its function, it is necessary to summarize the evolution and functional roles of miR398 and its target genes. In this review, we analyze the evolution of miR398 in plants and outline its involvement in abiotic and biotic stress responses, in growth and development and in model and non-model plants. We summarize recent functional analyses, highlighting the role of miR398 as a master regulator that coordinates growth and diverse responses to environmental factors. We also discuss the potential for fine-tuning miR398 to achieve the goal of simultaneously improving plant growth and stress tolerance.

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.

The MicroRNA397a-LACCASE17 module regulates lignin biosynthesis in Medicago ruthenica (L.)

Mechanical strength is essential for the upright growth habit, which is one of the most important characteristics of terrestrial plants. Lignin, a phenylpropanoid-derived polymer mainly present in secondary cell walls plays critical role in providing mechanical support. Here, we report that the prostrate-stem cultivar of the legume forage Medicago ruthenica cultivar 'Mengnong No. 1' shows compromised mechanical strength compared with the erect-stem cultivar 'Zhilixing'. The erect-stem cultivar, 'Zhilixing' has significantly higher lignin content, leading to higher mechanical strength than the prostrate-stem cultivar. The low abundance of miRNA397a in the Zhiixing cultivar causes reduced cleavage of MrLAC17 transcript, which results in enhanced expression level of MrLAC17 compared to that in the prostrate-stem cultivar Mengnong No. 1. Complementation of the Arabidopsis lac4 lac17 double mutants with MrLAC17 restored the lignin content to wild-type levels, confirming that MrLAC17 perform an exchangeable role with Arabidopsis laccases. LAC17-mediated lignin polymerization is therefore increased in the 'Zhilixing', causing the erect stem phenotype. Our data reveal the importance of the miR397a in the lignin biosynthesis and suggest a strategy for molecular breeding targeting plant architecture in legume forage.

microRNAs: Key Players in Plant Response to Metal Toxicity

Environmental metal pollution is a common problem threatening sustainable and safe crop production. Heavy metals (HMs) cause toxicity by targeting key molecules and life processes in plant cells. Plants counteract excess metals in the environment by enhancing defense responses, such as metal chelation, isolation to vacuoles, regulating metal intake through transporters, and strengthening antioxidant mechanisms. In recent years, microRNAs (miRNAs), as a small non-coding RNA, have become the central regulator of a variety of abiotic stresses, including HMs. With the introduction of the latest technologies such as next-generation sequencing (NGS), more and more miRNAs have been widely recognized in several plants due to their diverse roles. Metal-regulated miRNAs and their target genes are part of a complex regulatory network. Known miRNAs coordinate plant responses to metal stress through antioxidant functions, root growth, hormone signals, transcription factors (TF), and metal transporters. This article reviews the research progress of miRNAs in the stress response of plants to the accumulation of HMs, such as Cu, Cd, Hg, Cr, and Al, and the toxicity of heavy metal ions.

MicroRNA miR394 regulates flowering time in Arabidopsis thaliana

miR394 regulates Arabidopsis flowering time in a LCR-independent manner. Arabidopsis plants harboring mutations in theMIR394 genes exhibit early flowering, lower expression of floral repressor FLC and higher expression of floral integrators FT and SOC1. Plant development occurs throughout its entire life cycle and involves a phase transition between vegetative and reproductive phases, leading to the flowering process, fruit formation and ultimately seed production. It has been shown that the microRNA394 (miR394) regulates the accumulation of the transcript coding for LEAF CURLING RESPONSIVENESS, a member of a family of F-Box proteins. The miR394 pathway regulates several processes including leaf morphology and development of the shoot apical meristem during embryogenesis, as well as having been assigned a role in the response to biotic and abiotic stress in Arabidopsis thaliana and other species. Here, we characterized plants harboring mutations in MIR394 precursor genes and demonstrate that mir394a mir394b double mutants display an early flowering phenotype which correlates with a lower expression of FLOWERING LOCUS C earlier in development and higher expression of the floral integrators FLOWERING LOCUS T and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1. Consequently, mutant plants produce fewer branches and exhibit lower seed production. Our work reveals previously unknown developmental aspects regulated by the miR394 pathway, in an LCR-independent manner, contributing to the characterization of the multiple roles of this versatile plant regulatory miRNA.

Molecular and Biochemical Analysis of Duplicated Cytosolic CuZn Superoxide Dismutases of Rice and in silico Analysis in Plants

Superoxide dismutases (SODs, EC 1.15.1.1) are ubiquitous antioxidant metalloenzymes important for oxidative stress tolerance and cellular redox environment. Multiple factors have contributed toward the origin and diversity of SOD isoforms among different organisms. In plants, the genome duplication events, responsible for the generation of multiple gene copies/gene families, have also contributed toward the SOD diversity. However, the importance of such molecular events on the characteristics of SODs has not been studied well. This study investigated the effects of divergence on important characteristics of two block-duplicated rice cytosolic CuZn SODs (OsCSD1, OsCSD4), along with in silico assessment of similar events in other plants. The analysis revealed heterogeneity in gene length, regulatory regions, untranslated regions (UTRs), and coding regions of two OsCSDs. An inconsistency in the database-predicted OsCSD1 gene structure was also identified and validated experimentally. Transcript analysis showed differences in the basal levels and stress responsiveness of OsCSD1 and OsCSD4, and indicated the presence of two transcription start sites in the OsCSD1. At the amino acid level, the two OsCSDs showed differences at 18 sites; however, both exist as a homodimer, displaying typical CuZn SOD characteristics, and enhancing the oxidative stress tolerance of Escherichia coli cells. However, OsCSD4 showed higher specific activity as well as stability. The comparison of the two OsCSDs with reported thermostable CSDs from other plants identified regions likely to be associated with stability, while the homology modeling and superposition highlighted structural differences. The two OsCSDs displayed heteromeric interaction capability and forms an enzymatically active heterodimer (OsCSD1:OsCSD4) on co-expression, which may have significance as both are cytosolic. In silico analysis of 74 plant genomes revealed the prevalence of block duplications for multiple CSD copies (mostly cytosolic). The divergence and clustering analysis of CSDs suggested the possibility of an ancestral duplication event in monocots. Conserved SOD features indicating retention of SOD function among CSD duplicates were evident in few monocots and dicots. In most other species, the CSD copies lacked critical features and may not harbor SOD function; however, other feature-associated functions or novel functions might be present. These aspects of divergent CSD copies encoding co-localized CSDs may have implications in plant SOD functions in the cytosol and other organelles.

The miR528-D3 Module Regulates Plant Height in Rice by Modulating the Gibberellin and Abscisic Acid Metabolisms

Plant height, as one of the important agronomic traits of rice, is closely related to yield. In recent years, plant height-related genes have been characterized and identified, among which the DWARF3 (D3) gene is one of the target genes of miR528, and regulates rice plant height and tillering mainly by affecting strigolactone (SL) signal transduction. However, it remains unknown whether the miR528 and D3 interaction functions in controlling plant height, and the underlying regulatory mechanism in rice. In this study, we found that the plant height, internode length, and cell length of internodes of d3 mutants and miR528-overexpressing (OE-miR528) lines were greatly shorter than WT, D3-overexpressing (OE-D3), and miR528 target mimicry (OE-MIM528) transgenic plants. Knockout of D3 gene (d3 mutants) or miR528-overexpressing (OE-miR528) triggers a substantial reduction of gibberellin (GA) content, but a significant increase of abscisic acid (ABA) accumulation than in WT. The d3 and OE-miR528 transgenic plants were much more sensitive to GA, but less sensitive to ABA than WT. Moreover, the expression level of GA biosynthesis-related key genes, including OsCPS1, OsCPS2, OsKO2 and OsKAO was remarkably higher in OE-D3 plants, while the NECD2 expression, a key gene involved in ABA biosynthesis, was significantly higher in d3 mutants than in WT and OE-D3 plants. The results indicate that the miR528-D3 module negatively regulates plant height in rice by modulating the GA and ABA homeostasis, thereby further affecting the elongation of internodes, and resulting in lower plant height, which adds a new regulatory role to the D3-mediated plant height controlling in rice.

KNOX II transcription factor HOS59 functions in regulating rice grain size

Plant Knotted1-like homeobox (KNOX) genes encode homeodomain-containing transcription factors. In rice (Oryza sativa L.), little is known about the downstream target genes of KNOX Class II subfamily proteins. Here we generated chromatin immunoprecipitation (ChIP)-sequencing datasets for HOS59, a member of the rice KNOX Class II subfamily, and characterized the genome-wide binding sites of HOS59. We conducted trait ontology (TO) analysis of 9705 identified downstream target genes, and found that multiple TO terms are related to plant structure morphology and stress traits. ChIP-quantitative PCR (qPCR) was conducted to validate some key target genes. Meanwhile, our IP-MS datasets showed that HOS59 was closely associated with BELL family proteins, some grain size regulators (OsSPL13, OsSPL16, OsSPL18, SLG, etc.), and some epigenetic modification factors such as OsAGO4α and OsAGO4β, proteins involved in small interfering RNA-mediated gene silencing. Furthermore, we employed CRISPR/Cas9 editing and transgenic approaches to generate hos59 mutants and overexpression lines, respectively. Compared with wild-type plants, the hos59 mutants have longer grains and increased glume cell length, a loose plant architecture, and drooping leaves, while the overexpression lines showed smaller grain size, erect leaves, and lower plant height. The qRT-PCR results showed that mutation of the HOS59 gene led to upregulation of some grain size-related genes such as OsSPL13, OsSPL18, and PGL2. In summary, our results indicate that HOS59 may be a repressor of the downstream target genes, negatively regulating glume cell length, rice grain size, plant architecture, etc. The identified downstream target genes and possible interaction proteins of HOS59 improve our understanding of the KNOX regulatory networks.

AUXIN RESPONSE FACTOR 16 (StARF16) regulates defense gene StNPR1 upon infection with necrotrophic pathogen in potato

We demonstrate a new regulatory mechanism in the jasmonic acid (JA) and salicylic acid (SA) mediated crosstalk in potato defense response, wherein, miR160 target StARF16 (a gene involved in growth and development) binds to the promoter of StNPR1 (a defense gene) and negatively regulates its expression to suppress the SA pathway. Overall, our study establishes the importance of StARF16 in regulation of StNPR1 during JA mediated defense response upon necrotrophic pathogen interaction. Plants employ antagonistic crosstalk between salicylic acid (SA) and jasmonic acid (JA) to effectively defend them from pathogens. During biotrophic pathogen attack, SA pathway activates and suppresses the JA pathway via NONEXPRESSOR OF PATHOGENESIS-RELATED GENES 1 (NPR1). However, upon necrotrophic pathogen attack, how JA-mediated defense response suppresses the SA pathway, is still not well-understood. Recently StARF10 (AUXIN RESPONSE FACTOR), a miR160 target, has been shown to regulate SA and binds to the promoter of StGH3.6 (GRETCHEN HAGEN3), a gene proposed to maintain the balance between the free SA and auxin in plants. In the current study, we investigated the role of StARF16 (a miR160 target) in the regulation of the defense gene StNPR1 in potato upon activation of the JA pathway. We observed that a negative correlation exists between StNPR1 and StARF16 upon infection with the pathogen. The results were further confirmed through the exogenous application of SA and JA. Using yeast one-hybrid assay, we demonstrated that StARF16 binds to the StNPR1 promoter through putative ARF binding sites. Additionally, through protoplast transfection and chromatin immunoprecipitation experiments, we showed that StARF16 could bind to the StNPR1 promoter and regulate its expression. Co-transfection assays using promoter deletion constructs established that ARF binding sites are present in the 2.6 kb sequence upstream to the StNPR1 gene and play a key role in its regulation during infection. In summary, we demonstrate the importance of StARF16 in the regulation of StNPR1, and thus SA pathway, during JA-mediated defense response upon necrotrophic pathogen interaction.

MicroRNAs Are Involved in Regulating Plant Development and Stress Response through Fine-Tuning of TIR1/AFB-Dependent Auxin Signaling

Auxin, primarily indole-3-acetic acid (IAA), is a versatile signal molecule that regulates many aspects of plant growth, development, and stress response. Recently, microRNAs (miRNAs), a type of short non-coding RNA, have emerged as master regulators of the auxin response pathways by affecting auxin homeostasis and perception in plants. The combination of these miRNAs and the autoregulation of the auxin signaling pathways, as well as the interaction with other hormones, creates a regulatory network that controls the level of auxin perception and signal transduction to maintain signaling homeostasis. In this review, we will detail the miRNAs involved in auxin signaling to illustrate its in planta complex regulation.

Soybean miR159-GmMYB33 Regulatory Network Involved in Gibberellin-Modulated Resistance to Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is an obligate sedentary biotroph that poses major threats to soybean production globally. Recently, multiple miRNAome studies revealed that miRNAs participate in complicated soybean-SCN interactions by regulating their target genes. However, the functional roles of miRNA and target genes regulatory network are still poorly understood. In present study, we firstly investigated the expression patterns of miR159 and targeted GmMYB33 genes. The results showed miR159-3p downregulation during SCN infection; conversely, GmMYB33 genes upregulated. Furthermore, miR159 overexpressing and silencing soybean hairy roots exhibited strong resistance and susceptibility to H. glycines, respectively. In particular, miR159-GAMYB genes are reported to be involve in GA signaling and metabolism. Therefore, we then investigated the effects of GA application on the expression of miR159-GAMYB module and the development of H. glycines. We found that GA directly controls the miR159-GAMYB module, and exogenous GA application enhanced endogenous biologically active GA1 and GA3, the abundance of miR159, lowered the expression of GmMYB33 genes and delayed the development of H. glycines. Moreover, SCN infection also results in endogenous GA content decreased in soybean roots. In summary, the soybean miR159-GmMYB33 module was directly involved in the GA-modulated soybean resistance to H. glycines.

miRNAomic Approach to Plant Nitrogen Starvation

Nitrogen (N) is one of the indispensable nutrients required by plants for their growth, development, and survival. Being a limited nutrient, it is mostly supplied exogenously to the plants, to maintain quality and productivity. The increased use of N fertilizers is associated with high-cost inputs and negative environmental consequences, which necessitates the development of nitrogen-use-efficient plants for sustainable agriculture. Understanding the regulatory mechanisms underlying N metabolism in plants under low N is one of the prerequisites for the development of nitrogen-use-efficient plants. One of the important and recently discovered groups of regulatory molecules acting at the posttranscriptional and translational levels are microRNAs (miRNAs). miRNAs are known to play critical roles in the regulation of gene expression in plants under different stress conditions including N stress. Several classes of miRNAs associated with N metabolism have been identified so far. These nitrogen-responsive miRNAs may provide a platform for a better understanding of the regulation of N metabolism and pave a way for the development of genotypes for better N utilization. The current review presents a brief outline of miRNAs and their regulatory role in N metabolism.

Temporal expression study of miRNAs in the crown tissues of winter wheat grown under natural growth conditions

BACKGROUND

Winter wheat requires prolonged exposure to low temperature to initiate flowering (vernalization). Shoot apical meristem of the crown is the site of cold perception, which produces leaf primordia during vegetative growth before developing into floral primordia at the initiation of the reproductive phase. Although many essential genes for winter wheat cold acclimation and floral initiation have been revealed, the importance of microRNA (miRNA) meditated post-transcriptional regulation in crowns is not well understood. To understand the potential roles of miRNAs in crown tissues, we performed a temporal expression study of miRNAs in crown tissues at the three-leaf stage, winter dormancy stage, spring green-up stage, and jointing stage of winter wheat grown under natural growth conditions.

RESULTS

In total, 348 miRNAs belonging to 298 miRNA families, were identified in wheat crown tissues. Among them, 92 differentially expressed miRNAs (DEMs) were found to be significantly regulated from the three-leaf stage to the jointing stage. Most of these DEMs were highly expressed at the three-leaf stage and winter dormancy stage, and then declined in later stages. Six DEMs, including miR156a-5p were markedly induced during the winter dormancy stage. Eleven DEMs, including miR159a.1, miR390a-5p, miR393-5p, miR160a-5p, and miR1436, were highly expressed at the green-up stage. Twelve DEMs, such as miR172a-5p, miR394a, miR319b-3p, and miR9676-5p were highly induced at the jointing stage. Moreover, 14 novel target genes of nine wheat or Pooideae-specific miRNAs were verified using RLM-5' RACE assay. Notably, six mTERFs and two Rf1 genes, which are associated with mitochondrial gene expression, were confirmed as targets of three wheat-specific miRNAs.

CONCLUSIONS

The present study not only confirmed the known miRNAs associated with phase transition and floral development, but also identified a number of wheat or Pooideae-specific miRNAs critical for winter wheat cold acclimation and floral development. Most importantly, this study provided experimental evidence that miRNA could regulate mitochondrial gene expression by targeting mTERF and Rf1 genes. Our study provides valuable information for further exploration of the mechanism of miRNA mediated post-transcriptional regulation during winter wheat vernalization and inflorescent initiation.

Lateral root formation and nutrients: nitrogen in the spotlight

Lateral roots are important to forage for nutrients due to their ability to increase the uptake area of a root system. Hence, it comes as no surprise that lateral root formation is affected by nutrients or nutrient starvation, and as such contributes to the root system plasticity. Understanding the molecular mechanisms regulating root adaptation dynamics toward nutrient availability is useful to optimize plant nutrient use efficiency. There is at present a profound, though still evolving, knowledge on lateral root pathways. Here, we aimed to review the intersection with nutrient signaling pathways to give an update on the regulation of lateral root development by nutrients, with a particular focus on nitrogen. Remarkably, it is for most nutrients not clear how lateral root formation is controlled. Only for nitrogen, one of the most dominant nutrients in the control of lateral root formation, the crosstalk with multiple key signals determining lateral root development is clearly shown. In this update, we first present a general overview of the current knowledge of how nutrients affect lateral root formation, followed by a deeper discussion on how nitrogen signaling pathways act on different lateral root-mediating mechanisms for which multiple recent studies yield insights.

ZmmiR190 and its target regulate plant responses to drought stress through an ABA-dependent pathway

MicroRNAs (miRNAs) are small, non-coding regulatory RNAs that regulate gene expression by facilitating target mRNA cleavage in plants. They are crucial for responses to diverse stresses. The novel drought-responsive miRNA ZmmiR190 was previously identified during an analysis of the maize transcriptome. In this study, we revealed that transgenic Arabidopsis thaliana overexpressing ZmmiR190 is more sensitive to drought than the wild-type control. The transcript of a nuclear-localized gene, ZmCRP04, was identified as a likely target of ZmmiR190. Moreover, ZmmiR190 and ZmCRP04 had the opposite expression profiles following drought and salt treatments. Additionally, 5' RACE and coexpression analyses in A. thaliana provided evidence of the in vivo targeting of the ZmCRP04 transcript by ZmmiR190. Furthermore, the overexpression of ZmCRP04 in A. thaliana and rice significantly enhanced drought tolerance, with lower malonaldehyde contents and relative electrolyte leakage in the transgenic A. thaliana and rice plants than in the wild-type control. Transgenic plants overexpressing ZmmiR190 or ZmCRP04 were hypersensitive to abscisic acid. These results suggest that the ZmCRP04 transcript is targeted by ZmmiR190 and may encode a protein that positively regulates drought stress tolerance via an abscisic acid-dependent pathway. These findings may be relevant for future molecular breeding aimed at improving crop drought tolerance.

Identification and Characterization of Downy Mildew-Responsive microRNAs in Indian Vitis vinifera by High-Throughput Sequencing

Downy mildew (DM) is one of the most devastating diseases disturbing viticulture, mainly during temperate and humid climates. The DM pathogen can attack grapevine leaves and berries differentially, and the disease is managed with recurring applications of fungicides that direct pathogen pressure, develop of resistant strains, and lead to residual soil toxicity and increased pollution effects. Plant microRNAs (miRNAs) are important candidates in physiological regulatory roles in response to biotic stress in plants. In this study, high-throughput sequencing and MiRDeep-P were employed to identify miRNAs in Vitis vinifera. Altogether, 22,492,910, 25,476,471, and 22,448,438 clean reads from the sterile distilled water (SDW)-control, bio-pesticide Trichoderma harzianum (TriH_JSB36)-treated, and downy mildew Plasmopara viticola pathogen libraries, respectively, were obtained. On the basis of the sequencing results and analysis (differential expression analysis), we observed significant differences in 15 miRNAs (5 novel upregulated, and 10 known downregulated) in the pathogen-infected sample (Test) in comparison to the SDW-control sample, with majority of the reads beingin the range of 20-24 bp. This study involves the identification and characterization of vvi-miRNAs that are involved in resistance against downy mildew disease in grapes.

Investigations on Regulation of MicroRNAs in Rice Reveal [Ca2+]cyt Signal Transduction Regulated MicroRNAs

MicroRNAs (miRNAs) are critical components of the multidimensional regulatory networks in eukaryotic systems. Given their diverse spectrum of function, it is apparent that the transcription, processing, and activity of the miRNAs themselves, is very dynamically regulated. One of the most important and universally implicated signaling molecules is [Ca²⁺]_cyt. It is known to regulate a plethora of developmental and metabolic processes in both plants and animals; however, its impact on the regulation of miRNA expression is relatively less explored. The current study employed a combination of internal and external calcium channel inhibitors to establishing that [Ca²⁺]_cyt signatures actively regulate miRNA expression in rice. Involvement of [Ca²⁺]_cyt in the regulation of miRNA expression was further confirmed by treatment with calcimycin, the calcium ionophore. Modulation of the cytosolic calcium levels was also found to regulate the drought-responsive expression as well as ABA-mediated response of miRNA genes in rice seedlings. The study further establishes the role of calmodulins and Calmodulin-binding Transcription Activators (CAMTAs) as important components of the signal transduction schema that regulates miRNA expression. Yeast one-hybrid assay established that OsCAMTA4 & 6 are involved in the transcriptional regulation of miR156a and miR167h. Thus, the study was able to establish that [Ca²⁺]_cyt is actively involved in regulating the expression of miRNA genes both under control and stress conditions.

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.

Mdm-MIR393b-mediated adventitious root formation by targeted regulation of MdTIR1A expression and weakened sensitivity to auxin in apple rootstock

Adventitious root (AR) formation is of great significance for apple rootstock breeding. It is widely accepted that miR393 influences AR formation in many plant species; however, the molecular mechanism by which factors regulate AR formation remains insufficient. In this study, the evolutionary relationship of mdm-miR393 and candidate target genes MdTIR1/AFB was systematically identified, and the expression patterns were analysed. Multisequence alignment analysis of miR393 family members suggests that miR393 conservatively evolved between different species. The evolutionary relationship of the TIR1/AFBs can be divided into G1, G2 and G3 subgroups. During AR formation, the expression level of mdm-miR393a/b/c was significantly upregulated at 1 d and 7 d by exogenous auxin treatment. Furthermore, the expression levels of MdTIR1A, MdTIR1D, MdAFB1, MdAFB2, MdAFB3, MdAFB4 and MdAFB8 also appeared to be significantly changed by exogenous auxin induction. Subsequently, tissue-specific expression analysis showed that the expression levels of mdm-miR393 and MdTIR1/AFBs in different tissues exhibited significant differences. The promoter of mdm-miR393 contains multiple elements that respond to ABA, adversity and light signals; auxin treatment can activate the mdm-MIR393b promoter but is obviously inhibited by NPA treatment. The targeting relationship between mdm-MIR393b and MdTIR1A was verified by expression patterns, degradation group data, transient tobacco conversion results, and genes functions experiments. Heterologous overexpression of mdm-MIR393b (35S::mdm-MIR393b) decreased the number of ARs in the phenotype and reduced the expression level of the target gene NtTIR1 in tobacco. Compared to the wild type, the 35S::mdm-MIR393b transgenic plants demonstrated insensitivity to auxin. Furthermore, tir1 mutant exhibited reduced root system structure relative to the control. The above results illustrated that mdm-MIR393b is involved in mediating AR formation by targeted regulation of MdTIR1A expression in apple rootstock.

Expression Alteration of Candidate Rice MiRNAs in Response to Sheath Blight Disease

Background:

MicroRNAs, as small non-coding RNAs, are recently reported to be involved in plant defense system against pathogens including fungi.

Objective:

In this research, it was intended to investigate candidate susceptible rice (Oryza Sativa) Osa-miRNA expression alteration following the infection by Rhizoctonia solani.

Materials and Methods:

To this aim, literature review suggested eight conserved plant miRNAs that are involved in other plant-pathogen interactions. Then, sixty days old rice plants (Hashemi, susceptible cultivar) were inoculated with R. solani and candidate miRNA expression alterations were investigated 2 hpi (hours post inoculation), 2 dpi (days post inoculation) and 6 dpi.

Results:

RT-qPCR analysis suggested four subgroups of candidate miRNAs based on the time of their responses to the pathogenesis of R. solani. While Osa-miR-156 was early-responsive, Osa-miR159 was the last-responsive and Osa-miR167, Osa-miR171, Osa-miR408, and Osa-miR444 were late responsive to R. solani infection. Osa-miR166 and Osa-miR393 were non-responsive to this infection, compared to the mock-inoculated control group. Consistently, Os-SPL3 and Os-MADS known target genes were expressed in reverse correlation to Osa-miR156 and Osa-miR444, respectively.

Conclusions:

From these data, it is suggested that both early (Osa-miR-156) and late (Osa-miR167, Osa-miR171, Osa- miR408, Osa-miR444) responsive miRNAs might be involved in R. solani infection in rice plants.

MicroRNAs Involved in Regulatory Cytoplasmic Male Sterility by Analysis RNA-seq and Small RNA-seq in Soybean

Cytoplasmic male sterility (CMS) is an important plant characteristic for exploiting heterosis to enhance crop traits during breeding. However, the CMS regulatory network remains unclear in plants, even though researchers have attempted to isolate genes associated with CMS. In this study, we performed high-throughput sequencing and degradome analyses to identify microRNAs (miRNAs) and their targets in a soybean CMS line (JLCMS9A) and its maintainer line (JLCMS9B). Additionally, the differentially expressed genes during reproductive development were identified using RNA-seq data. A total of 280 miRNAs matched soybean miRNA sequences in miRBase, including mature miRNAs and pre-miRNAs. Of the 280 miRNAs, 30, 23, and 21 belonged to the miR166, miR156, and miR171 families, respectively. Moreover, 410 novel low-abundant miRNAs were identified in the JLCMS9A and JLCMS9B flower buds. Furthermore, 303 and 462 target genes unique to JLCMS9A and JLCMS9B, respectively, as well as 782 common targets were predicted based on the degradome analysis. Target genes differentially expressed between the CMS line and the maintainer line were revealed by an RNA-seq analysis. Moreover, all target genes were annotated with diverse functions related to biological processes, cellular components, and molecular functions, including transcriptional regulation, the nucleus, meristem maintenance, meristem initiation, cell differentiation, auxin-activated signaling, plant ovule development, and anther development. Finally, a network was built based on the interactions. Analyses of the miRNA, degradome, and transcriptome datasets generated in this study provided a comprehensive overview of the reproductive development of a CMS soybean line. The data presented herein represent useful information for soybean hybrid breeding. Furthermore, the study results indicate that miRNAs might contribute to the soybean CMS regulatory network by modulating the expression of CMS-related genes. These findings lay the foundation for future studies on the molecular mechanisms underlying soybean CMS.

Splice Variants of Superoxide Dismutases in Rice and Their Expression Profiles under Abiotic Stresses

The superoxide dismutases (SODs) play vital roles in controlling cellular reactive oxygen species (ROS) that are generated both under optimal as well as stress conditions in plants. The rice genome harbors seven SOD genes (CSD1, CSD2, CSD3, CSD4, FSD1, FSD2, and MSD) that encode seven constitutive transcripts. Of these, five (CSD2, CSD3, CSD4, FSD1, and MSD) utilizes an alternative splicing (AS) strategy and generate seven additional splice variants (SVs) or mRNA variants, i.e., three for CSD3, and one each for CSD2, CSD4, FSD1, and MSD. The exon-intron organization of these SVs revealed variations in the number and length of exons and/or untranslated regions (UTRs). We determined the expression patterns of SVs along with their constitutive forms of SODs in rice seedlings exposed to salt, osmotic, cold, heavy metal (Cu⁺²) stresses, as well as copper-deprivation. The results revealed that all seven SVs were transcriptionally active in both roots and shoots. When compared to their corresponding constitutive transcripts, the profiles of five SVs were almost similar, while two specific SVs (CSD3-SV4 and MSD-SV2) differed significantly, and the differences were also apparent between shoots and roots suggesting that the specific SVs are likely to play important roles in a tissue-specific and stress-specific manner. Overall, the present study has provided a comprehensive analysis of the SVs of SODs and their responses to stress conditions in shoots and roots of rice seedlings.

MicroRNAs: Potential Targets for Developing Stress-Tolerant Crops

Crop yield is challenged every year worldwide by changing climatic conditions. The forecasted climatic scenario urgently demands stress-tolerant crop varieties to feed the ever-increasing global population. Molecular breeding and genetic engineering approaches have been frequently exploited for developing crops with desired agronomic traits. Recently, microRNAs (miRNAs) have emerged as powerful molecules, which potentially serve as expression markers during stress conditions. The miRNAs are small non-coding endogenous RNAs, usually 20-24 nucleotides long, which mediate post-transcriptional gene silencing and fine-tune the regulation of many abiotic- and biotic-stress responsive genes in plants. The miRNAs usually function by specifically pairing with the target mRNAs, inducing their cleavage or repressing their translation. This review focuses on the exploration of the functional role of miRNAs in regulating plant responses to abiotic and biotic stresses. Moreover, a methodology is also discussed to mine stress-responsive miRNAs from the enormous amount of transcriptome data available in the public domain generated using next-generation sequencing (NGS). Considering the functional role of miRNAs in mediating stress responses, these molecules may be explored as novel targets for engineering stress-tolerant crop varieties.

MicroRNAs modulate ethylene induced retrograde signal for rice endosperm starch biosynthesis by default expression of transcriptome

Control of stage specific spike in ethylene production at anthesis has been a vauable route to potentially enhance genetic ceiling for grain filling of rice spikelet. A number of genes controlling ethylene homeostasis and starch synthesis have been identified so long, but lack of credible information on master modulation of gene expression by miRNAs and their target genes associated with hormonal dynamics obfuscate mechanisms controlling genotype difference in quantum of grain filling. The confusion accounts for consequent shrinkage of options for yield manipulation. In a two by two factorial design, miRNA regulation of spikelet specific grain development in low against high sterile recombinant inbred lines of rice Oryza sativa L. namely CR 3856-62-11-3-1-1-1-1-1-1 (SR 157) and CR 3856-63-1-1-1-1-1-1-1 (SR 159) respectively, and inferior verses superior spikelets were compared during first 10 days after anthesis. Grain filling was poorer in SR159 than SR157 and inferior spikelets in the former were most vulnerable. Between the cultivars, overall expression of unique miRNAs with targets on ethylene pathway genes was higher in SR159 than SR157 and the situation was opposite for auxin pathway genes. Precision analysis in psTarget server database identified up-regulation of MIR2877 and MIR530-5p having Os11t0141000-02 and Os07t0239400-01 (PP2A regulatory subunit-like protein and ethylene-responsive small GTP-binding proteins) and MIR396h having Os01t0643300-02 (an auxin efflux carrier protein) and Os01t0643300-01 (a PIN1-like auxin transport protein), as targets with highest probability at anthesis and 5 days after anthesis respectively, in the inferior spikelet and the fold change values of DGE matched with pattern of gene expression (relative transcript level) in the qRT-PCR studies conducted for relevant miRNAs and protein factors for ethylene and auxin signalling. In conclusion, epigenetic regulation of both auxin and ethylene homeostasis control grain filling of rice spikelet was established, but evidences were more robust for the latter.

Genome-wide miRNA analysis and integrated network for flavonoid biosynthesis in Osmanthus fragrans

Background

Osmanthus fragrans is an important economical plant containing multiple secondary metabolites including flavonoids and anthocyanins. During the past years, the roles of miRNAs in regulating the biosynthesis of secondary metabolites in plants have been widely investigated. However, few studies on miRNA expression profiles and the potential roles in regulating flavonoid biosynthesis have been reported in O. fragrans.

Results

In this study, we used high-throughput sequencing technology to analyze the expression profiles of miRNAs in leaf and flower tissues of O. fragrans. As a result, 106 conserved miRNAs distributed in 47 families and 88 novel miRNAs were identified. Further analysis showed there were 133 miRNAs differentially expressed in leaves and flowers. Additionally, the potential target genes of miRNAs as well as the related metabolic pathways were predicted. In the end, flavonoid content was measured in flower and leaf tissues and potential role of miR858 in regulating flavonoid synthesis was illustrated in O. fragrans.

Conclusions

This study not only provided the genome-wide miRNA profiles in the flower and leaf tissue of O. fragrans, but also investigated the potential regulatory role of miR858a in flavonoid synthesis in O. fragrans. The results specifically indicated the connection of miRNAs to the regulation of secondary metabolite biosynthesis in non-model economical plant.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07439-y.

Suppression of rice miR168 improves yield, flowering time and immunity

MicroRNA168 (miR168) is a key miRNA that targets Argonaute1 (AGO1), a major component of the RNA-induced silencing complex^1,2. Previously, we reported that miR168 expression was responsive to infection by Magnaporthe oryzae, the causal agent of rice blast disease³. However, how miR168 regulates immunity to rice blast and whether it affects rice development remains unclear. Here, we report our discovery that the suppression of miR168 by a target mimic (MIM168) not only improves grain yield and shortens flowering time in rice but also enhances immunity to M. oryzae. These results were validated through repeated tests in rice fields in the absence and presence of rice blast pressure. We found that the miR168-AGO1 module regulates miR535 to improve yield by increasing panicle number, miR164 to reduce flowering time, and miR1320 and miR164 to enhance immunity. Our discovery demonstrates that changes in a single miRNA enhance the expression of multiple agronomically important traits.

Degradome sequencing-based identification of phasiRNAs biogenesis pathways in Oryza sativa

BACKGROUND

The microRNAs(miRNA)-derived secondary phased small interfering RNAs (phasiRNAs) participate in post-transcriptional gene silencing and play important roles in various bio-processes in plants. In rice, two miRNAs, miR2118 and miR2275, were mainly responsible for triggering of 21-nt and 24-nt phasiRNAs biogenesis, respectively. However, relative fewer phasiRNA biogenesis pathways have been discovered in rice compared to other plant species, which limits the comprehensive understanding of phasiRNA biogenesis and the miRNA-derived regulatory network.

RESULTS

In this study, we performed a systematical searching for phasiRNA biogenesis pathways in rice. As a result, five novel 21-nt phasiRNA biogenesis pathways and five novel 24-nt phasiRNA biogenesis pathways were identified. Further investigation of their regulatory function revealed that eleven novel phasiRNAs in 21-nt length recognized forty-one target genes. Most of these genes were involved in the growth and development of rice. In addition, five novel 24-nt phasiRNAs targeted to the promoter of an OsCKI1 gene and thereafter resulted in higher level of methylation in panicle, which implied their regulatory function in transcription of OsCKI1,which acted as a regulator of rice development.

CONCLUSIONS

These results substantially extended the information of phasiRNA biogenesis pathways and their regulatory function in rice.

Profiling of MicroRNAs Involved in Mepiquat Chloride-Mediated Inhibition of Internode Elongation in Cotton (Gossypium hirsutum L.) Seedlings

Mepiquat chloride (MC) is the most important plant growth retardant that is widely used in cotton (Gossypium hirsutum L.) production to suppress excessive vegetative growth and improve plant architecture. MicroRNAs (miRNAs) are important gene expression regulators that control plant growth and development. However, miRNA-mediated post-transcriptional regulation in MC-induced growth inhibition remains unclear. In this study, the dynamic expression profiles of miRNAs responsive to MC in cotton internodes were investigated. A total of 508 known miRNAs belonging to 197 families and five novel miRNAs were identified. Among them, 104 miRNAs were differentially expressed at 48, 72, or 96 h post MC treatment compared with the control (0 h); majority of them were highly conserved miRNAs. The number of differentially expressed miRNAs increased with time after treatment. The expression of 14 known miRNAs was continuously suppressed, whereas 12 known miRNAs and one novel miRNA were continuously induced by MC. The expression patterns of the nine differentially expressed miRNAs were verified using qRT-PCR. The targets of the known and novel miRNAs were predicted. Four conserved and six novel targets were validated using the RLM-5' RACE assay. This study revealed that miRNAs play crucial regulatory roles in the MC-induced inhibition of internode elongation. It can improve our understanding of post-transcriptional gene regulation in MC-mediated growth inhibition and could potentially facilitate the breeding of dwarf cotton.

A Cotton Lignin Biosynthesis Gene, GhLAC4, Fine-Tuned by ghr-miR397 Modulates Plant Resistance Against Verticillium dahliae

Plant lignin is a component of the cell wall, and plays important roles in the transport potential of water and mineral nutrition and plant defence against biotic stresses. Therefore, it is necessary to identify lignin biosynthesis-related genes and dissect their functions and underlying mechanisms. Here, we characterised a cotton LAC, GhLAC4, which participates in lignin biosynthesis and plant resistance against Verticillium dahliae. According to degradome sequencing and GUS reporter analysis, ghr-miR397 was identified to directedly cleave the GhLAC4 transcript through base complementary. GhLAC4 knockdown and ghr-miR397 overexpression significantly reduced basal lignin content compared to the control, whereas ghr-miR397 silencing significantly increased basal lignin levels. Based on staining patterns and GC/MS analysis, GhLAC4 acted in G-lignin biosynthesis. Under V. dahliae infection, we found that G-lignin content in ghr-miR397-knockdowned plants significantly increased, compared to these plants under the mock treatment, while G-lignin contents in GhLAC4-silenced plants and ghr-miR397-overexpressed plants treated with pathogen were comparable with these plants treated with mock, indicating that GhLAC4 participates in defence-induced G-lignin biosynthesis in the cell wall. Knockdown of ghr-miR397 in plants inoculated with V. dahliae promoted lignin accumulation and increased plant resistance. The overexpression of ghr-miR397 and knockdown of GhLAC4 reduced lignin content and showed higher susceptibility of plants to the fungal infection compared to the control. The extract-free stems of ghr-miR397-knockdowned plants lost significantly less weight when treated with commercial cellulase and V. dahliae secretion compared to the control, while the stems of ghr-miR397-overexpressed and GhLAC4-silenced plants showed significantly higher loss of weight. These results suggest that lignin protects plant cell walls from degradation mediated by cellulase or fungal secretions. In summary, the ghr-miR397-GhLAC4 module regulates both basal lignin and defence-induced lignin biosynthesis and increases plant resistance against infection by V. dahliae.

In planta transcriptome analysis reveals tissue-specific expression of pathogenicity genes and microRNAs during rice-Magnaporthe interactions

Transcriptional re-programming in host and pathogen upon leaf and neck infection is an evolving area of research for the rice blast community. Analysis of in planta rice transcriptome in leaf and neck tissues revealed tissue-specific and infection-specific expression of rice and Magnaporthe oryzae genes in host and pathogen. The glycosyl hydrolase, isocitrate lyase, cupin domain containing protein, TF2, CMPG1, CHIT17 and OsCML14 genes were uniquely expressed in leaf infection. Genes like cytochrome P450, inhibitor I family protein, GSTU6, abscisic stress ripening, and cupin domain containing protein were up-regulated during neck infection. In our microRNA sequencing study, Osa-miR166n-3p was highly expressed in upon Magnaporthe leaf infection, whereas osa-miR1661-3p, osa-miR166n-3p and osa-miR159b were overexpressed in neck infection. Here we report several transcripts being targeted by up and down regulated microRNAs during infection. The putative genes expressed upon infection in leaf and neck could be used in understanding the dual-epidemics of blast disease.

Integration of small RNA, degradome and proteome sequencing in Oryza sativa reveals a delayed senescence network in tetraploid rice seed

Seed of rice is an important strategic resource for ensuring the security of China's staple food. Seed deterioration as a result of senescence is a major problem during seed storage, which can cause major economic losses. Screening among accessions in rice germplasm resources for traits such as slow senescence and increased seed longevity during storage is, therefore, of great significance. However, studies on delayed senescence in rice have been based mostly on diploid rice seed to date. Despite better tolerance have been verified by the artificial aging treatment for polyploid rice seed, the delayed senescence properties and delayed senescence related regulatory mechanisms of polyploid rice seed are rarely reported, due to the lack of polyploid rice materials with high seed set. High-throughput sequencing was applied to systematically investigate variations in small RNAs, the degradome, and the proteome between tetraploid and diploid rice seeds. Degradome sequencing analysis of microRNAs showed that expression of miR-164d, which regulates genes encoding antioxidant enzymes, was changed significantly, resulting in decreased miRNA-mediated cleavage of target genes in tetraploid rice. Comparisons of the expression levels of small RNAs (sRNAs) in the tetraploid and diploid libraries revealed that 12 sRNAs changed significantly, consistent with the findings from degradome sequencing. Furthermore, proteomics also showed that antioxidant enzymes were up-regulated in tetraploid rice seeds, relative to diploids.