Evolution of plant microRNAs and their targets

The role of microRNAs in NBS-LRR gene expression and its implications for plant immunity and crop development

Plants evolved, over millions of years, complex defense systems against pathogens. Once infected, the interaction between pathogen effector molecules and host receptors triggers plant immune responses, which include apoptosis, systemic immune response, among others. An important protein family responsible for pathogen effector recognition is the nucleotide binding site-leucine repeat rich (NBS-LRR) proteins. The NBS-LRR gene family is the largest disease resistance gene class in plants. These proteins are widely distributed in vascular plants and have a complex multigenic cluster distribution in plant genomes. To counteract the genetic load of such a large gene family on fitness cost, plants evolved a mechanism using post transcriptional gene silencing induced by small RNAs, particularly microRNAs. For the NBS-LRR gene family, the small RNAs involved in this silencing mechanism are mainly the microRNA482/2118 superfamily. This suppression mechanism is relieved upon pathogen infection, thus allowing increased NBS-LRR expression and triggering plant immunity. In this review, we will discuss the biogenesis of microRNAs and secondary RNAs involved in this silencing mechanism, biochemical and structural features of NBS-LRR proteins in response to pathogen effectors and the evolution of microRNA-based silencing mechanism with a focus on the miR482/2118 family. Furthermore, the biotechnological manipulation of microRNA expression, using both transgenic or genome editing approaches to improve cultivated plants will be discussed, with a focus on the miR482/2118 family in soybean.

Improving plant miRNA-target prediction with self-supervised k-mer embedding and spectral graph convolutional neural network

Deciphering the targets of microRNAs (miRNAs) in plants is crucial for comprehending their function and the variation in phenotype that they cause. As the highly cell-specific nature of miRNA regulation, recent computational approaches usually utilize expression data to identify the most physiologically relevant targets. Although these methods are effective, they typically require a large sample size and high-depth sequencing to detect potential miRNA-target pairs, thereby limiting their applicability in improving plant breeding. In this study, we propose a novel miRNA-target prediction framework named kmerPMTF (k-mer-based prediction framework for plant miRNA-target). Our framework effectively extracts the latent semantic embeddings of sequences by utilizing k-mer splitting and a deep self-supervised neural network. We construct multiple similarity networks based on k-mer embeddings and employ graph convolutional networks to derive deep representations of miRNAs and targets and calculate the probabilities of potential associations. We evaluated the performance of kmerPMTF on four typical plant datasets: Arabidopsis thaliana, Oryza sativa, Solanum lycopersicum, and Prunus persica. The results demonstrate its ability to achieve AUPRC values of 84.9%, 91.0%, 80.1%, and 82.1% in 5-fold cross-validation, respectively. Compared with several state-of-the-art existing methods, our framework achieves better performance on threshold-independent evaluation metrics. Overall, our study provides an efficient and simplified methodology for identifying plant miRNA-target associations, which will contribute to a deeper comprehension of miRNA regulatory mechanisms in plants.

Convergent and divergent evolution of microRNA-mediated regulation in metazoans

The evolution of microRNAs (miRNAs) has been studied extensively to understand their roles in gene regulation and evolutionary processes. This review focuses on how miRNA-mediated regulation has evolved in bilaterian animals, highlighting both convergent and divergent evolution. Since animals and plants display significant differences in miRNA biogenesis and target recognition, the 'independent origin' hypothesis proposes that miRNA pathways in these groups independently evolved from the RNA interference (RNAi) pathway, leading to modern miRNA repertoires through convergent evolution. However, recent evidence raises the alternative possibility that the miRNA pathway might have already existed in the last common ancestor of eukaryotes, and that the differences in miRNA pathway and miRNA repertoires among animal and plant lineages arise from lineage-specific innovations and losses of miRNA pathways, miRNA acquisition, and loss of miRNAs after eukaryotic divergence. The repertoire of miRNAs has considerably expanded during bilaterian evolution, primarily through de novo creation and duplication processes, generating new miRNAs. Although ancient functionally established miRNAs are rarely lost, many newly emerged miRNAs are transient and lineage specific, following a birth-death evolutionary pattern aligning with the 'out-of-the-testis' and 'transcriptional control' hypotheses. Our focus then shifts to the convergent molecular evolution of miRNAs. We summarize how miRNA clustering and seed mimicry contribute to this phenomenon, and we review how miRNAs from different sources converge to degrade maternal messenger RNAs (mRNAs) during animal development. Additionally, we describe how miRNAs evolve across species due to changes in sequence, seed shifting, arm switching, and spatiotemporal expression patterns, which can result in variations in target sites among orthologous miRNAs across distant strains or species. We also provide a summary of the current understanding regarding how the target sites of orthologous miRNAs can vary across strains or distantly related species. Although many paralogous miRNAs retain their seed or mature sequences after duplication, alterations can occur in the seed or mature sequences or expression patterns of paralogous miRNAs, leading to functional diversification. We discuss our current understanding of the functional divergence between duplicated miRNAs, and illustrate how the functional diversification of duplicated miRNAs impacts target site evolution. By investigating these topics, we aim to enhance our current understanding of the functions and evolutionary dynamics of miRNAs. Additionally, we shed light on the existing challenges in miRNA evolutionary studies, particularly the complexity of deciphering the role of miRNA-mediated regulatory network evolution in shaping gene expression divergence and phenotypic differences among species.

OsmiR319-OsPCF5 modulate resistance to brown planthopper in rice through association with MYB proteins

BACKGROUND

The brown planthopper (BPH) is a kind of piercing-sucking insect specific to rice, with the damage tops the list of pathogens and insects in recent years. microRNAs (miRNAs) are pivotal regulators of plant-environment interactions, while the mechanism underlying their function against insects is largely unknown.

RESULTS

Here, we confirmed that OsmiR319, an ancient and conserved miRNA, negatively regulated resistance to BPHs, with overexpression of OsmiR319 susceptible to BPH, while suppression of OsmiR319 resistant to BPH in comparison with wild type. Meanwhile, we identified several targets of OsmiR319 that may mediate BPH resistance. Among them, OsPCF5 was the most obviously induced by BPH feeding, and over expression of OsPCF5 was resistance to BPH. In addition, various biochemical assays verified that OsPCF5 interacted with several MYB proteins, such as OsMYB22, OsMYB30, and OsMYB30C.Genetically, we revealed that both OsMYB22 and OsMYB30C positively regulated BPH resistance. Genetic interaction analyses confirmed that OsMYB22 and OsMYB30C both function in the same genetic pathway with OsmiR319b to mediate BPH resistance.

CONCLUSIONS

Altogether, we revealed that OsPCF5 regulates BPH resistance via association with several MYB proteins downstream of OsmiR319, these MYB proteins might function as regulators of BPH resistance through regulating the phenylpropane synthesis.

Dicer-like2b suppresses the wiry leaf phenotype in tomato induced by tobacco mosaic virus

Dicer-like (DCL) proteins are principal components of RNA silencing, a major defense mechanism against plant virus infections. However, their functions in suppressing virus-induced disease phenotypes remain largely unknown. Here, we identified a role for tomato (Solanum lycopersicum) DCL2b in regulating the wiry leaf phenotype during defense against tobacco mosaic virus (TMV). Knocking out SlyDCL2b promoted TMV accumulation in the leaf primordium, resulting in a wiry phenotype in distal leaves. Biochemical and bioinformatics analyses showed that 22-nt virus-derived small interfering RNAs (vsiRNAs) accumulated less abundantly in slydcl2b mutants than in wild-type plants, suggesting that SlyDCL2b-dependent 22-nt vsiRNAs are required to exclude virus from leaf primordia. Moreover, the wiry leaf phenotype was accompanied by upregulation of Auxin Response Factors (ARFs), resulting from a reduction in trans-acting siRNAs targeting ARFs (tasiARFs) in TMV-infected slydcl2b mutants. Loss of tasiARF production in the slydcl2b mutant was in turn caused by inhibition of miRNA390b function. Importantly, silencing SlyARF3 and SlyARF4 largely restored the wiry phenotype in TMV-infected slydcl2b mutants. Our work exemplifies the complex relationship between RNA viruses and the endogenous RNA silencing machinery, whereby SlyDCL2b protects the normal development of newly emerging organs by excluding virus from these regions and thus maintaining developmental silencing.

