PmiREN: a comprehensive encyclopedia of plant miRNAs

LettuceGDB: The community database for lettuce genetics and omics

As a globally popular leafy vegetable and a representative plant of the Asteraceae family, lettuce has great economic and academic significance. In the last decade, high-throughput sequencing, phenotyping, and other multi-omics data in lettuce have accumulated on a large scale, thus increasing the demand for an integrative lettuce database. Here, we report the establishment of a comprehensive lettuce database, LettuceGDB (https://www.lettucegdb.com/). As an omics data hub, the current LettuceGDB includes two reference genomes with detailed annotations; re-sequencing data from over 1000 lettuce varieties; a collection of more than 1300 worldwide germplasms and millions of accompanying phenotypic records obtained with manual and cutting-edge phenomics technologies; re-analyses of 256 RNA sequencing datasets; a complete miRNAome; extensive metabolite information for representative varieties and wild relatives; epigenetic data on the genome-wide chromatin accessibility landscape; and various lettuce research papers published in the last decade. Five hierarchically accessible functions (Genome, Genotype, Germplasm, Phenotype, and O-Omics) have been developed with a user-friendly interface to enable convenient data access. Eight built-in tools (Assembly Converter, Search Gene, BLAST, JBrowse, Primer Design, Gene Annotation, Tissue Expression, Literature, and Data) are available for data downloading and browsing, functional gene exploration, and experimental practice. A community forum is also available for information sharing, and a summary of current research progress on different aspects of lettuce is included. We believe that LettuceGDB can be a comprehensive functional database amenable to data mining and database-driven exploration, useful for both scientific research and lettuce breeding.

Beyond the Primary Structure of Nucleic Acids: Potential Roles of Epigenetics and Noncanonical Structures in the Regulations of Plant Growth and Stress Responses

Epigenetics deals with changes in gene expression that are not caused by modifications in the primary sequence of nucleic acids. These changes beyond primary structures of nucleic acids not only include DNA/RNA methylation, but also other reversible conversions, together with histone modifications or RNA interference. In addition, under particular conditions (such as specific ion concentrations or protein-induced stabilization), the right-handed double-stranded DNA helix (B-DNA) can form noncanonical structures commonly described as "non-B DNA" structures. These structures comprise, for example, cruciforms, i-motifs, triplexes, and G-quadruplexes. Their formation often leads to significant differences in replication and transcription rates. Noncanonical RNA structures have also been documented to play important roles in translation regulation and the biology of noncoding RNAs. In human and animal studies, the frequency and dynamics of noncanonical DNA and RNA structures are intensively investigated, especially in the field of cancer research and neurodegenerative diseases. In contrast, noncanonical DNA and RNA structures in plants have been on the fringes of interest for a long time and only a few studies deal with their formation, regulation, and physiological importance for plant stress responses. Herein, we present a review focused on the main fields of epigenetics in plants and their possible roles in stress responses and signaling, with special attention dedicated to noncanonical DNA and RNA structures.

sRNAtoolbox: Dockerized Analysis of Small RNA Sequencing Data in Model and Non-model Species

The current versions of the microRNA databases MiRgeneDB, miRBase, and PmiREN contain annotations for a total of 358 different species. Public repositories, however, host small RNA sequencing data for over 800 species. This discrepancy implies that microRNA research is also very active in species that neither have an available high-quality genome assembly nor annotations for microRNAs or other types of noncoding genes. These cases are particularly challenging to analyze because reference sequences need to be collected from different sources and processed and formatted appropriately so that the dedicated small RNA analysis tools can make use of them. In this protocol we describe how small RNA sequencing data can be easily analyzed by means of a dockerized version of the well-established sRNAtoolbox/sRNAbench small RNA tools. We outline the analysis of two publicly available datasets to demonstrate basic aspects like the preparation of the local database, expression profiling, or differential expression analysis as well as more advanced features such as quantification of exogenous RNA content and data analysis in non-model species.

Plant MicroRNA Identification and Annotation Using Deep Sequencing Data

MicroRNAs (miRNAs) are endogenous non-coding small RNAs, which regulate gene expression at the post-transcriptional level. A large number of studies have revealed that they play key roles in diverse life activities, such as growth and development. In the last decade, deep sequencing technology has generated substantial small RNA sequencing (sRNA-Seq) data. Meanwhile, numerous tools have been developed to identify miRNAs from these sRNA-Seq data, resulting in a surge of miRNA annotations. Among these tools, the series of miRDeep-P and miRDeep-P2 have been widely used in plant miRNA annotation. Here, we employed miRDeep-P2 to demonstrate the plant miRNA annotation processes step by step using the deep sequencing data.

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.

Surface-localized glycoproteins act through class C ARFs to fine-tune gametophore initiation in Physcomitrium patens

Arabinogalactan proteins are functionally diverse cell wall structural glycoproteins that have been implicated in cell wall remodeling, although the mechanistic actions remain elusive. Here, we identify and characterize two AGP glycoproteins, SLEEPING BEAUTY (SB) and SB-like (SBL), that negatively regulate the gametophore bud initiation in Physcomitrium patens by dampening cell wall loosening/softening. Disruption of SB and SBL led to accelerated gametophore formation and altered cell wall compositions. The function of SB is glycosylation dependent and genetically connected with the class C auxin response factor (ARF) transcription factors PpARFC1B and PpARFC2. Transcriptomics profiling showed that SB upregulates PpARFC2, which in turn suppresses a range of cell wall-modifying genes that are required for cell wall loosening/softening. We further show that PpARFC2 binds directly to multiple AuxRE motifs on the cis-regulatory sequences of PECTIN METHYLESTERASE to suppress its expression. Hence, our results demonstrate a mechanism by which the SB modulates the strength of intracellular auxin signaling output, which is necessary to fine-tune the timing of gametophore initials formation.

Transcriptional Regulation of zma-MIR528a by Action of Nitrate and Auxin in Maize

In recent years, miR528, a monocot-specific miRNA, has been assigned multifaceted roles during development and stress response in several plant species. However, the transcription regulation and the molecular mechanisms controlling MIR528 expression in maize are still poorly explored. Here we analyzed the zma-MIR528a promoter region and found conserved transcription factor binding sites related to diverse signaling pathways, including the nitrate (TGA1/4) and auxin (AuxRE) response networks. Accumulation of both pre-miR528a and mature miR528 was up-regulated by exogenous nitrate and auxin treatments during imbibition, germination, and maize seedling establishment. Functional promoter analyses demonstrated that TGA1/4 and AuxRE sites are required for transcriptional induction by both stimuli. Overall, our findings of the nitrogen- and auxin-induced zma-MIR528a expression through cis-regulatory elements in its promoter contribute to the knowledge of miR528 regulome.

Evolutionary characterization of miR396s in Poaceae exemplified by their genetic effects in wheat and maize

MiR396s play important roles in regulating plant growth and stress response, and great potential for crop yield promotion was anticipated. For more comprehensive and precise understanding of miR396s in Poaceae, we analyzed the phylogenetic linkage, gene expression, and chromosomal distribution of miR396s in this study. Although the mature miR396s' sequences were mostly conserved, differential expression patterns and chromosomal distribution were found among Poaceae species including the major cereal crops rice, wheat, and maize. Consistently, in comparison with rice, wheat and maize plants transformed with the target mimicry construct of miR396 (MIM396) exhibited differential effects on grain size and disease resistance. While the TaMIM396 plants showed increased grain size, panicle length and sensitivity to B. graminis, the ZmMIM396 plants didn't show obvious changes in grain size and disease resistance. In Addition, several GROWTH-REGULATING FACTOR (GRF) genes in wheat and maize were repressed by miR396s, which could be reversed by MIM396, confirming the conserved regulatory roles of miR396 on GRFs. While providing new solution to enhance grain yield in wheat and revealing potential regulatory variations of miR396s in controlling grain size and disease resistance in different crops, this study gives clues to further explore miR396s' functions in other Poaceae species.

Identification of Species-Specific MicroRNAs Provides Insights into Dynamic Evolution of MicroRNAs in Plants

MicroRNAs (miRNAs) are an important class of regulatory small RNAs that program gene expression, mainly at the post-transcriptional level. Although sporadic examples of species-specific miRNAs (termed SS-miRNAs) have been reported, a genome-scale study across a variety of distant species has not been assessed. Here, by comprehensively analyzing miRNAs in 81 plant species phylogenetically ranging from chlorophytes to angiosperms, we identified 8048 species-specific miRNAs from 5499 families, representing over 61.2% of the miRNA families in the examined species. An analysis of the conservation from different taxonomic levels supported the high turnover rate of SS-miRNAs, even over short evolutionary distances. A comparison of the intrinsic features between SS-miRNAs and NSS-miRNAs (non-species-specific miRNAs) indicated that the AU content of mature miRNAs was the most striking difference. Our data further illustrated a significant bias of the genomic coordinates towards SS-miRNAs lying close to or within genes. By analyzing the 125,267 putative target genes for the 7966 miRNAs, we found the preferentially regulated functions of SS-miRNAs related to diverse metabolic processes. Collectively, these findings underscore the dynamic evolution of miRNAs in the species-specific lineages.

Degradome Assisted Plant MicroRNA Prediction Under Alternative Annotation Criteria

Current microRNA (miRNA) prediction methods are generally based on annotation criteria that tend to miss potential functional miRNAs. Recently, new miRNA annotation criteria have been proposed that could lead to improvements in miRNA prediction methods in plants. Here, we investigate the effect of the new criteria on miRNA prediction in Arabidopsis thaliana and present a new degradome assisted functional miRNA prediction approach. We investigated the effect by applying the new criteria, and a more permissive criteria on miRNA prediction using existing miRNA prediction tools. We also developed an approach to miRNA prediction that is assisted by the functional information extracted from the analysis of degradome sequencing. We demonstrate the improved performance of degradome assisted miRNA prediction compared to unassisted prediction and evaluate the approach using miRNA differential expression analysis. We observe how the miRNA predictions fit under the different criteria and show a potential novel miRNA that has been missed within Arabidopsis thaliana. Additionally, we introduce a freely available software 'PAREfirst' that employs the degradome assisted approach. The study shows that some miRNAs could be missed due to the stringency of the former annotation criteria, and combining a degradome assisted approach with more permissive miRNA criteria can expand confident miRNA predictions.

RiceNCexp: a rice non-coding RNA co-expression atlas based on massive RNA-seq and small-RNA seq data

Non-coding RNAs (ncRNAs) play important roles in regulating expression of protein-coding genes. Although gene expression databases have emerged in a timely manner, a comprehensive expression database for ncRNAs is still lacking. Herein, we constructed a rice ncRNA co-expression atlas (RiceNCexp), based on 491 RNA-seq and 274 small RNA (sRNA)-seq datasets. RiceNCexp hosts four types of ncRNAs, namely lncRNAs, PHAS genes, miRNAs, and phasiRNAs. RiceNCexp provides comprehensive expression information for rice ncRNAs in 22 tissues/organs, an efficient tau-based mining tool for tissue-specific ncRNAs, and the robust co-expression analysis among ncRNAs or between ncRNAs and protein-coding genes, based on 116 pairs of RNA-seq and sRNA-seq libraries from the same experiments. In summary, RiceNCexp is a user-friendly and comprehensive rice ncRNA co-expression atlas and can be freely accessed at https://cbi.njau.edu.cn/RiceNCexp/.

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.

Regulatory non-coding RNA: The core defense mechanism against plant pathogens

Plant pathogens damage crops and threaten global food security. Plants have evolved complex defense networks against pathogens, using crosstalk among various signaling pathways. Key regulators conferring plant immunity through signaling pathways include protein-coding genes and non-coding RNAs (ncRNAs). The discovery of ncRNAs in plant transcriptomes was first considered "transcriptional noise". Recent reviews have highlighted the importance of non-coding RNAs. However, understanding interactions among different types of noncoding RNAs requires additional research. This review attempts to consider how long-ncRNAs, small-ncRNAs and circular RNAs interact in response to pathogenic diseases within different plant species. Developments within genomics and bioinformatics could lead to the further discovery of plant ncRNAs, knowledge of their biological roles, as well as an understanding of their importance in exploiting the recent molecular-based technologies for crop protection.

Genome-Wide Identification and Analysis of Lbd Transcription Factor Genes in Jatropha curcas and Related Species

Lateral organ boundaries domain (LBD) proteins are plant-specific transcription factors that play important roles in organ development and stress response. However, the function of LBD genes has not been reported in Euphorbiaceae. In this paper, we used Jatropha curcas as the main study object and added rubber tree (Hevea brasiliensis), cassava (Manihot esculenta Crantz) and castor (Ricinus communis L.) to take a phylogenetic analysis of LBD genes. Of LBD, 33, 58, 54 and 30 members were identified in J. curcas, rubber tree, cassava and castor, respectively. The phylogenetic analysis showed that LBD members of Euphorbiaceae could be classified into two major classes and seven subclasses (Ia-Ie,IIa-IIb), and LBD genes of Euphorbiaceae tended to cluster in the same branch. Further analysis showed that the LBD genes of Euphorbiaceae in the same clade usually had similar protein motifs and gene structures, and tissue expression patterns showed that they also have similar expression profiles. JcLBDs in class Ia and Ie are mainly expressed in male and female flowers, and there are multiple duplication genes with similar expression profiles in these clades. It was speculated that they are likely to play important regulatory roles in flower development. Our study provided a solid foundation for further investigation of the role of LBD genes in the sexual differentiaion of J. curcas.

Molecular evolution and functional modification of plant miRNAs with CRISPR

Gene editing using clustered regularly interspaced short palindromic repeat/CRISPR-associated proteins (CRISPR/Cas) has revolutionized biotechnology and provides genetic tools for medicine and life sciences. However, the application of this technology to miRNAs, with the function as negative gene regulators, has not been extensively reviewed in plants. Here, we summarize the evolution, biogenesis, and structure of miRNAs, as well as their interactions with mRNAs and computational models for predicting target genes. In addition, we review current advances in CRISPR/Cas for functional analysis and for modulating miRNA genes in plants. Extending our knowledge of miRNAs and their manipulation with CRISPR will provide fundamental understanding of the functions of plant miRNAs and facilitate more sustainable and publicly acceptable genetic engineering of crops.

Characteristics of members of IGT family genes in controlling rice root system architecture and tiller development

Root system architecture (RSA) and tiller are important agronomic traits. However, the mechanisms of the IGT family genes regulate RSA and tiller development in different rice varieties remain unclear. In this study, we demonstrated that 38 rice varieties obtained from Yuanyang Hani's terraced fields with different RSA and could be classified into six groups based on the ratio of root length and width. We found a positive correlation between RSA (including root width, length, and area) and tiller number in most of rice varieties. Furthermore, the IGT family genes Deeper Rooting 1 (DRO1), LAZY1, TAC1, and qSOR1 showed different expression patterns when rice grown under irrigation and drought conditions. Moreover, the qSOR1 gene had higher levels in the roots and tillers, and accompanied with higher levels of PIN1b gene in roots when rice grown under drought environmental condition. DRO1 gene had two single nucleotide polymorphisms (SNPs) in the exon 3 sequences and showed different expression patterns in the roots and tillers of the 38 rice varieties. Overexpression of DRO1 with a deletion of exon 5 caused shorter root length, less lateral roots and lower levels of LAZY1, TAC1, and qSOR1. Further protein interaction network, microRNA targeting and co-expression analysis showed that DRO1 plays a critical role in the root and tiller development associated with auxin transport. These data suggest that the RSA and tiller development are regulated by the IGT family genes in an intricate network way, which is tightly related to rice genetic background in rice adapting to different environmental conditions.

The miR166–mRNA network regulates vascular tissue differentiation in Moso bamboo

miR166s play an important role in plant tissue differentiation. However, the functions of miR166s in the differentiation of vascular tissue in bamboo have not yet been elucidated. Here, we showed that five miR166s are overexpressed (tags per million reads > 2,000) in underground shoot samples of wild-type (WT) Moso bamboo (Phyllostachys edulis) and a thick-walled variant (P. edulis “Pachyloen”) throughout the developmental process. Potential targets of these miR166s include some genes encoding homeodomain-leucine zipper (HD-ZIP) transcription factors and protein kinases. Cleavage sites for miR166s were identified in seven PeHD-ZIP homologs and a protein kinase gene via degradome sequencing (p < 0.05). Dual-luciferase and transient expression assays confirmed the binding of miR166s to PeHOXs. Fluorescence in situ hybridization revealed that miR166s were localized to the xylem of the leaf, root, and internode of 2-month-old pot seedlings of WT Moso bamboo. Overall, these findings reveal that miR166s are regulators of vascular tissue differentiation in bamboo. The miR166s identified in our study provide novel targets for bamboo breeding.

Biogenesis, conservation, and function of miRNA in liverworts

MicroRNAs (miRNAs) are small non-coding endogenous RNA molecules, 18-24 nucleotides long, that control multiple gene regulatory pathways via post-transcriptional gene silencing in eukaryotes. To develop a comprehensive picture of the evolutionary history of miRNA biogenesis and action in land plants, studies on bryophyte representatives are needed. Here, we review current understanding of liverwort MIR gene structure, miRNA biogenesis, and function, focusing on the simple thalloid Pellia endiviifolia and the complex thalloid Marchantia polymorpha. We review what is known about conserved and non-conserved miRNAs, their targets, and the functional implications of miRNA action in M. polymorpha and P. endiviifolia. We note that most M. polymorpha miRNAs are encoded within protein-coding genes and provide data for 23 MIR gene structures recognized as independent transcriptional units. We identify M. polymorpha genes involved in miRNA biogenesis that are homologous to those identified in higher plants, including those encoding core microprocessor components and other auxiliary and regulatory proteins that influence the stability, folding, and processing of pri-miRNAs. We analyzed miRNA biogenesis proteins and found similar domain architecture in most cases. Our data support the hypothesis that almost all miRNA biogenesis factors in higher plants are also present in liverworts, suggesting that they emerged early during land plant evolution.

Interactive Web-based Annotation of Plant MicroRNAs with iwa-miRNA

MicroRNAs (miRNAs) are important regulators of gene expression. The large-scale detection and profiling of miRNAs have been accelerated with the development of high-throughput small RNA sequencing (sRNA-Seq) techniques and bioinformatics tools. However, generating high-quality comprehensive miRNA annotations remains challenging due to the intrinsic complexity of sRNA-Seq data and inherent limitations of existing miRNA prediction tools. Here, we present iwa-miRNA, a Galaxy-based framework that can facilitate miRNA annotation in plant species by combining computational analysis and manual curation. iwa-miRNA is specifically designed to generate a comprehensive list of miRNA candidates, bridging the gap between already annotated miRNAs provided by public miRNA databases and new predictions from sRNA-Seq datasets. It can also assist users in selecting promising miRNA candidates in an interactive mode, contributing to the accessibility and reproducibility of genome-wide miRNA annotation. iwa-miRNA is user-friendly and can be easily deployed as a web application for researchers without programming experience. With flexible, interactive, and easy-to-use features, iwa-miRNA is a valuable tool for the annotation of miRNAs in plant species with reference genomes. We also illustrate the application of iwa-miRNA for miRNA annotation using data from plant species with varying genomic complexity. The source codes and web server of iwa-miRNA are freely accessible at http://iwa-miRNA.omicstudio.cloud/.

An Insight Into Pentatricopeptide-Mediated Chloroplast Necrosis via microRNA395a During Rhizoctonia solani Infection

Sheath blight (ShB) disease, caused by Rhizoctonia solani, is one of the major biotic stress-oriented diseases that adversely affect the rice productivity worldwide. However, the regulatory mechanisms are not understood yet comprehensively. In the current study, we had investigated the potential roles of miRNAs in economically important indica rice variety Pusa Basmati-1 upon R. solani infection by carrying out in-depth, high-throughput small RNA sequencing with a total data size of 435 million paired-end raw reads from rice leaf RNA samples collected at different time points. Detailed data analysis revealed a total of 468 known mature miRNAs and 747 putative novel miRNAs across all the libraries. Target prediction and Gene Ontology functional analysis of these miRNAs were found to be unraveling various cellular, molecular, and biological functions by targeting various plant defense-related genes. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was performed to validate the miRNAs and their putative target genes. Out of the selected miRNA-specific putative target genes, miR395a binding and its cleavage site on pentatricopeptide were determined by 5’ RACE-PCR. It might be possible that R. solani instigated chloroplast degradation by modulating the pentatricopeptide which led to increased susceptibility to fungal infection.

Construction and characterization of a de novo draft genome of garden cress (Lepidium sativum L.)

Garden cress (Lepidium sativum L.) is a Brassicaceae crop recognized as a healthy vegetable and a medicinal plant. Lepidium is one of the largest genera in Brassicaceae, yet, the genus has not been a focus of extensive genomic research. In the present work, garden cress genome was sequenced using the long read high-fidelity sequencing technology. A de novo, draft genome assembly that spans 336.5 Mb was produced, corresponding to 88.6% of the estimated genome size and representing 90% of the evolutionarily expected orthologous gene content. Protein coding gene content was structurally predicted and functionally annotated, resulting in the identification of 25,668 putative genes. A total of 599 candidate disease resistance genes were identified by predicting resistance gene domains in gene structures, and 37 genes were detected as orthologs of heavy metal associated protein coding genes. In addition, 4289 genes were assigned as "transcription factor coding." Six different machine learning algorithms were trained and tested for their performance in classifying miRNA coding genomic sequences. Logistic regression proved the best performing trained algorithm, thus utilized for pre-miRNA coding loci identification in the assembly. Repetitive DNA analysis involved the characterization of transposable element and microsatellite contents. L. sativum chloroplast genome was also assembled and functionally annotated. Data produced in the present work is expected to constitute a foundation for genomic research in garden cress and contribute to genomics-assisted crop improvement and genome evolution studies in the Brassicaceae family.

Regulatory role of microRNAs (miRNAs) in the recent development of abiotic stress tolerance of plants

MicroRNAs (miRNAs) are a distinct groups of single-stranded non-coding, tiny regulatory RNAs approximately 20-24 nucleotides in length. miRNAs negatively influence gene expression at the post-transcriptional level and have evolved considerably in the development of abiotic stress tolerance in a number of model plants and economically important crop species. The present review aims to deliver the information on miRNA-mediated regulation of the expression of major genes or Transcription Factors (TFs), as well as genetic and regulatory pathways. Also, the information on adaptive mechanisms involved in plant abiotic stress responses, prediction, and validation of targets, computational tools, and databases available for plant miRNAs, specifically focus on their exploration for engineering abiotic stress tolerance in plants. The regulatory function of miRNAs in plant growth, development, and abiotic stresses consider in this review, which uses high-throughput sequencing (HTS) technologies to generate large-scale libraries of small RNAs (sRNAs) for conventional screening of known and novel abiotic stress-responsive miRNAs adds complexity to regulatory networks in plants. The discoveries of miRNA-mediated tolerance to multiple abiotic stresses, including salinity, drought, cold, heat stress, nutritional deficiency, UV-radiation, oxidative stress, hypoxia, and heavy metal toxicity, are highlighted and discussed in this review.

Genome-wide identification, classification and expression analysis of the JmjC domain-containing histone demethylase gene family in Jatropha curcas L

JmjC domain-containing proteins, an important family of histone lysine demethylase, play significant roles in maintaining the homeostasis of histone methylation. In this study, we comprehensively analyzed the JmjC domain-containing gene family in Jatropha curcas and found 20 JmjC domain-containing genes (JcJMJ genes). Phylogenetic analysis revealed that these JcJMJ genes can be classified into five major subgroups, and genes in each subgroup had similar motif and domain composition. Cis-regulatory element analysis showed that the number and types of cis-regulatory elements owned by the promoter of JcJMJ genes in different subgroup were significantly different. Moreover, miRNA target prediction result revealed a complicated miRNA-mediated post-transcriptional regulatory network, in which JcJMJ genes were regulated by different numbers and types of miRNAs. Further analysis of the tissue and stress expression profiles showed that many JcJMJ genes had tissue and stress expression specificity. All these results provided valuable information for understanding the evolution of JcJMJ genes and the complex transcriptional and post transcriptional regulation involved, and laid the foundation for further functional analysis of JcJMJ genes.

Recent advances and potential applications of cross-kingdom movement of miRNAs in modulating plant's disease response

In the recent past, cross-kingdom movement of miRNAs, small (20–25 bases), and endogenous regulatory RNA molecules has emerged as one of the major research areas to understand the potential implications in modulating the plant’s biotic stress response. The current review discussed the recent developments in the mechanism of cross-kingdom movement (long and short distance) and critical cross-talk between host’s miRNAs in regulating gene function in bacteria, fungi, viruses, insects, and nematodes, and vice-versa during host-pathogen interaction and their potential implications in crop protection. Moreover, cross-kingdom movement during symbiotic interaction, the emerging role of plant’s miRNAs in modulating animal’s gene function, and feasibility of spray-induced gene silencing (SIGS) in combating biotic stresses in plants are also critically evaluated. The current review article analysed the horizontal transfer of miRNAs among plants, animals, and microbes that regulates gene expression in the host or pathogenic organisms, contributing to crop protection. Further, it highlighted the challenges and opportunities to harness the full potential of this emerging approach to mitigate biotic stress efficiently.

Walnut N-Acetylserotonin Methyltransferase Gene Family Genome-Wide Identification and Diverse Functions Characterization During Flower Bud Development

Melatonin widely mediates multiple developmental dynamics in plants as a vital growth stimulator, stress protector, and developmental regulator. N-acetylserotonin methyltransferase (ASMT) is the key enzyme that catalyzes the final step of melatonin biosynthesis in plants and plays an essential role in the plant melatonin regulatory network. Studies of ASMT have contributed to understanding the mechanism of melatonin biosynthesis in plants. However, AMST gene is currently uncharacterized in most plants. In this study, we characterized the JrASMT gene family using bioinformatics in a melatonin-rich plant, walnut. Phylogenetic, gene structure, conserved motifs, promoter elements, interacting proteins and miRNA analyses were also performed. The expansion and differentiation of the ASMT family occurred before the onset of the plant terrestrialization. ASMT genes were more differentiated in dicotyledonous plants. Forty-six ASMT genes were distributed in clusters on 10 chromosomes of walnut. Four JrASMT genes had homologous relationships both within walnut and between species. Cis-regulatory elements showed that JrASMT was mainly induced by light and hormones, and targeted cleavage of miRNA172 and miR399 may be an important pathway to suppress JrASMT expression. Transcriptome data showed that 13 JrASMT were differentially expressed at different periods of walnut bud development. WGCNA showed that JrASMT1/10/13/23 were coexpressed with genes regulating cell fate and epigenetic modifications during early physiological differentiation of walnut female flower buds. JrASMT12/28/37/40 were highly expressed during morphological differentiation of flower buds, associated with altered stress capacity of walnut flower buds, and predicted to be involved in the regulatory network of abscisic acid, salicylic acid, and cytokinin in walnut. The qRT-PCR validated the results of differential expression analysis and further provided three JrASMT genes with different expression profiles in walnut flower bud development. Our study explored the evolutionary relationships of the plant ASMT gene family and the functional characteristics of walnut JrASMT. It provides a valuable perspective for further understanding the complex melatonin mechanisms in plant developmental regulation.

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

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

Transcriptome and miRNAome analysis reveals components regulating tissue differentiation of bamboo shoots

Primary thickening determines bamboo yield and wood property. However, little is known about the regulatory networks involved in this process. This study identified a total of 58,652 genes and 150 miRNAs via transcriptome and small RNA sequencing using the underground thickening shoot samples of wild-type (WT) Moso bamboo (Phyllostachys edulis) and a thick wall (TW) variant (P. edulis "Pachyloen") at five developmental stages (WTS1/TWS1-WTS5/TWS5). A total of 14,029 (65.17%) differentially expressed genes and 68 (45.33%) differentially expressed miRNAs were identified from the WT, TW, and WTTW groups. The first two groups were composed of four pairwise combinations, each between two successive stages (WTS2/TWS2_versus_WTS1/TWS1, WTS3/TWS3_versus_WTS2/TWS2, WTS4/TWS4_versus_WTS3/TWS3, and WTS5/TWS5_versus_WTS4/TWS4), and the WTTW group was composed of five combinations, each between two relative stages (TWS1-5_versus_WTS1-5). Additionally, among the phytohormones, zeatin showed more remarkable changes in concentrations than indole-3-acetic acid, gibberellic acid, and abscisic acid throughout the five stages in the WT and the TW groups. Moreover, 125 cleavage sites were identified for 387 miRNA-mRNA pairs via degradome sequencing (P < 0.05). The dual-luciferase reporter assay confirmed that 13 miRNAs bound to 12 targets. Fluorescence in situ hybridization localized miR166 and miR160 in the shoot apical meristem and the procambium of Moso bamboo shoots at the S1 stage. Thus, primary thickening is a complex process regulated by miRNA-gene-phytohormone networks, and the miRNAome and transcriptome dynamics regulate phenotypic plasticity. These findings provide insights into the molecular mechanisms underlying wood formation and properties and propose targets for bamboo breeding.

Genome wide study of cysteine rich receptor like proteins in Gossypium sp

Cysteine-rich receptor-like-kinases (CRKs), a transmembrane subfamily of receptor-like kinase, play crucial roles in plant adaptation. As such cotton is the major source of fiber for the textile industry, but environmental stresses are limiting its growth and production. Here, we have performed a deep computational analysis of CRKs in five Gossypium species, including G. arboreum (60 genes), G. raimondii (74 genes), G. herbaceum (65 genes), G. hirsutum (118 genes), and G. barbadense (120 genes). All identified CRKs were classified into 11 major classes and 43 subclasses with the finding of several novel CRK-associated domains including ALMT, FUSC_2, Cript, FYVE, and Pkinase. Of these, DUF26_DUF26_Pkinase_Tyr was common and had elevated expression under different biotic and abiotic stresses. Moreover, the 35 land plants comparison identified several new CRKs domain-architectures. Likewise, several SNPs and InDels were observed in CLCuD resistant G. hirsutum. The miRNA target side prediction and their expression profiling in different tissues predicted miR172 as a major CRK regulating miR. The expression profiling of CRKs identified multiple clusters with co-expression under certain stress conditions. The expression analysis under CLCuD highlighted the role of GhCRK057, GhCRK059, GhCRK058, and GhCRK081 in resistant accession. Overall, these results provided primary data for future potential functional analysis as well as a reference study for other agronomically important crops.

Genome-wide identification of MITE-derived microRNAs and their targets in bread wheat

Background

Plant miRNAs are a class of small non-coding RNAs that can repress gene expression at the post-transcriptional level by targeting RNA degradation or promoting translational repression. There is increasing evidence that some miRNAs can derive from a group of non-autonomous class II transposable elements called Miniature Inverted-repeat Transposable Elements (MITEs).

Results

We used public small RNA and degradome libraries from Triticum aestivum to screen for microRNAs production and predict their cleavage target sites. In parallel, we also created a comprehensive wheat MITE database by identifying novel elements and compiling known ones. When comparing both data sets, we found high homology between MITEs and 14% of all the miRNAs production sites detected. Furthermore, we show that MITE-derived miRNAs have preference for targeting degradation sites with MITE insertions in the 3’ UTR regions of the transcripts.

Conclusions

Our results revealed that MITE-derived miRNAs can underlay the origin of some miRNAs and potentially shape a regulatory gene network. Since MITEs are found in millions of insertions in the wheat genome and are closely linked to genic regions, this kind of regulatory network could have a significant impact on the post-transcriptional control of gene expression.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08364-4).