The miR408a-BBP-LAC3/CSD1 module regulates anthocyanin biosynthesis mediated by crosstalk between copper homeostasis and ROS homeostasis during light induction in Malus plants

INTRODUCTION

The expression of miR408 is affected by copper (Cu) conditions and positively regulates anthocyanin biosynthesis in Arabidopsis. However, the underlying mechanisms by which miR408 regulates anthocyanin biosynthesis mediated by Cu homeostasis and reactive oxygen species (ROS) homeostasis remain unclear in Malus plants.

OBJECTIVES

Our study aims to elucidate how miR408a and its target, basic blue protein (BBP) regulate Cu homeostasis and ROS homeostasis, and anthocyanin biosynthesis in Malus plants.

METHODS

The roles of miR408a and its target BBP in regulating anthocyanin biosynthesis, Cu homeostasis, and ROS homeostasis were mainly identified in Malus plants.

RESULTS

We found that the BBP protein interacted with the copper-binding proteins LAC3 (laccase) and CSD1 (Cu/Zn SOD superoxide dismutase), indicating a potential crosstalk between Cu homeostasis and ROS homeostasis might be mediated by miR408 to regulate the anthocyanin accumulation. Further studies showed that overexpressing miR408a or suppressing BBP transiently significantly increased the expression of genes related to Cu binding and Cu transport, leading to anthocyanin accumulation under light induction in apple fruit and Malus plantlets. Consistently, opposite results were obtained when repressing miR408a or overexpressing BBP. Moreover, light induction significantly increased the expression of miR408a, CSD1, and LAC3, but significantly reduced the BBP expression, resulting in increased Cu content and anthocyanin accumulation. Furthermore, excessive Cu significantly increased the anthocyanin accumulation, accompanied by reduced expression of miR408a and Cu transport genes, and upregulated expression of Cu binding proteins including BBP, LAC3, and CSD1 to maintain the Cu homeostasis and ROS homeostasis in Malus plantlets.

CONCLUSION

Our findings provide new insights into the mechanism by which the miR408a-BBP-LAC3/CSD1 module perceives light and Cu signals regulating Cu and ROS homeostasis, ultimately affecting anthocyanin biosynthesis in Malus plants.

Manipulating microRNA miR408 enhances both biomass yield and saccharification efficiency in poplar

The conversion of lignocellulosic feedstocks to fermentable sugar for biofuel production is inefficient, and most strategies to enhance efficiency directly target lignin biosynthesis, with associated negative growth impacts. Here we demonstrate, for both laboratory- and field-grown plants, that expression of Pag-miR408 in poplar (Populus alba × P. glandulosa) significantly enhances saccharification, with no requirement for acid-pretreatment, while promoting plant growth. The overexpression plants show increased accessibility of cell walls to cellulase and scaffoldin cellulose-binding modules. Conversely, Pag-miR408 loss-of-function poplar shows decreased cell wall accessibility. Overexpression of Pag-miR408 targets three Pag-LACCASES, delays lignification, and modestly reduces lignin content, S/G ratio and degree of lignin polymerization. Meanwhile, the LACCASE loss of function mutants exhibit significantly increased growth and cell wall accessibility in xylem. Our study shows how Pag-miR408 regulates lignification and secondary growth, and suggest an effective approach towards enhancing biomass yield and saccharification efficiency in a major bioenergy crop.

microRNA408 and its encoded peptide regulate sulfur assimilation and arsenic stress response in Arabidopsis

MicroRNAs (miRNAs) are small non-coding RNAs that play a central role in regulating various developmental and biological processes. The expression of miRNAs is differentially modulated in response to various biotic and abiotic stresses. Recent findings have shown that some pri-miRNAs encode small regulatory peptides known as microRNA-encoded peptides (miPEPs). miPEPs regulate the growth and development of plants by modulating corresponding miRNA expression; however, the role of these peptides under different stress conditions remains unexplored. Here, we report that pri-miR408 encodes a small peptide, miPEP408, that regulates the expression of miR408, its targets, and associated phenotype in Arabidopsis. We also report that miR408, apart from Plantacyanin (ARPN) and Laccase3 (LAC3), targets a glutathione S-transferase (GSTU25) that plays a role in sulfur assimilation and exhibits a range of detoxification activities with the environmental pollutant. Plants overexpressing miR408 showed severe sensitivity under low sulfur (LS), arsenite As(III), and LS + As(III) stress, while miR408 mutants developed using the CRISPR/Cas9 approach showed tolerance. Transgenic lines showed phenotypic alteration and modulation in the expression of genes involved in the sulfur reduction pathway and affect sulfate and glutathione accumulation. Similar to miR408 overexpressing lines, the exogenous application of synthetic miPEP408 and miPEP408OX lines led to sensitivity in plants under LS, As(III), and combined LS + As(III) stress compared to the control. This study suggests the involvement of miR408 and miPEP408 in heavy metal and nutrient deficiency responses through modulation of the sulfur assimilation pathway.

Plant non-coding RNAs function in pollen development and male sterility

Male sterility is classified as either cytoplasmic male sterility (CMS) or genic male sterility (GMS). Generally, CMS involves mitochondrial genomes interacting with the nuclear genome, while GMS is caused by nuclear genes alone. Male sterility is regulated by multilevel mechanisms in which non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and phased small interfering RNAs (phasiRNAs), which have been proven to be critical elements. The development of high-throughput sequencing technology offers new opportunities to evaluate the genetic mechanism of ncRNAs in plant male sterility. In this review, we summarize the critical ncRNAs that regulate gene expression in ways dependent on or independent of hormones, which involve the differentiation of the stamen primordia, degradation of the tapetum, formation of microspores, and the release of pollen. In addition, the key mechanisms of the miRNA-lncRNA-mRNA interaction networks mediating male sterility in plants are elaborated. We present a different perspective on exploring the ncRNA-mediated regulatory pathways that control CMS in plants and create male-sterile lines through hormones or genome editing. A refined understanding of the ncRNA regulatory mechanisms in plant male sterility for the development of new sterile lines would be conducive to improve hybridization breeding.

Convergent evolution of the annual life history syndrome from perennial ancestors

Despite most angiosperms being perennial, once-flowering annuals have evolved multiple times independently, making life history traits among the most labile trait syndromes in flowering plants. Much research has focused on discerning the adaptive forces driving the evolution of annual species, and in pinpointing traits that distinguish them from perennials. By contrast, little is known about how 'annual traits' evolve, and whether the same traits and genes have evolved in parallel to affect independent origins of the annual syndrome. Here, we review what is known about the distribution of annuals in both phylogenetic and environmental space and assess the evidence for parallel evolution of annuality through similar physiological, developmental, and/or genetic mechanisms. We then use temperate grasses as a case study for modeling the evolution of annuality and suggest future directions for understanding annual-perennial transitions in other groups of plants. Understanding how convergent life history traits evolve can help predict species responses to climate change and allows transfer of knowledge between model and agriculturally important species.

A novel tasi RNA-based micro RNA-induced gene silencing strategy to tackle multiple pests and pathogens in cotton (Gossypium hirsutum L.)

This study demonstrates the combinatorial management of multiple pests through a trans-acting siRNA (tasiRNA)-based micro RNA-induced gene silencing (MIGS) strategy. Transgenic cotton events demonstrated improved efficacy against cotton leaf curl disease, cotton leaf hopper and root-knot nematode. Cotton (Gossypium hirsutum L.), an important commercial crop grown worldwide is confronted by several pests and pathogens, thus reiterating interventions for their management. In this study, we report, the utility of a novel Arabidopsis miRNA173-directed trans-acting siRNA (tasiRNA)-based micro RNA-induced gene silencing (MIGS) strategy for the simultaneous management of cotton leaf curl disease (CLCuD), cotton leaf hopper (CLH; Amrasca biguttula biguttula) and root-knot nematode (RKN, Meloidogyne incognita). Cotton transgenics were developed with the MIGS construct targeting a total of 7 genes by an apical meristem-targeted in planta transformation strategy. Stable transgenics were selected using stringent selection pressure, molecular characterization and stress-specific bio-efficacy studies. We identified 8 superior events with 50-100% resistance against CLCuD, while reduction in the root-knot nematode multiplication factor in the range of 35-75% confirmed resistance to RKN. These transgenic cotton events were also detrimental to the growth and development of CLH, as only 43.3-62.5% of nymphs could survive. Based on the corroborating evidences obtained by all the bioefficacy analyses, 3 events viz., L-75-1, E-27-11, E-27-7 were found to be consistent in tackling the target pests. To the best of our knowledge, this report is the first of its kind demonstrating the possibility of combinatorial management of pests/diseases in cotton using MIGS approach. These identified events demonstrate immense utility of the strategy towards combinatorial stress management in cotton improvement programs.

Analysis of Homologous Regions of Small RNAs MIR397 and MIR408 Reveals the Conservation of Microsynteny among Rice Crop-Wild Relatives

MIRNAs are small non-coding RNAs that play important roles in a wide range of biological processes in plant growth and development. MIR397 (involved in drought, low temperature, and nitrogen and copper (Cu) starvation) and MIR408 (differentially expressed in response to environmental stresses such as copper, light, mechanical stress, dehydration, cold, reactive oxygen species, and drought) belong to conserved MIRNA families that either negatively or positively regulate their target genes. In the present study, we identified the homologs of MIR397 and MIR408 in Oryza sativa and its six wild progenitors, three non-Oryza species, and one dicot species. We analyzed the 100 kb segments harboring MIRNA homologs from 11 genomes to obtain a comprehensive view of their community evolution around these loci in the farthest (distant) relatives of rice. Our study showed that mature MIR397 and MIR408 were highly conserved among all Oryza species. Comparative genomics analyses also revealed that the microsynteny of the 100 kb region surrounding MIRNAs was only conserved in Oryza spp.; disrupted in Sorghum, maize, and wheat; and completely lost in Arabidopsis. There were deletions, rearrangements, and translocations within the 100 kb segments in Oryza spp., but the overall microsynteny of the region was maintained. The phylogenetic analyses of the precursor regions of all MIRNAs under study revealed a bimodal clade of common origin. This comparative analysis of miRNA involved in abiotic stress tolerance in plants provides a powerful tool for future Oryza research. Crop wild relatives (CWRs) offer multiple traits with potential to decrease the amount of yield loss owing to biotic and abiotic stresses. Using a comparative genomics approach, the exploration of CWRs as a source of tolerance to these stresses by understanding their evolution can be further used to leverage their yield potential.

Integrated miRNA, target mRNA, and metabolome profiling of Tinospora cordifolia with reference to berberine biosynthesis

MicroRNAs play a central role in gene regulation and emerge as novel targets for secondary metabolites improvement in plants. The crops thus can be improved through knowledge obtained by the study of miRNAs because of their conserved nature in gene regulation. The present study has been carried out on Tinospora cordifolia (T. cordifolia), because of its illimitable application for the treatment of various diseases. This plant has tremendous medicinal properties, yet unexplored at the molecular level, and has not received much recognition in the scientific field. Thus, here computational analysis was performed to identify T. cordifolia miRNAs using EST database. Using these miRNAs, we predicted their targets which were found to be associated with the regulation of diverse gene networks including 433 berberine biosynthesis genes in T. cordifolia. Further, selected miRNAs were validated and their expression was detected in different T. cordifolia tissues followed by expression analysis of their target mRNAs. These data were then compared with the metabolic profile of T. cordifolia with an emphasis on therapeutically important compound berberine. In this study, we did simultaneous miRNA/target gene expression and metabolome analysis which opens a new way for initiating new proposition and prioritization of miRNAs/genes/metabolites for targeted follow‑up metabolic engineering experimentations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03342-9.

Tomato MicroRNAs and Their Functions

MicroRNAs (miRNAs) define an essential class of non-coding small RNAs that function as posttranscriptional modulators of gene expression. They are coded by MIR genes, several hundreds of which exist in the genomes of Arabidopsis and rice model plants. The functional analysis of Arabidopsis and rice miRNAs indicate that their miRNAs regulate a wide range of processes including development, reproduction, metabolism, and stress. Tomato serves as a major model crop for the study of fleshy fruit development and ripening but until recently, information on the identity of its MIR genes and their coded miRNAs was limited and occasionally contradictory. As a result, the majority of tomato miRNAs remained uncharacterized. Recently, a comprehensive annotation of tomato MIR genes has been carried out by several labs and us. In this review, we curate and organize the resulting partially overlapping MIR annotations into an exhaustive and non-redundant atlas of tomato MIR genes. There are 538 candidate and validated MIR genes in the atlas, of which, 169, 18, and 351 code for highly conserved, Solanaceae-specific, and tomato-specific miRNAs, respectively. Furthermore, a critical review of functional studies on tomato miRNAs is presented, highlighting validated and possible functions, creating a useful resource for future tomato miRNA research.

Genome-wide identification of miRNAs and targets associated with cell wall biosynthesis: Differential roles of dlo-miR397a and dlo-miR408-3p during early somatic embryogenesis in longan

The dynamic alterations in cell wall (CW) biosynthesis play an essential role in physiological isolation during the plant somatic embryogenesis (SE). However, the mechanisms underlying the functions of cell wall-associated miRNAs (CW-miRNA) remain poorly understood in plant SE. Here, we have identified 36 distinct candidate miRNAs associated with CW biosynthesis from longan third-generation genome as well as miRNA transcriptome, and modified RLM-RACE validated four distinct miRNA, which specifically targeted four CW-related genes. More importantly, we found that the dlo-miR397a-antagomir significantly enhanced DlLAC7 expression and improved laccase activity. Interestingly, inhibition of dlo-miR397a increased CW lignin deposition and promoted the tightening of protodermal cell by miRNA-mimic technology during early SE. Moreover, overexpression of dlo-miR408-3p (dlo-miR408-3p-agomir) markedly decreased DlLAC12 expression. dlo-miR408-3p-agomir activated rapid cell division, thus promoting the globular embryo (GE) development, which might be due to high DNA synthesis activity in protoepidermal cells, rather than affecting lignin synthesis. The subcellular location also indicated that both DlLAC7 and DlLAC12 proteins were primarily localized in CW and regulated CW biosynthesis. Overall, our findings provided new insight on the molecular regulatory networks comprising various miRNAs associated with cell wall, and established that dlo-miR397a and dlo-miR408-3p played differential roles during early SE in longan. The findings also shed some light on the potential role of miRNA target DlLAC regulating in vivo embryonic development of plant.

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.

The renaissance and enlightenment of Marchantia as a model system

The liverwort Marchantia polymorpha has been utilized as a model for biological studies since the 18th century. In the past few decades, there has been a Renaissance in its utilization in genomic and genetic approaches to investigating physiological, developmental, and evolutionary aspects of land plant biology. The reasons for its adoption are similar to those of other genetic models, e.g. simple cultivation, ready access via its worldwide distribution, ease of crossing, facile genetics, and more recently, efficient transformation, genome editing, and genomic resources. The haploid gametophyte dominant life cycle of M. polymorpha is conducive to forward genetic approaches. The lack of ancient whole-genome duplications within liverworts facilitates reverse genetic approaches, and possibly related to this genomic stability, liverworts possess sex chromosomes that evolved in the ancestral liverwort. As a representative of one of the three bryophyte lineages, its phylogenetic position allows comparative approaches to provide insights into ancestral land plants. Given the karyotype and genome stability within liverworts, the resources developed for M. polymorpha have facilitated the development of related species as models for biological processes lacking in M. polymorpha.

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

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