Comprehensive Annotation and Functional Exploration of MicroRNAs in Lettuce

MicroRNA (miRNA) is an important endogenous post-transcriptional regulator, while lettuce (Lactuca sativa) is a leafy vegetable of global economic significance. However, there are few studies on miRNAs in lettuce, and research on miRNA regulatory network in lettuce is absent. In this study, through deep sequencing of small RNAs in different tissues, together with a reference genome, 157 high-confidence miRNA loci in lettuce were comprehensively identified, and their expression patterns were determined. Using a combination of computational prediction and high-throughput experimental verification, a set of reliable lettuce miRNA targets were obtained. Furthermore, through RNA-Seq, the expression profiles of these targets and a comprehensive view of the negative regulatory relationship between miRNAs and their targets was acquired based on a correlation analysis. To further understand miRNA functions, a miRNA regulatory network was constructed, with miRNAs at the core and combining transcription factors and miRNA target genes. This regulatory network, mainly composed of feed forward loop motifs, greatly increases understanding of the potential functions of miRNAs, and many unknown potential regulatory links were discovered. Finally, considering its specific expression pattern, Lsa-MIR408 as a hub gene was employed to illustrate the function of the regulatory network, and genetic experiments revealed its ability to increase the fresh weight and achene size of lettuce. In short, this work lays a solid foundation for the study of miRNA functions and regulatory networks in lettuce.