Evaluation of the Possible Contribution of Various Regulatory Genes to Determination of Carpel Number as a Potential Mechanism for Optimal Agricultural Yield

Various regulatory genes encoding transcription factors and miRNAs regulate carpel number. Multicarpelly is normally associated with increased size of the floral meristem, and several genetic factors have been discovered that influence this characteristic. A fundamental understanding of the regulatory genes affecting carpel number can facilitate strategies for agricultural yield improvement, which is crucial, given that the global population is growing rapidly. A multicarpellate plant may provide a significantly higher yield than a plant bearing fewer carpels. Higher yields can be achieved via various means; in this review, we provide an overview of the current knowledge of the various regulatory factors that contribute to multicarpelly and the potential of increasing carpel number to achieve an increased yield.

Evolutionary conservation of nested MIR159 structural microRNA genes and their promoter characterization in Arabidopsis thaliana

MicroRNAs (miRNAs) are endogenous small RNAs, that are vital for gene expression regulation in eukaryotes. Whenever a pri-miRNA precursor includes another miRNA precursor, and both of these precursors may generate independent, non-overlapping mature miRNAs, we named them nested miRNAs. However, the extent of nested miR159 structural evolutionary conservation and its promoter characterization remains unknown. In this study, the sequence alignment and phylogenetic analysis reveal that the MIR159 family is ancient, and its nested miR159 structures are evolutionary conserved in different plant species. The overexpression of ath-MIR159a, including the 1.2 kb downstream region, has no effect on rescuing the mir159ab phenotype. The promoter truncation results revealed that the 1.0 kb promoter of ath-MIR159a is sufficient for rescuing the mir159ab phenotype. The cis-regulatory elements in the ath-miR159a promoters indicated functions related to different phytohormones, abiotic stresses, and transcriptional activation. While the MybSt1 motif-containing region is not responsible for activating the regulation of the miR159a promoter. The qRT-PCR results showed that overexpression of ath-MIR159a led to high expression levels of miR159a.1-5 and miR159a.1-3 and complemented the growth defect of mir159ab via downregulation of MYB33 and MYB65. Furthermore, continuously higher expression of the miR159a.2 duplex in transgenic lines with the curly leaf phenotype indicates that miR159a.2 is functional in Arabidopsis and suggests that it is possible for a miRNA precursor to encode several regulatory small RNAs in plants. Taken together, our study demonstrates that the nested miR159 structure is evolutionary conserved and miRNA-mediated gene regulation is more complex than previously thought.

Differential Responses of Wheat (Triticum aestivum L.) and Cotton (Gossypium hirsutum L.) to Nitrogen Deficiency in the Root Morpho-Physiological Characteristics and Potential MicroRNA-Mediated Mechanisms

Understanding the mechanism of crop response to nitrogen (N) deficiency is very important for developing sustainable agriculture. In addition, it is unclear if the microRNA-mediated mechanism related to root growth complies with a common mechanism in monocots and dicots under N deficiency. Therefore, the root morpho-physiological characteristics and microRNA-mediated mechanisms were studied under N deficiency in wheat (Triticum aestivum L.) and cotton (Gossypium hirsutum L.). For both crops, shoot dry weight, plant dry weight and total leaf area as well as some physiological traits, i.e., the oxygen consuming rate in leaf and root, the performance index based on light energy absorption were significantly decreased after 8 days of N deficiency. Although N deficiency did not significantly impact the root biomass, an obvious change on the root morphological traits was observed in both wheat and cotton. After 8 days of treatment with N deficiency, the total root length, root surface area, root volume of both crops showed an opposite trend with significantly decreasing in wheat but significantly increasing in cotton, while the lateral root density was significantly increased in wheat but significantly decreased in cotton. At the same time, the seminal root length in wheat and the primary root length in cotton were increased after 8 days of N deficiency treatment. Additionally, the two crops had different root regulatory mechanisms of microRNAs (miRNAs) to N deficiency. In wheat, the expressions of miR167, miR319, miR390, miR827, miR847, and miR165/166 were induced by N treatment; these miRNAs inhibited the total root growth but promoted the seminal roots growth and lateral root formation to tolerate N deficiency. In cotton, the expressions of miR156, miR167, miR171, miR172, miR390, miR396 were induced and the expressions of miR162 and miR393 were inhibited; which contributed to increasing in the total root length and primary root growth and to decreasing in the lateral root formation to adapt the N deficiency. In conclusion, N deficiency significantly affected the morpho-physiological characteristics of roots that were regulated by miRNAs, but the miRNA-mediated mechanisms were different in wheat and cotton.

The genetic basis of natural variation in the timing of vegetative phase change in Arabidopsis thaliana

The juvenile-to-adult transition in plants is known as vegetative phase change and is marked by changes in the expression of leaf traits in response to a decrease in the level of miR156 and miR157. To determine whether this is the only mechanism of vegetative phase change, we measured the appearance of phase-specific leaf traits in 70 natural accessions of Arabidopsis thaliana. We found that leaf shape was poorly correlated with abaxial trichome production (two adult traits), that variation in these traits was not necessarily correlated with the level of miR156, and that there was little to no correlation between the appearance of adult-specific vegetative traits and flowering time. We identified eight quantitative trait loci controlling phase-specific vegetative traits from a cross between the Columbia (Col-0) and Shakdara (Sha) accessions. Only one of these quantitative trait loci includes genes known to regulate vegetative phase change (MIR156A and TOE1), which were expressed at levels consistent with the precocious phenotype of Sha. Our results suggest that vegetative phase change is regulated both by the miR156/SPL module and by genes specific to different vegetative traits, and that natural variation in vegetative phase change can arise from either source.

Verticillium dahliae Secretes Small RNA to Target Host MIR157d and Retard Plant Floral Transition During Infection

Bidirectional trans-kingdom RNA silencing [or RNA interference (RNAi)] plays a key role in plant-pathogen interactions. It has been shown that plant hosts export specific endogenous miRNAs into pathogens to inhibit their virulence, whereas pathogens deliver small RNAs (sRNAs) into plant cells to disturb host immunity. Here, we report a trans-kingdom fungal sRNA retarding host plant floral transition by targeting a miRNA precursor. From Arabidopsis plants infected with Verticillium dahliae, a soil-borne hemibiotrophic pathogenic fungus that causes wilt diseases in a wide range of plant hosts, we obtained a number of possible trans-kingdom V. dahliae sRNAs (VdsRNAs) by sequencing AGO1-immunoprecipitated sRNAs. Among these, a 24-nt VdsRNA derived from V. dahliae rRNA, VdrsR-1, was shown to be an actual trans-kingdom VdsRNA that targets the miR157d precursor MIR157d, resulting in increased rather than reduced miR157d accumulation in V. dahliae-infected plants. Consistent with the miR157 family in the regulation of vegetative and floral transitions by targeting SPL genes in several plant species, we detected two SPL genes, SPL13A/B, that were notably reduced in V. dahliae-infected and VdrsR-1-expressing plants compared with control plants. Furthermore, V. dahliae-infected and VdrsR-1-expressing plants also displayed delayed vegetative phase change and floral transition compared to control plants. Taken together, we disclosed a novel mode of action for a trans-kingdom fungal sRNA, VdrsR-1, which was secreted into host cells to modulate plant floral transition by employing the miR157d/SPL13A/B regulatory module, leading to prolonged host vegetative growth that would undoubtedly benefit fungal propagation.

Advances in the regulation of plant salt-stress tolerance by miRNA