MicroRNA-mediated bioengineering for climate-resilience in crops

Global projections on the climate change and the dynamic environmental perturbations indicate severe impacts on food security in general, and crop yield, vigor and the quality of produce in particular. Sessile plants respond to environmental challenges such as salt, drought, temperature, heavy metals at transcriptional and/or post-transcriptional levels through the stress-regulated network of pathways including transcription factors, proteins and the small non-coding endogenous RNAs. Amongs these, the miRNAs have gained unprecedented attention in recent years as key regulators for modulating gene expression in plants under stress. Hence, tailoring of miRNAs and their target pathways presents a promising strategy for developing multiple stress-tolerant crops. Plant stress tolerance has been successfully achieved through the over expression of microRNAs such as Os-miR408, Hv-miR82 for drought tolerance; OsmiR535A and artificial DST miRNA for salinity tolerance; and OsmiR535 and miR156 for combined drought and salt stress. Examples of miR408 overexpression also showed improved efficiency of irradiation utilization and carbon dioxide fixation in crop plants. Through this review, we present the current understanding about plant miRNAs, their roles in plant growth and stress-responses, the modern toolbox for identification, characterization and validation of miRNAs and their target genes including in silico tools, machine learning and artificial intelligence. Various approaches for up-regulation or knock-out of miRNAs have been discussed. The main emphasis has been given to the exploration of miRNAs for development of bioengineered climate-smart crops that can withstand changing climates and stressful environments, including combination of stresses, with very less or no yield penalties.

MirGeneDB 2.1: toward a complete sampling of all major animal phyla

We describe an update of MirGeneDB, the manually curated microRNA gene database. Adhering to uniform and consistent criteria for microRNA annotation and nomenclature, we substantially expanded MirGeneDB with 30 additional species representing previously missing metazoan phyla such as sponges, jellyfish, rotifers and flatworms. MirGeneDB 2.1 now consists of 75 species spanning over ∼800 million years of animal evolution, and contains a total number of 16 670 microRNAs from 1549 families. Over 6000 microRNAs were added in this update using ∼550 datasets with ∼7.5 billion sequencing reads. By adding new phylogenetically important species, especially those relevant for the study of whole genome duplication events, and through updating evolutionary nodes of origin for many families and genes, we were able to substantially refine our nomenclature system. All changes are traceable in the specifically developed MirGeneDB version tracker. The performance of read-pages is improved and microRNA expression matrices for all tissues and species are now also downloadable. Altogether, this update represents a significant step toward a complete sampling of all major metazoan phyla, and a widely needed foundation for comparative microRNA genomics and transcriptomics studies. MirGeneDB 2.1 is part of RNAcentral and Elixir Norway, publicly and freely available at http://www.mirgenedb.org/.

Identification of GROWTH-REGULATING FACTOR transcription factors in lettuce (Lactuca sativa) genome and functional analysis of LsaGRF5 in leaf size regulation