Salt stress significantly affects the growth, development, yield, and quality of plants. MicroRNAs (miRNAs) are involved in various stress responses via target gene regulation. Their role in regulating salt stress has also received significant attention from researchers. Various transcription factor families are the common target genes of plant miRNAs. Thus, regulating the expression of miRNAs is a novel method for developing salt-tolerant crops. This review summarizes plant miRNAs that mediate salt tolerance, specifically miRNAs that have been utilized in genetic engineering to modify plant salinity tolerance. The molecular mechanism by which miRNAs mediate salt stress tolerance merits elucidation, and this knowledge will promote the development of miRNA-mediated salt-tolerant crops and provide new strategies against increasingly severe soil salinization.

An Integration of MicroRNA and Transcriptome Sequencing Analysis Reveal Regulatory Roles of miRNAs in Response to Chilling Stress in Wild Rice

A chromosome single segment substitution line (CSSL) DC90, which was generated by introgressing CTS-12, a locus derived from common wild rice (Oryza rufipogon Griff.), into the 9311 (Oryza sativa L. ssp. indica) background, exhibits a chilling tolerance phenotype under chilling stress. Here, an integration of microRNA (miRNA) deep sequencing and transcriptomic sequencing analysis was performed to explore the expression profiles of miRNAs and their target genes mediated by CTS-12 under chilling stress, and to reveal the possible regulatory mechanisms of miRNAs that are involved in chilling tolerance. Integration analysis revealed that a number of differentially expressed miRNAs (DEMs) and putative target genes with different expression patterns and levels were identified in 9311 and DC90 under chilling stress. KEGG enrichment analysis revealed that the target genes that are regulated by chilling-induced miRNAs are involved in the regulation of various biological processes/pathways, including protein biosynthesis, redox process, photosynthetic process, and chloroplast development in two genotypes. CRISPR/Cas9 editing of the target genes of the key DEMs in a chilling tolerant rice variety Zhonghua 11 (ZH11) found that LOC_Os11g48020 (OsGL1-11), one of the putative target genes of osa-miR1846a/b-5p and encoding a wax synthesis protein, is correlated with a chilling stress tolerance phenotype, implying osa-miR1846a/b-5p/OsGL1-11 plays an important role in CTS-12-mediated chilling stress tolerance regulatory pathway(s). Therefore, we speculate that the CTS-12 may regulate the key miRNA target genes in response to chilling stress by differential regulation of miRNAs in wild rice, thereby resulting in the variation of chilling tolerance phenotype between 9311 and DC90.

The carbon economics of vegetative phase change

Across plant species and biomes, a conserved set of leaf traits govern the economic strategy used to assimilate and invest carbon. As plants age, they face new challenges that may require shifts in this leaf economic strategy. In this study, we investigate the role of the developmental transition, vegetative phase change (VPC), in altering carbon economics as plants age. We used overexpression of microRNA 156 (miR156), the master regulator of VPC, to modulate the timing of VPC in Populus tremula x alba, Arabidopsis thaliana and Zea mays to understand the impact of this transition on leaf economic traits, including construction cost, payback time and return on investment. Here, we find that VPC causes a shift from a low-cost, quick return juvenile strategy to a high-cost, high-return adult strategy. The juvenile strategy is advantageous in light-limited conditions, whereas the adult strategy provides greater returns in high light. The transition between these strategies is correlated with the developmental decline in the level of miR156, suggesting that is regulated by the miR156/SPL pathway. Our results provide an ecophysiological explanation for the existence of juvenile and adult leaf types and suggest that natural selection for these alternative economic strategies could be an important factor in plant evolution.

Rice microRNA171f/SCL6 module enhances drought tolerance by regulation of flavonoid biosynthesis genes

Plants have evolved sophisticated defense systems to enhance drought tolerance. These include the microRNA (miRNA) group of small noncoding RNAs that act as post-transcriptional regulators; however, details of the mechanisms by which they confer drought tolerance are not well understood. Here, we show that osa-MIR171f, a member of osa-MIR171 gene family, is mainly expressed in response to drought stress and regulates the transcript levels of SCARECROW-LIKE6-I (SCL6-I) and SCL6-II in rice (Oryza sativa). The SCL6 genes are known to be involved in shoot branching and flag leaf morphology. Osa-MIR171f-overexpressing (osa-MIR171f-OE) transgenic plants showed reduced drought symptoms compared with non-transgenic (NT) control plants under both field drought and polyethylene glycol (PEG)-mediated dehydration stress conditions. Transcriptome analysis of osa-MIR171f-OE plants and osa-mir171f-knockout (K/O) lines generated by clustered regularly interspaced short palindromic repeats (CRISPR/Cas9) revealed that osa-mature-miR171a-f (osa-miR171) regulates the expression of flavonoid biosynthesis genes, consequently leading to drought tolerance. This upregulation in the osa-MIR171f-OE plants, which did not occur in NT control plants, was observed under both normal and drought conditions. Our findings indicate that osa-miR171 plays a role in drought tolerance by regulating SCL6-I and SCL6-II transcript levels.

Plant RNA-mediated gene regulatory network

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

CRISPR/Cas9 mutants of tomato MICRORNA164 genes uncover their functional specialization in development

Plant MICRORNA164 (miR164) plays diverse regulatory functions by post-transcriptional repression of certain NAM/ATAF/CUC-domain transcription factors. However, the involvement of miR164 in fleshy fruit development and ripening remains poorly understood. Here, de novo prediction of tomato (Solanum lycopersicum) MIR164 genes identified four genes (SlMIR164a-d), of which SlMIR164d has an atypically long pre-miRNA. The roles of the fruit expressed SlMIR164a, b, and d were studied by analysis of their Clustered Regularly Interspaced Short Palindromic Repeats mutants. The slmir164bCR mutant plants exhibited shoot and flower abnormalities characteristic of ectopic boundary specification, whereas the shoot and flower development of slmir164aCR and slmir164dCR mutants were indistinguishable from wild-type. Strikingly, the knockout of SlMIR164a practically eliminated sly-miR164 from the developing and ripening fruit pericarp. The sly-miR164-deficient slmir164aCR fruits were smaller than the wild-type, due to reduced pericarp cell division and expansion, and displayed intense red color and matte, instead of glossy appearance, upon ripening. We found that the fruit skin phenotypes were associated with morphologically abnormal outer epidermis and thicker cuticle. Quantitation of sly-miR164 target transcripts in slmir164aCR ripening fruits demonstrated the upregulation of SlNAM3 and SlNAM2. Specific expression of their miR164-resistant versions in the pericarp resulted in the formation of extremely small fruits with abnormal epidermis, highlighting the importance of their negative regulation by sly-miR164a. Taken together, our results demonstrate that SlMIR164a and SlMIR164b play specialized roles in development: SlMIR164b is required for shoot and flower boundary specification, and SlMIR164a is required for fruit growth including the expansion of its outer epidermis, which determines the properties of the fruit skin.

PpGRAS12 acts as a positive regulator of meristem formation in Physcomitrium patens

This study focused on the key regulatory function of Physcomitrium patens GRAS12 gene underlying an increasing plant complexity, an important step in plant terrestrialization and the evolutionary history of life. The miR171-GRAS module has been identified as a key player in meristem maintenance in angiosperms. PpGRAS12 is a member of the GRAS family and a validated target for miR171 in Physcomitrium (Physcomitrella) patens. Here we show a regulatory function of miR171 at the gametophytic vegetative growth stage and targeted deletion of the PpGRAS12 gene adversely affects sporophyte production since fewer sporophytes were produced in ΔPpGRAS12 knockout lines compared to wild type moss. Furthermore, highly specific and distinct growth arrests were observed in inducible PpGRAS12 overexpression lines at the protonema stage. Prominent phenotypic aberrations including the formation of multiple apical meristems at the gametophytic vegetative stage in response to elevated PpGRAS12 transcript levels were discovered via scanning electron microscopy. The production of multiple buds in the PpGRAS12 overexpression lines similar to ΔPpCLV1a/1b disruption mutants is accompanied by an upregulation of PpCLE and downregulation of PpCLV1, PpAPB, PpNOG1, PpDEK1, PpRPK2 suggesting that PpGRAS12 acts upstream of these genes and negatively regulates the proposed pathway to specify simplex meristem formation. As CLV signaling pathway components are not present in the chlorophytic or charophytic algae and arose with the earliest land plants, we identified a key regulatory function of PpGRAS12 underlying an increasing plant complexity, an important step in plant terrestrialization and the evolutionary history of life.