BACKGROUND

GROWTH-REGULATING FACTORs (GRFs), a type of plant-specific transcription factors, play important roles in regulating plant growth and development. Although GRF gene family has been identified in various plant species, a genome-wide analysis of this family in lettuce (Lactuca sativa L.) has not been reported yet.

RESULTS

Here we identified 15 GRF genes in lettuce and performed comprehensive analysis of them, including chromosomal locations, gene structures, and conserved motifs. Through phylogenic analysis, we divided LsaGRFs into six groups. Transactivation assays and subcellular localization of LsaGRF5 showed that this protein is likely to act as a transcriptional factor in the cell nucleus. Furthermore, transgenic lettuce lines overexpressing LsaGRF5 exhibited larger leaves, while smaller leaves were observed in LsaMIR396a overexpression lines, in which LsaGRF5 was down-regulated.

CONCLUSIONS

These results in lettuce provide insight into the molecular mechanism of GRF gene family in regulating leaf growth and development and foundational information for genetic improvement of the lettuce variations specialized in leaf character.

Mining plant endogenous target mimics from miRNA-lncRNA interactions based on dual-path parallel ensemble pruning method

The interactions between microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) play important roles in biological activities. Specially, lncRNAs as endogenous target mimics (eTMs) can bind miRNAs to regulate the expressions of target messenger RNAs (mRNAs). A growing number of studies focus on animals, but the studies on plants are scarce and many functions of plant eTMs are unknown. This study proposes a novel ensemble pruning protocol for predicting plant miRNA-lncRNA interactions at first. It adaptively prunes the base models based on dual-path parallel ensemble method to meet the challenge of cross-species prediction. Then potential eTMs are mined from predicted results. The expression levels of RNAs are identified through biological experiment to construct the lncRNA-miRNA-mRNA regulatory network, and the functions of potential eTMs are inferred through enrichment analysis. Experiment results show that the proposed protocol outperforms existing methods and state-of-the-art predictors on various plant species. A total of 17 potential eTMs are verified by biological experiment to involve in 22 regulations, and 14 potential eTMs are inferred by Gene Ontology enrichment analysis to involve in 63 functions, which is significant for further research.