Cytokinin regulates vegetative phase change in Arabidopsis thaliana through the miR172/TOE1-TOE2 module

During vegetative growth plants pass from a juvenile to an adult phase causing changes in shoot morphology. This vegetative phase change is primarily regulated by the opposite actions of two microRNAs, the inhibitory miR156 and the promoting miR172 as well as their respective target genes, constituting the age pathway. Here we show that the phytohormone cytokinin promotes the juvenile-to-adult phase transition through regulating components of the age pathway. Reduction of cytokinin signalling substantially delayed the transition to the adult stage. tZ-type cytokinin was particularly important as compared to iP- and the inactive cZ-type cytokinin, and root-derived tZ influenced the phase transition significantly. Genetic and transcriptional analyses indicated the requirement of SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors and miR172 for cytokinin activity. Two miR172 targets, TARGET OF EAT1 (TOE1) and TOE2 encoding transcriptional repressors were necessary and sufficient to mediate the influence of cytokinin on vegetative phase change. This cytokinin pathway regulating plant aging adds to the complexity of the regulatory network controlling the juvenile-to-adult phase transition and links cytokinin to miRNA action.

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

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

The Copper-microRNA Pathway Is Integrated with Developmental and Environmental Stress Responses in Arabidopsis thaliana

As an essential nutrient, copper (Cu) scarcity causes a decrease in agricultural production. Cu deficiency responses include the induction of several microRNAs, known as Cu-miRNAs, which are responsible for degrading mRNAs from abundant and dispensable cuproproteins to economize copper when scarce. Cu-miRNAs, such as miR398 and miR408 are conserved, as well as the signal transduction pathway to induce them under Cu deficiency. The Arabidopsis thaliana SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) family member SPL7 binds to the cis-regulatory motifs present in the promoter regions of genes expressed under Cu deficiency, including Cu-miRNAs. The expression of several other SPL transcription factor family members is regulated by miR156. This regulatory miR156-SPL module plays a crucial role in developmental phase transitions while integrating internal and external cues. Here, we show that Cu deficiency also affects miR156 expression and that SPL3 overexpressing plants, resistant to miR156 regulation, show a severe decrease in SPL7-mediated Cu deficiency responses. These include the expression of Cu-miRNAs and their targets and is probably due to competition between SPL7 and miR156-regulated SPL3 in binding to cis-regulatory elements in Cu-miRNA promoters. Thus, the conserved SPL7-mediated Cu-miRNA pathway could generally be affected by the miR156-SPL module, thereby underscoring the integration of the Cu-miRNA pathway with developmental and environmental stress responses in Arabidopsis thaliana.

Knockdown of MicroRNA160a/b by STTM leads to root architecture changes via auxin signaling in Solanum tuberosum

The root phenotype is an important aspect of plant architecture and plays a critical role in plant facilitation of the extraction of water and nutrition from the soil. MicroRNAs (miRNAs) are classes of small RNAs with important roles in regulating endogenous gene expression at the post-transcriptional level that function in a range of plant development processes and in the response to abiotic stresses. However, little is known concerning the molecular mechanism of miRNAs in regulating the generation and development of plant root architecture. Herein, we demonstrated that potato miR160a/b acted as a critical regulator and affected plant root architecture by targeting the mRNA of StARF10 and StARF16 for cleavage. The miR160a/b precursor was cloned from potato. Quantitative PCR assays showed that the expression levels of miR160 and its targets were down- or up-regulated with the development of potato roots, respectively. Moreover, transgenic lines with suppressed stu-miR160 expression were established with the short tandem targets mimic (STTM), and the results showed that the ectopic expression of miR160a/b altered the levels of auxin and the expression of auxin signaling-related genes and caused drastic change in root architecture compared with that in control plants. Suppressing the expression of miR160 led to a severe reduction in root length, an increase in the number of lateral roots, and a decrease in fresh root weight in potato. Collectively, our data established a key role of miR160 in modulating plant root architecture in potato.

MicroRNA156-mediated changes in leaf composition lead to altered photosynthetic traits during vegetative phase change

Plant morphology and physiology change with growth and development. Some of these changes are due to change in plant size and some are the result of genetically programmed developmental transitions. In this study we investigate the role of the developmental transition, vegetative phase change (VPC), on morphological and photosynthetic changes. We used overexpression of microRNA156, the master regulator of VPC, to modulate the timing of VPC in Populus tremula × alba, Zea mays, and Arabidopsis thaliana to determine its role in trait variation independent of changes in size and overall age. Here, we find that juvenile and adult leaves in all three species photosynthesize at different rates and that these differences are due to phase-dependent changes in specific leaf area (SLA) and leaf N but not photosynthetic biochemistry. Further, we found juvenile leaves with high SLA were associated with better photosynthetic performance at low light levels. This study establishes a role for VPC in leaf composition and photosynthetic performance across diverse species and environments. Variation in leaf traits due to VPC are likely to provide distinct benefits under specific environments; as a result, selection on the timing of this transition could be a mechanism for environmental adaptation.

Comparative analysis of sRNAs, degradome and transcriptomics in sweet sorghum reveals the regulatory roles of miRNAs in Cd accumulation and tolerance

Key miRNAs including sbi-miR169p/q, sbi-miR171g/j, sbi-miR172a/c/d, sbi-miR172e, sbi-miR319a/b, sbi-miR396a/b, miR408, sbi-miR5384, sbi-miR5565e and nov_23 were identified to function in the regulation of Cd accumulation and tolerance. As an energy plant, sweet sorghum shows great potential in the phytoremediation of Cd-contaminated soils. However, few studies have focused on the regulatory roles of miRNAs and their targets under Cd stress. In this study, comparative analysis of sRNAs, degradome and transcriptomics was conducted in high-Cd accumulation (H18) and low-Cd accumulation (L69) genotypes of sweet sorghum. A total of 38 conserved and 23 novel miRNAs with differential expressions were identified under Cd stress or between H18 and L69, and 114 target genes of 41 miRNAs were validated. Furthermore, 25 miRNA-mRNA pairs exhibited negatively correlated expression profiles and sbi-miR172e together with its target might participate in the distinct Cd tolerance between H18 and L69 as well as sbi-miR172a/c/d. Additionally, two groups of them: miR169p/q-nov_23 and miR408 were focused through the co-expression analysis, which might be involved in Cd uptake and tolerance by regulating their targets associated with transmembrane transportation, cytoskeleton activity, cell wall construction and ROS (reactive oxygen species) homeostasis. Further experiments exhibited that cell wall components of H18 and L69 were different when exposed to cadmium, which might be regulated by miR169p/q, miR171g/j, miR319a/b, miR396a/b, miR5384 and miR5565e through their targets. Through this study, we aim to reveal the potential miRNAs involved in sweet sorghum in response to Cd stress and provide references for developing high-Cd accumulation or high Cd-resistant germplasm of sweet sorghum that can be used in phytoremediation.

Efficient deletion of multiple circle RNA loci by CRISPR-Cas9 reveals Os06circ02797 as a putative sponge for OsMIR408 in rice

CRISPR-Cas9 is an emerging genome editing tool for reverse genetics in plants. However, its application for functional study of non-coding RNAs in plants is still at its infancy. Despite being a major class of non-coding RNAs, the biological roles of circle RNAs (circRNAs) remain largely unknown in plants. Previous plant circRNA studies have focused on identification and annotation of putative circRNAs, with their functions largely uninvestigated by genetic approaches. Here, we applied a multiplexed CRISPR-Cas9 strategy to efficiently acquire individual null mutants for four circRNAs in rice. We showed each of these rice circRNA loci (Os02circ25329, Os06circ02797, Os03circ00204 and Os05circ02465) can be deleted at 10% or higher efficiency in both protoplasts and stable transgenic T0 lines. Such high efficiency deletion enabled the generation of circRNA null allele plants without the CRISPR-Cas9 transgene in the T1 generation. Characterization of the mutants reveals these circRNAs' participation in salt stress response during seed germination and in particular the Os05circ02465 null mutant showed high salt tolerance. Notably, the seedlings of the Os06circ02797 mutant showed rapid growth phenotype after seed germination with the seedlings containing higher chlorophyll A/B content. Further molecular and computational analyses suggested a circRNA-miRNA-mRNA regulatory network where Os06circ02797 functions to bind and sequester OsMIR408, an important and conserved microRNA in plants. This study not only presents genetic evidence for the first time in plants that certain circRNAs may serve as sponges to negatively regulate miRNAs, a phenomenon previously demonstrated in mammalian cells, but also provides important insights for improving agronomic traits through gene editing of circRNA loci in crops.

Hierarchical Action of Mulberry miR156 in the Vegetative Phase Transition

The vegetative phase transition is a prerequisite for flowering in angiosperm plants. Mulberry miR156 has been confirmed to be a crucial factor in the vegetative phase transition in Arabidopsis thaliana. The over-expression of miR156 in transgenic Populus × canadensis dramatically prolongs the juvenile phase. Here, we find that the expression of mno-miR156 decreases with age in all tissues in mulberry, which led us to study the hierarchical action of miR156 in mulberry. Utilizing degradome sequencing and dual-luciferase reporter assays, nine MnSPLs were shown to be directly regulated by miR156. The results of yeast one-hybrid and dual-luciferase reporter assays also revealed that six MnSPLs could recognize the promoter sequences of mno-miR172 and activate its expression. Our results demonstrate that mno-miR156 performs its role by repressing MnSPL/mno-miR172 pathway expression in mulberry. This work uncovered a miR156/SPLs/miR172 regulation pathway in the development of mulberry and fills a gap in our knowledge about the molecular mechanism of vegetative phase transition in perennial woody plants.

MIR172d Is Required for Floral Organ Identity and Number in Tomato

MicroRNA172 (miR172) functions as a central regulator of flowering time and flower development by post-transcriptional repression of APETALA2-LIKE transcription factors. In the model crop Solanum lycopersicum (tomato), the miR172 family is still poorly annotated and information about the functions of specific members is lacking. Here, de-novo prediction of tomato miR172 coding loci identified seven genes (SlMIR172a-g), that code for four unique species of miR172 (sly-miR172). During reproductive development, sly-miR172s are differentially expressed, with sly-miR172c and sly-miR172d being the most abundant members in developing flowers, and are predicted to guide the cleavage of eight APETALA2-LIKE transcription factors. By CRISPR-Cas9 co-targeting of SlMIR172c and SlMIR172d we have generated a battery of loss-of-function and hypomorphic mutants (slmir172c-d^CR). The slmir172c-d^CR plants developed normal shoot but their flowers displayed graded floral organ abnormalities. Whereas slmir172c^CR loss-of-function caused only a slight greening of petals and stamens, hypomorphic and loss-of-function slmir172d^CR alleles were associated with the conversion of petals and stamens to sepaloids, which were produced in excess. Interestingly, the degrees of floral organ identity alteration and proliferation were directly correlated with the reduction in sly-miR172d activity. These results suggest that sly-miR172d regulates in a dose-dependent manner floral organ identity and number, likely by negatively regulating its APETALA2-like targets.

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

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

In silico identification of conserved miRNAs in the genome of fibre biogenesis crop Corchorus capsularis

Corchorus capsularis, commonly known as jute occupies the leading position in the production of natural fibre alongside lower environmental threat. Small noncoding ~21 to 24 nucleotides long microRNAs play significant roles in regulating the gene expression as well as different functions in cellular growth and development. Here, the study adopted a comprehensive in silico approach to identify and characterize the conserved miRNAs in the genome of C. capsularis including functional annotation of specific gene targets. Expressed Sequence Tags (ESTs) based homology search of 3350 known miRNAs of dicotyledons were allowed against 763 non-redundant ESTs of jute genome, resulted in the prediction of 5 potential miRNA candidates belonging five different miRNA families (miR1536, miR9567-3p, miR4391, miR11300, and miR8689). The putative miRNAs were composed of 18 nucleotides having a range of -0.49 to -1.56 MFEI values and 55%–61% of (A + U) content in their pre-miRNAs. A total of 1052 gene targets of putative miRNAs were identified and their functions were extensively analyzed. Most of the gene targets were involved in plant growth, cell cycle regulation, organelle synthesis, developmental process and environmental responses. Five gene targets, namely, NAC Domain Containing Protein, WRKY DNA binding protein, 3-dehydroquinate synthase, S-adenosyl-L-Met–dependent methyl transferase and Vascular-related NAC-Domain were found to be involved in the lignin biosynthesis, phenylpropanoid pathways and secondary wall formation. The present study might accelerate the more miRNA discovery, strengthening the complete understanding of miRNAs association in the cellular basis of lignin biosynthesis towards the production of high standard jute products.

Redundant and specific roles of individual MIR172 genes in plant development

Evolutionarily conserved microRNAs (miRNAs) usually have high copy numbers in the genome. The redundant and specific roles of each member of a multimember miRNA gene family are poorly understood. Previous studies have shown that the miR156-SPL-miR172 axis constitutes a signaling cascade in regulating plant developmental transitions. Here, we report the feasibility and utility of CRISPR-Cas9 technology to investigate the functions of all 5 MIR172 family members in Arabidopsis. We show that an Arabidopsis plant devoid of miR172 is viable, although it displays pleiotropic morphological defects. MIR172 family members exhibit distinct expression pattern and exert functional specificity in regulating meristem size, trichome initiation, stem elongation, shoot branching, and floral competence. In particular, we find that the miR156-SPL-miR172 cascade is bifurcated into specific flowering responses by matching pairs of coexpressed SPL and MIR172 genes in different tissues. Our results thus highlight the spatiotemporal changes in gene expression that underlie evolutionary novelties of a miRNA gene family in nature. The expansion of MIR172 genes in the Arabidopsis genome provides molecular substrates for the integration of diverse floral inductive cues, which ensures that plants flower at the optimal time to maximize seed yields.

This study uses CRISPR-Cas9 technology to investigate the functions of all five miR172 genes in Arabidopsis, finding that miRNA172 family members exhibit distinct expression pattern and exert functional specificity in regulating meristem size, trichome initiation, stem elongation, shoot branching and floral competence.

A Malvaceae-specific miRNA targeting the newly duplicated GaZIP1L to regulate Zn2+ ion transporter capacity in cotton ovules

MicroRNAs (miRNAs) play critical roles in regulating gene expression in plants, yet their functions underlying cultivated diploid Gossypium arboreum cotton ovule development are largely unknown. Here, we acquired small RNA profiles from G. arboreum ovules and fibers collected at different growth stages, and identified 46 novel miRNAs that accounted for 23.7% of all miRNAs in G. arboreum reported in the latest plant sRNA database. Through analysis of 84 (including 38 conserved) differentially expressed G. arboreum miRNAs, we detected 215 putative protein-coding genes in 26 biological processes as their potential targets. A Malvaceae-specific novel miRNA named gar-miRN44 was found to likely regulate cotton ovule growth by targeting to a newly duplicated Zn²⁺ ion transporter gene GaZIP1L. During cotton ovule development, gar-miRN44 transcript level decreased sharply after 10 to 15 days post-anthesis (DPA), while that of the GaZIP1L increased significantly, with a concomitant increase of Zn²⁺ ion concentration in late ovule developmental stages. Molecular dynamics simulation and ion absorption analysis showed that GaZIP1L has stronger Zn²⁺ ion binding ability than the original GaZIP1, indicating that the newly evolved GaZIP1L may be more suitable for maintaining high Zn²⁺ ion transport capacity that is likely required for cotton ovule growth via enhanced cellulose synthase activities. Our systematic miRNA profiling in G. arboreum and characterization of gar-miRN44 not only contribute to the understanding of miRNA function in cotton, but also provide potential targets for plant breeding.