PmiREN2.0: from data annotation to functional exploration of plant microRNAs

Nearly 200 plant genomes have been sequenced over the last two years, and new functions of plant microRNAs (miRNAs) have been revealed. Therefore, timely update of the plant miRNA databases by incorporating miRNAs from the newly sequenced species and functional information is required to provide useful resources for advancing plant miRNA research. Here we report the update of PmiREN2.0 (https://pmiren.com/) with an addition of 19 363 miRNA entries from 91 plants, doubling the amount of data in the original version. Meanwhile, abundant regulatory information centred on miRNAs was added, including predicted upstream transcription factors through binding motifs scanning and elaborate annotation of miRNA targets. As an example, a genome-wide regulatory network centred on miRNAs was constructed for Arabidopsis. Furthermore, phylogenetic trees of conserved miRNA families were built to expand the understanding of miRNA evolution across the plant lineages. These data are helpful to deduce the regulatory relationships concerning miRNA functions in diverse plants. Beside the new data, a suite of design tools was incorporated to facilitate experimental practice. Finally, a forum named ‘PmiREN Community’ was added for discussion and resource and new discovery sharing. With these upgrades, PmiREN2.0 should serve the community better and accelerate miRNA research in plants.

Genome-wide identification of the BURP domain-containing genes in Phaseolus vulgaris

Plant-specific BURP domain-containing proteins have an essential role in the plant's development and stress responses. Although BURP domain-containing proteins have been identified in several plant species, genome-wide analysis of the BURP gene family has not been investigated in the common bean. In the present study, we identified 11 BURP family members in the common bean (Phaseolus vulgaris) genome with a comprehensive in silico analysis. Pairwise alignment and phylogenetic analyses grouped PvBURP members into four subfamilies [RD-22 like (3), PG1β-like (4), BNM2-like (3), and USP-like (1)] according to their amino acid motifs, protein domains and intron-exon structure. The physical and biochemical characteristics of amino acids, motif and intron-exon structure, and cis-regulatory elements of BURPs members were determined. Promoter regions of BURP members included stress, light, and hormone response-related cis-elements. Therefore, expression profiles of PvBURP genes were identified with in silico tools and qRT-PCR analyses under stress (salt and drought) and hormone treatment (ABA, IAA) in the current study. While significant activity changes were not observed in BURP genes in RNA-seq data sets related to salt stress, it was determined that some BURP genes were expressed differently in those with drought stress. We identified 12 different miRNA, including miRNA395, miRNA156, miRNA169, miRNA171, miRNA319, and miRNA390, targeting the nine PvBURP genes using two different in silico tools based on perfect or near-perfect complementarity to their targets. Here we present the first study to identify and characterize the BURP genes in common bean using whole-genome analysis, and the findings may serve as a reference for future functional research in common bean.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01052-9.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

Knockdown of miR393 Promotes the Growth and Biomass Production in Poplar

Short tandem target mimic (STTM), which is composed of two short sequences mimicking small RNA target sites, separated by a linker of optimal size, can block the functions of all members in a miRNA family. microRNA393 (miR393), which is one of the conserved miRNA families in plants, can regulate plant root growth, leaf development, plant architecture, and stress resistance. In order to verify the role of miR393 in the secondary growth of trees, we created its STTM transgenic poplar lines (STTM393). The expression of miR393 in STTM393 lines was reduced by over 10 times compared with the control plants. STTM393 lines showed promoted growth with about 20% higher, 15% thicker, and 2-4 more internodes than the control plants after 3 months of growth. The cross-section of the stems showed that STTM393 lines had wider phloem, xylem, and more cambium cell layers than control plants, and the lignin content in STTM393 lines was also higher as revealed by staining and chemical determination. Based on the transcriptome analysis, the genes related to the auxin signaling pathway, cell cyclin, cell expansion, and lignin synthesis had higher expression in STTM393 lines than that in control plants. The higher expression levels of FBL family members suggested that the auxin signaling pathway was strengthened in STTM393 lines to promote plant growth. Therefore, the knockdown of miR393 using the STTM approach provides a way to improve poplar growth and biomass production.

The PIF1-miR408-PLANTACYANIN repression cascade regulates light-dependent seed germination

Light-dependent seed germination is a vital process for many seed plants. A decisive event in light-induced germination is degradation of the central repressor PHYTOCHROME INTERACTING FACTOR 1 (PIF1). The balance between gibberellic acid (GA) and abscisic acid (ABA) helps to control germination. However, the cellular mechanisms linking PIF1 turnover to hormonal balancing remain elusive. Here, employing far-red light-induced Arabidopsis thaliana seed germination as the experimental system, we identified PLANTACYANIN (PCY) as an inhibitor of germination. It is a blue copper protein associated with the vacuole that is both highly expressed in mature seeds and rapidly silenced during germination. Molecular analyses showed that PIF1 binds to the miR408 promoter and represses miR408 accumulation. This in turn posttranscriptionally modulates PCY abundance, forming the PIF1-miR408-PCY repression cascade for translating PIF1 turnover to PCY turnover during early germination. Genetic analysis, RNA-sequencing, and hormone quantification revealed that PCY is necessary and sufficient to maintain the PIF1-mediated seed transcriptome and the low-GA-high-ABA state. Furthermore, we found that PCY domain organization and regulation by miR408 are conserved features in seed plants. These results revealed a cellular mechanism whereby PIF1-relayed external light signals are converted through PCY turnover to internal hormonal profiles for controlling seed germination.

Identification and Expression Analysis of miR160 and Their Target Genes in Cucumber

miR160 plays a crucial role in various biological processes by regulating their target gene auxin response factor (ARF) in plants. However, little is known about miR160 and ARF in cucumber fruit expansion. Here, 4 Csa-MIR160 family members and 17 CsARFs were identified through a genome-wide search. Csa-miR160 showed a closer relationship with those in melon. Phylogenetic analysis revealed that CsARFs were divided into four classes and most of CsARFs presented a closer evolutionary relationship with those from tomato. Putative cis-elements analysis predicted that Csa-MIR160 and CsARFs were involved in light, phytohormone and stress response, which proved that they might take part in light, phytohormone and stress signaling pathway by the miR160-ARF module. In addition, CsARF5, CsARF11, CsARF13 and CsARF14 were predicted as the target genes of Csa-miR160. qRT-PCR revealed that Csa-miR160 and their target gene CsARFs were differentially expressed in differential cucumber tissues and developmental stages. Csa-miR160d was only expressed in the expanded cucumber fruit. CsARF5, CsARF11 and CsARF13 exhibited the lower expression in the expanded fruit than those in the ovary, while, CsARF14 showed the reverse trend. Our results suggested that Csa-miR160d might play a crucial role in cucumber fruit expansion by negatively targeting CsARF5, CsARF11 and CsARF13. This is the first genome-wide analysis of miR160 in cucumber. These findings provide useful information and resources for further studying the role of miR160 and ARF in cucumber fruit expansion.

A prescient evolutionary model for genesis, duplication and differentiation of MIR160 homologs in Brassicaceae

MicroRNA160 is a class of nitrogen-starvation responsive genes which governs establishment of root system architecture by down-regulating AUXIN RESPONSE FACTOR genes (ARF10, ARF16 and ARF17) in plants. The high copy number of MIR160 variants discovered by us from land plants, especially polyploid crop Brassicas, posed questions regarding genesis, duplication, evolution and function. Absence of studies on impact of whole genome and segmental duplication on retention and evolution of MIR160 homologs in descendent plant lineages prompted us to undertake the current study. Herein, we describe ancestry and fate of MIR160 homologs in Brassicaceae in context of polyploidy driven genome re-organization, copy number and differentiation. Paralogy amongst Brassicaceae MIR160a, MIR160b and MIR160c was inferred using phylogenetic analysis of 468 MIR160 homologs from land plants. The evolutionarily distinct MIR160a was found to represent ancestral form and progenitor of MIR160b and MIR160c. Chronology of evolutionary events resulting in origin and diversification of genomic loci containing MIR160 homologs was delineated using derivatives of comparative synteny. A prescient model for causality of segmental duplications in establishment of paralogy in Brassicaceae MIR160, with whole genome duplication accentuating the copy number increase, is being posited in which post-segmental duplication events viz. differential gene fractionation, gene duplications and inversions are shown to drive divergence of chromosome segments. While mutations caused the diversification of MIR160a, MIR160b and MIR160c, duplicated segments containing these diversified genes suffered gene rearrangements via gene loss, duplications and inversions. Yet the topology of phylogenetic and phenetic trees were found congruent suggesting similar evolutionary trajectory. Over 80% of Brassicaceae genomes and subgenomes showed a preferential retention of single copy each of MIR160a, MIR160b and MIR160c suggesting functional relevance. Thus, our study provides a blue-print for reconstructing ancestry and phylogeny of MIRNA gene families at genomics level and analyzing the impact of polyploidy on organismal complexity. Such studies are critical for understanding the molecular basis of agronomic traits and deploying appropriate candidates for crop improvement.

TarDB: an online database for plant miRNA targets and miRNA-triggered phased siRNAs

Background

In plants, microRNAs (miRNAs) are pivotal regulators of plant development and stress responses. Different computational tools and web servers have been developed for plant miRNA target prediction; however, in silico prediction normally contains false positive results. In addition, many plant miRNA target prediction servers lack information for miRNA-triggered phased small interfering RNAs (phasiRNAs). Creating a comprehensive and relatively high-confidence plant miRNA target database is much needed.

Results

Here, we report TarDB, an online database that collects three categories of relatively high-confidence plant miRNA targets: (i) cross-species conserved miRNA targets; (ii) degradome/PARE (Parallel Analysis of RNA Ends) sequencing supported miRNA targets; (iii) miRNA-triggered phasiRNA loci. TarDB provides a user-friendly interface that enables users to easily search, browse and retrieve miRNA targets and miRNA initiated phasiRNAs in a broad variety of plants. TarDB has a comprehensive collection of reliable plant miRNA targets containing previously unreported miRNA targets and miRNA-triggered phasiRNAs even in the well-studied model species. Most of these novel miRNA targets are relevant to lineage-specific or species-specific miRNAs. TarDB data is freely available at http://www.biosequencing.cn/TarDB.

Conclusions

In summary, TarDB serves as a useful web resource for exploring relatively high-confidence miRNA targets and miRNA-triggered phasiRNAs in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07680-5.

MicroRNA775 regulates intrinsic leaf size and reduces cell wall pectin levels by targeting a galactosyltransferase gene in Arabidopsis

Plants possess unique primary cell walls made of complex polysaccharides that play critical roles in determining intrinsic cell and organ size. How genes responsible for synthesizing and modifying the polysaccharides in the cell wall are regulated by microRNAs (miRNAs) to control plant size remains largely unexplored. Here we identified 23 putative cell wall-related miRNAs, termed as CW-miRNAs, in Arabidopsis thaliana and characterized miR775 as an example. We showed that miR775 post-transcriptionally silences GALT9, which encodes an endomembrane-located galactosyltransferase belonging to the glycosyltransferase 31 family. Over-expression of miR775 and deletion of GALT9 led to significantly enlarged leaf-related organs, primarily due to increased cell size. Monosaccharide quantification, confocal Raman imaging, and immunolabeling combined with atomic force microscopy revealed that the MIR775A-GALT9 circuit modulates pectin levels and the elastic modulus of the cell wall. We also showed that MIR775A is directly repressed by the transcription factor ELONGATED HYPOCOTYL5 (HY5). Genetic analysis confirmed that HY5 is a negative regulator of leaf size that acts through the HY5-MIR775A-GALT9 repression cascade to control pectin levels. These findings demonstrate that miR775-regulated cell wall remodeling is an integral determinant of intrinsic leaf size in A. thaliana. Studying other CW-miRNAs would provide more insights into cell wall biology.

Interplay between miRNAs and lncRNAs: Mode of action and biological roles in plant development and stress adaptation

Plants employ sophisticated mechanisms to control developmental processes and to cope with environmental changes at transcriptional and post-transcriptional levels. MicroRNAs (miRNAs) and long noncoding RNAs (lncRNAs), two classes of endogenous noncoding RNAs, are key regulators of gene expression in plants. Recent studies have identified the interplay between miRNAs and lncRNAs as a novel regulatory layer of gene expression in plants. On one hand, miRNAs target lncRNAs for the production of phased small interfering RNAs (phasiRNAs). On the other hand, lncRNAs serve as origin of miRNAs or regulate the accumulation or activity of miRNAs at transcription and post-transcriptional levels. Theses lncRNA-miRNA interplays are crucial for plant development, physiology and responses to biotic and abiotic stresses. In this review, we summarize recent advances in the biological roles, interaction mechanisms and computational predication methods of the interplay between miRNAs and lncRNAs in plants.

Transcriptomic Profiling of Fe-Responsive lncRNAs and Their Regulatory Mechanism in Rice

Iron (Fe) deficiency directly affects crop growth and development, ultimately resulting in reduced crop yield and quality. Recently, long non-coding RNAs (lncRNAs) have been demonstrated to play critical regulatory roles in a multitude of pathways across numerous species. However, systematic screening of lncRNAs responding to Fe deficiency and their regulatory mechanism in plants has not been reported. In this work, 171 differently expressed lncRNAs (DE-lncRNAs) were identified based on analysis of strand-specific RNA-seq data from rice shoots and roots under Fe-deficient conditions. We also found several lncRNAs, which could generate miRNAs or act as endogenous target mimics to regulate expression of Fe-related genes. Analysis of interaction networks and gene ontology enrichment revealed that a number of DE-lncRNAs were associated with iron transport and photosynthesis, indicating a possible role of lncRNAs in regulation of Fe homeostasis. Moreover, we identified 76 potential lncRNA targets of OsbHLH156, a key regulator for transcriptional response to Fe deficiency. This study provides insight into the potential functions and regulatory mechanism of Fe-responsive lncRNAs and would be an initial and reference for any further studies regarding lncRNAs involved in Fe deficiency in plants.

The role of miRNA in plant-virus interaction: a review

Plant viruses affect crop production both quantitatively and qualitatively. The viral genome consists of either DNA or RNA. However, most plant viruses are positive single-strand RNA viruses. MicroRNAs are involved in gene regulation and affect development as well as host-virus interaction. They are non-coding short with 20-24 nucleotides long capable of regulating gene expression. The miRNA gene is transcribed by RNA polymerase II to form pri-miRNA which will later cleaved by Dicer-like 1 to produce pre-miRNA with the help of HYPONASTIC LEAVES1 and SERRATE which finally methylated and exported via nucleopore with the help of HASTY. The outcome of plant virus interaction depends on the effectiveness of host defense and the ability of a virus counter-defense mechanism. In plants, miRNAs are involved in the repression of gene expression through transcript cleavage. On the other hand, viruses use viral suppressors of RNA silencing (VSRs) which affect RISC assembly and subsequent mRNA degradation. Passenger strands, miRNA*, have a significant biological function in plant defense response as well as plant development.

MicroRNA annotation in plants: current status and challenges

Last two decades, the studies on microRNAs (miRNAs) and the numbers of annotated miRNAs in plants and animals have surged. Herein, we reviewed the current progress and challenges of miRNA annotation in plants. Via the comparison of plant and animal miRNAs, we pinpointed out the difficulties on plant miRNA annotation and proposed potential solutions. In terms of recalling the history of methods and criteria in plant miRNA annotation, we detailed how the major progresses made and evolved. By collecting and categorizing bioinformatics tools for plant miRNA annotation, we surveyed their advantages and disadvantages, especially for ones with the principle of mimicking the miRNA biogenesis pathway by parsing deeply sequenced small RNA (sRNA) libraries. In addition, we summarized all available databases hosting plant miRNAs, and posted the potential optimization solutions such as how to increase the signal-to-noise ratio (SNR) in these databases. Finally, we discussed the challenges and perspectives of plant miRNA annotations, and indicated the possibilities offered by an all-in-one tool and platform according to the integration of artificial intelligence.

Small Non-Coding RNAs at the Crossroads of Regulatory Pathways Controlling Somatic Embryogenesis in Seed Plants

Small non-coding RNAs (sncRNAs) are molecules with important regulatory functions during development and environmental responses across all groups of terrestrial plants. In seed plants, the development of a mature embryo from the zygote follows a synchronized cell division sequence, and growth and differentiation events regulated by highly regulated gene expression. However, given the distinct features of the initial stages of embryogenesis in gymnosperms and angiosperms, it is relevant to investigate to what extent such differences emerge from differential regulation mediated by sncRNAs. Within these, the microRNAs (miRNAs) are the best characterized class, and while many miRNAs are conserved and significantly represented across angiosperms and other seed plants during embryogenesis, some miRNA families are specific to some plant lineages. Being a model to study zygotic embryogenesis and a relevant biotechnological tool, we systematized the current knowledge on the presence and characterization of miRNAs in somatic embryogenesis (SE) of seed plants, pinpointing the miRNAs that have been reported to be associated with SE in angiosperm and gymnosperm species. We start by conducting an overview of sncRNA expression profiles in the embryonic tissues of seed plants. We then highlight the miRNAs described as being involved in the different stages of the SE process, from its induction to the full maturation of the somatic embryos, adding references to zygotic embryogenesis when relevant, as a contribution towards a better understanding of miRNA-mediated regulation of SE.

Genome-wide profiling of miRNAs during seed development reveals their functional relevance in seed size/weight determination in chickpea

MicroRNAs (miRNAs) are non-coding small RNAs that regulate gene expression at transcriptional and post-transcriptional levels. The role of miRNAs in seed development and seed size/weight determination is poorly understood in legumes. In this study, we profiled miRNAs at seven successive stages of seed development in a small-seeded and a large-seeded chickpea cultivar via small RNA sequencing. In total, 113 known and 243 novel miRNAs were identified. Gene ontology analysis revealed the enrichment of seed/reproductive/post-embryonic development and signaling pathways processes among the miRNA target genes. A large fraction of the target genes exhibited antagonistic correlation with miRNA expression. The sets of co-expressed miRNAs showing differential expression between the two cultivars were recognized. Known transcription factor (TF) encoding genes involved in seed size/weight determination, including SPL, GRF, MYB, ARF, HAIKU1, SHB1, KLUH/CYP78A5, and E2Fb along with novel genes were found to be targeted by the predicted miRNAs. Differential expression analysis revealed higher transcript levels of members of SPL and REVOLUTA TF families and lower expression of their corresponding miRNAs in the large-seeded cultivar. At least 19 miRNAs known to be involved in seed development or differentially expressed between small-seeded and large-seeded cultivars at late-embryogenesis and/or mid-maturation stages were located within known quantitative trait loci (QTLs) associated with seed size/weight determination. Moreover, 41 target genes of these miRNAs were also located within these QTLs. Altogether, we revealed important roles of miRNAs in seed development and identified candidate miRNAs and their target genes that have functional relevance in determining seed size/weight in chickpea.

sRNAanno-a database repository of uniformly annotated small RNAs in plants

Small RNAs (sRNAs) are essential regulatory molecules, and there are three major sRNA classes in plants: microRNAs (miRNAs), phased small interfering RNAs (phased siRNAs or phasiRNAs), and heterochromatic siRNAs (hc-siRNAs). Excluding miRNAs, the other two classes are not well annotated or available in public databases for most sequenced plant genomes. We performed a comprehensive sRNA annotation of 143 plant species that have fully sequenced genomes and next-generation sequencing sRNA data publicly available. The results are available via an online repository called sRNAanno ( www.plantsRNAs.org ). Compared with other public plant sRNA databases, we obtained was much more miRNA annotations, which are more complete and reliable because of the consistent and highly stringent criteria used in our miRNA annotations. sRNAanno also provides free access to genomic information for >22,721 PHAS loci and >22 million hc-siRNA loci annotated from these 143 plant species. Both miRNA and PHAS loci can be easily browsed to view their main features, and a collection of archetypal trans-acting siRNA 3 (TAS3) genes were annotated separately for quick access. To facilitate the ease of sRNA annotation, sRNAanno provides free service for sRNA annotations to the community. In summary, the sRNAanno database is a great resource to facilitate genomic and genetic research on plant small RNAs.