MiR529a controls plant height, tiller number, panicle architecture and grain size by regulating SPL target genes in rice (Oryza sativa L.)

Rice is one of the most important food crops in the world. Breeding high-yield, multi-resistant and high-quality varieties has always been the goal of rice breeding. Rice tiller, panicle architecture and grain size are the constituent factors of yield, which are regulated by both genetic and environmental factors, including miRNAs, transcription factors, and downstream target genes. Previous studies have shown that SPL (SQUAMOSA PROMOTER BINDING-LIKE) transcription factors can control rice tiller, panicle architecture and grain size, which were regulated by miR156, miR529 and miR535. In this study, we obtained miR529a target mimicry (miR529a-MIMIC) transgenic plants to investigate plant phenotypes, physiological and molecular characteristics together with miR529a overexpression (miR529a-OE) and wild type (WT) to explore the function of miR529a and its SPL target genes in rice. We found that OsSPL2, OsSPL17 and OsSPL18 at seedling stage were regulated by miR529a, but there had complicated mechanism to control plant height. OsSPL2, OsSPL16, OsSPL17 and SPL18 at tillering stage were regulated by miR529a to control plant height and tiller number. And panicle architecture and grain size were controlled by miR529a through altering the expression of all five target genes OsSPL2, OsSPL7, OsSPL14, OsSPL16, OsSPL17 and OsSPL18. Our study suggested that miR529a might control rice growth and development by regulating different SPL target genes at different stages, which could provide a new method to improve rice yield by regulating miR529a and its SPL target genes.

The Crosstalk between MicroRNAs and Gibberellin Signaling in Plants

Gibberellin (GA) is an integral phytohormone that plays prominent roles in controlling seed germination, stem elongation, leaf development and floral induction. It has been shown that GA regulates these diverse biological processes mainly through overcoming the suppressive effects of the DELLA proteins, a family of nuclear repressors of GA response. MicroRNAs (miRNAs), which have been identified as master regulators of gene expression in eukaryotes, are also involved in a wide range of plant developmental events through the repression of their target genes. The pathways of GA biosynthesis and signaling, as well as the pathways of miRNA biogenesis and regulation, have been profoundly delineated in the past several decades. Growing evidence has shown that miRNAs and GAs are coordinated in regulating plant development, as several components in GA pathways are targeted by miRNAs, and GAs also regulate the expression of miRNAs or their target genes vice versa. Here, we review the recent advances in our understanding of the molecular connections between miRNAs and GA, with an emphasis on the two miRNAs, miR156 and miR159.

MiR319a-targeted PtoTCP20 regulates secondary growth via interactions with PtoWOX4 and PtoWND6 in Populus tomentosa

Secondary growth is a key characteristic of trees, which requires the coordination of multiple regulatory mechanisms including transcriptional regulators and microRNAs (miRNAs). However, the roles of microRNAs in the regulation of secondary growth need to be explored in depth. Here, the role of miR319a and its target, PtoTCP20, in the secondary growth of Populus tomentosa stem was investigated using genetic and molecular analyses. The expression level of miR319a gradually decreased from primary to secondary growth in P. tomentosa, while that of PtoTCP20 gradually increased. MiR319a overexpression in seedlings resulted in delayed secondary growth and decreased xylem production, while miR319a knockdown and PtoTCP20 overexpression promoted secondary growth and increased xylem production. Further analysis showed that PtoTCP20 interacted with PtoWOX4a and activated PtoWND6 transcription in vitro and in vivo. Our data show that PtoTCP20 controls vascular cambium proliferation by binding to PtoWOX4a, and promotes secondary xylem differentiation by activating PtoWND6 transcription, thereby regulating secondary growth in P. tomentosa. Our findings provide insight into the molecular mechanisms underlying secondary growth in trees.

WRKY transcription factors WRKY12 and WRKY13 interact with SPL10 to modulate age-mediated flowering

WRKY12 and WRKY13 are two WRKY transcription factors that play important roles in the control of flowering time under short-day (SD) conditions. The temporally regulated expression of WRKY12 and WRKY13 indicates that they may be involved in the age-mediated flowering pathway. However, their roles in this pathway are poorly understood. Here, we show that the transcription of WRKY12 and WRKY13 is directly regulated by SQUAMOSA PROMOTER BINDING-LIKE 10 (SPL10), a transcription factor downstream of the age pathway. Binding and activation analyses revealed that SPL10 functions as a positive regulator of WRKY12 and a negative regulator of WRKY13. Further mechanistic investigation revealed that WRKY12 and WRKY13 physically interact with SPL10 and that both of them bind to the promoter of miR172b. Thus, the WRKY12-SPL10 and WRKY13-SPL10 interactions facilitate and inhibit SPL10 transcriptional function, respectively, to regulate miR172b expression. Together, our results show that WRKY12 and WRKY13 participate in the control of age-mediated flowering under SD conditions though physically interacting with SPLs and co-regulating the target gene miR172b.

miR319a/TCP module and DELLA protein regulate trichome initiation synergistically and improve insect defenses in Populus tomentosa

Trichomes are specialized epidermal cells that contribute to plant resistance against herbivores. Their formation is controlled precisely by multiple genetic and environmental signals. Previous studies have shown that microRNA319 (miR319) and gibberellin (GA) signaling are involved in trichome development in Arabidopsis, but little is known about their interaction between these factors. Here we reported that the miR319a/TEOSINTE BRANCHED/CYCLOIDEA/PCF (TCP) module participates in trichome initiation synergistically with GA signaling in Populus tomentosa. We demonstrated that overexpression of miR319a decreased transcription levels of its targeted TCPs and significantly elevated leaf trichome density in transgenic poplar, resulting in decreasing insect herbivory. Conversely, repressing miR319a by short tandem target mimics (STTM) elevated TCP expression levels and decreased trichome density in transgenic plants. The trichome phenotype of 35S:miR319a plants could be abolished by introducing a miR319a-resistant form of TCP19. Furthermore, the miR319a-targeted TCP19 interacted directly with REPRESSOR OF ga1-3 (RGA), a downstream repressor of GA signaling. TCP19 and RGA synergistically inhibited the GLABROUS1 (GL1)-induced expression of trichome marker gene GLABRA2 (GL2), thereby repressing leaf trichome initiation. Our results provide an insight into the molecular mechanism by which miR319/TCP19 module and GA signaling coordinated regulating trichome initiation in P. tomentosa.

A conserved sequence signature is essential for robust plant miRNA biogenesis

Micro (mi)RNAs are 20–22nt long non-coding RNA molecules involved in post-transcriptional silencing of targets having high base-pair complementarity. Plant miRNAs are processed from long Pol II-transcripts with specific stem-loop structures by Dicer-like (DCL) 1 protein. Although there were reports indicating how a specific region is selected for miRNA biogenesis, molecular details were unclear. Here, we show that the presence of specific GC-rich sequence signature within miRNA/miRNA* region is required for the precise miRNA biogenesis. The involvement of GC-rich signatures in precise processing and abundance of miRNAs was confirmed through detailed molecular and functional analysis. Consistent with the presence of the miRNA-specific GC signature, target RNAs of miRNAs also possess conserved complementary sequence signatures in their miRNA binding motifs. The selection of these GC signatures was dependent on an RNA binding protein partner of DCL1 named HYL1. Finally, we demonstrate a direct application of this discovery for enhancing the abundance and efficiency of artificial miRNAs that are popular in plant functional genomic studies.