The evolution of microRNAs in plants

PmiR-Select® - a computational approach to plant pre-miRNA identification in genomes

Precursors of microRNAs (pre-miRNAs) are less used in silico to mine miRNAs. This study developed PmiR-Select® based on covariance models (CMs) to identify new pre-miRNAs, detecting conserved secondary structural features across RNA sequences and eliminating the redundancy. The pipeline preceded PmiR-Select® filtered 20% plant pre-miRNAs (from 38589 to 8677) from miRBase. The second filter reduced pre-miRNAs by 7% (from 8677 to 8045) through length limit to pre-miRNAs (70-300 nt) and miRNAs (20-24 nt). The 80% redundancy threshold was statistically the best, eliminating 55% pre-miRNAs (from 8045 to 3608). Angiosperms retained the highest number of pre-miRNAs and their families (2981 and 2202), followed by gymnosperms (362 and 271), bryophytes (183 and 119), and algae (82 and 78). Thirty-seven conserved pre-miRNA families happened among plant land clades, but none with algae. The PmiR-Select® was applied to the rice genome, producing 8536 pre-miRNAs from 36 families. The 80% redundancy threshold retained 3% pre-miRNAs (n = 264) from 36 families, valuable experimental and computational research resources. 14% (n = 1216) of 8536 were new pre-miRNAs from 19 new families in rice. Only 16 new sequences from six families overlapped (39 to 54% identities) with rice pre-miRNAs and five species on miRBase. The validation against mature miRNAs identified 8086 pre-miRNAs from 13 families. Eleven ones have already been recorded, but two new and abundant pre-miRNAs [miR437 (n = 296) and miR1435 (n = 725)] scattered in all 12-rice chromosomes. PmiR-Select® identified pre-miRNAs, decreased the redundancy, and discovered new miRNAs. These findings pave the way to delineating benchtop and computational experiments.

Regulation of abscisic acid receptor mRNA stability: Involvement of microRNA5628 in PYL6 transcript decay

Phytohormone signaling is fine-tuned by regulatory feedback loops. The phytohormone abscisic acid (ABA) plays key roles in plant development and abiotic stress tolerance. PYRABACTIN RESISTENCE 1/PYR1-LIKE/REGULATORY COMPONENT OF ABA RECEPTOR (PYR/PYL/RCAR) receptors sense ABA, and in turn, ABA represses their expression. Conversely, ABA induces expression of Type 2C PROTEIN PHOSPHATASES (PP2C) genes, which negatively regulate the ABA signaling pathway. This regulatory feedback scheme is likely important for modulating ABA signaling. Here, we provide insight into the mechanisms underlying the ABA-induced repression of PYR/PYL/RCAR expression in Arabidopsis (Arabidopsis thaliana). ABA time course analyses revealed strong and sustained repression of PYR/PYL/RCARs, suggesting that receptor gene regulation is an important step in resetting the ABA signaling pathway. Cordycepin-induced transcription inhibition showed that PYL1/4/5/6 mRNA destabilization is involved in the ABA-induced repression of these genes. Furthermore, genetic evidence indicated that decapping may play a role in PYL4/5/6 mRNA decay. We also provide evidence that the Arabidopsis-specific microRNA5628 (miR5628), which is transiently induced by the ABA core signaling pathway, guides PYL6 transcript cleavage in response to ABA. After cleavage, the resulting 5'- and 3'-cleaved fragments of PYL6 mRNA may be degraded by the XRN4 exoribonuclease. miR5628 is an evolutionary novelty that may enhance PYL6 mRNA degradation, along with decapping and XRN4 activity. Thus, regulating the stability of PYR/PYL/RCAR transcripts maintains ABA signaling homeostasis.

Integrative analysis of miRNA profile and degradome reveals post-transcription regulation involved in fragrance formation of Rosa rugosa

Rosa rugosa is renowned for its fragrant essential oils (EOs) including the primary volatile compounds such as terpenes (geraniol and citronellol) and 2-phenylethanol. While the role of miRNAs in plant secondary metabolism has been explored, their involvement in EOs metabolism remains largely unknown. Sequencing of the petals of R. rugosa identified 383 conserved miRNAs and 625 novel miRNAs including 53 miRNAs differentially expressed in a strong fragrance variety R. rugosa 'White Purple Branch'. Degradome sequencing predicted 1969 targets enriched in GO terms involved in the negative regulation of macromolecule metabolic process. Furthermore, 122 targets of differentially expressed miRNAs were enriched in phenylalanine metabolism and other KEGG pathways. A post-transcriptional regulation network of 52 miRNAs and 70 miRNA-transcription factor modules target terpene and 2-phenylethanol biosynthesis pathways. Six interactions including miR535f-RrHMGR, NOV146-RrNUDX1, miR166l-RrHY5 and miR156c-RrSPL2 were validated using RNA ligase-mediated RACE. Sequence alignment revealed that the NOV146-RrNUDX1 was conserved in the Rosa genus. Moreover, weaker silencing of RrNUDX1 by NOV146 contributed to the stronger fragrance of R. rugosa. These findings offer a comprehensive understanding of the post-transcriptional regulation involved in essential oil biosynthesis and identify candidate miRNAs for further genetic improvement of EO yields in R. rugosa.

Genome-Wide Transcriptional Response of Avocado to Fusarium sp. Infection

The avocado crop is relevant for its economic importance and because of its unique evolutionary history. However, there is a lack of information regarding the molecular processes during the defense response against fungal pathogens. Therefore, using a genome-wide approach in this work, we investigated the transcriptional response of the Mexican horticultural race of avocado (Persea americana var. drymifolia), including miRNAs profile and their possible targets. For that, we established an avocado-Fusarium hydroponic pathosystem and studied the response for 21 days. To guarantee robustness in the analysis, first, we improved the avocado genome assembly available for this variety, resulting in 822.49 Mbp in length with 36,200 gene models. Then, using an RNA-seq approach, we identified 13,778 genes differentially expressed in response to the Fusarium infection. According to their expression profile across time, these genes can be clustered into six groups, each associated with specific biological processes. Regarding non-coding RNAs, 8 of the 57 mature miRNAs identified in the avocado genome are responsive to infection caused by Fusarium, and the analysis revealed a total of 569 target genes whose transcript could be post-transcriptionally regulated. This study represents the first research in avocados to comprehensively explore the role of miRNAs in orchestrating defense responses against Fusarium spp. Also, this work provides valuable data about the genes involved in the intricate response of the avocado during fungal infection.

Emerging Orchestrator of Ecological Adaptation: m6A Regulation of Post-Transcriptional Mechanisms

Genetic mechanisms have been at the forefront of our exploration into the substrate of adaptive evolution and phenotypic diversification. However, genetic variation only accounts for a fraction of phenotypic variation. In the last decade, the significance of RNA modification mechanisms has become more apparent in the context of organismal adaptation to rapidly changing environments. RNA m6A methylation, the most abundant form of RNA modification, is emerging as a potentially significant player in various biological processes. Despite its fundamental function to regulate other major post-transcriptional mechanisms such as microRNA and alternative splicing, its role in ecology and evolution has been understudied. This review highlights the potential importance of m6A RNA methylation in ecological adaptation, emphasising the need for further research, especially in natural systems. We focus on how m6A not only affects mRNA fate but also influences miRNA-mediated gene regulation and alternative splicing, potentially contributing to organismal adaptation. The aim of this review is to synthesise key background information to enhance our understanding of m6A mechanisms driving species survival in dynamic environments and motivate future research into the dynamics of adaptive RNA methylation.

MicroRNA166: Old Players and New Insights into Crop Agronomic Traits Improvement

MicroRNA (miRNA), a type of non-coding RNA, is crucial for controlling gene expression. Among the various miRNA families, miR166 stands out as a highly conserved group found in both model and crop plants. It plays a key role in regulating a wide range of developmental and environmental responses. In this review, we explore the diverse sequences of MIR166s in major crops and discuss the important regulatory functions of miR166 in plant growth and stress responses. Additionally, we summarize how miR166 interacts with other miRNAs and highlight the potential for enhancing agronomic traits by manipulating the expression of miR166 and its targeted HD-ZIP III genes.

Genome-Wide Identification and Characterization of miRNAs and Natural Antisense Transcripts Show the Complexity of Gene Regulatory Networks for Secondary Metabolism in Aristolochia contorta

Aristolochia contorta Bunge is an academically and medicinally important plant species. It belongs to the magnoliids, with an uncertain phylogenetic position, and is one of the few plant species lacking a whole-genome duplication (WGD) event after the angiosperm-wide WGD. A. contorta has been an important traditional Chinese medicine material. Since it contains aristolochic acids (AAs), chemical compounds with nephrotoxity and carcinogenicity, the utilization of this plant has attracted widespread attention. Great efforts are being made to increase its bioactive compounds and reduce or completely remove toxic compounds. MicroRNAs (miRNAs) and natural antisense transcripts (NATs) are two classes of regulators potentially involved in metabolism regulation. Here, we report the identification and characterization of 223 miRNAs and 363 miRNA targets. The identified miRNAs include 51 known miRNAs belonging to 20 families and 172 novel miRNAs belonging to 107 families. A negative correlation between the expression of miRNAs and their targets was observed. In addition, we identified 441 A. contorta NATs and 560 NAT-sense transcript (ST) pairs, of which 12 NATs were targets of 13 miRNAs, forming 18 miRNA-NAT-ST modules. Various miRNAs and NATs potentially regulated secondary metabolism through the modes of miRNA-target gene-enzyme genes, NAT-STs, and NAT-miRNA-target gene-enzyme genes, suggesting the complexity of gene regulatory networks in A. contorta. The results lay a solid foundation for further manipulating the production of its bioactive and toxic compounds.

The roles of non-coding RNAs in male reproductive development and abiotic stress responses during this unique process in flowering plants

Successful male reproductive development is the guarantee for sexual reproduction of flowering plants. Male reproductive development is a complicated and multi-stage process that integrates physiological processes and adaptation and tolerance to a myriad of environmental stresses. This well-coordinated process is governed by genetic and epigenetic machineries. Non-coding RNAs (ncRNAs) play pleiotropic roles in the plant growth and development. The identification, characterization and functional analysis of ncRNAs and their target genes have opened a new avenue for comprehensively revealing the regulatory network of male reproductive development and its response to environmental stresses in plants. This review briefly addresses the types, origin, biogenesis and mechanisms of ncRNAs in plants, highlights important updates on the roles of ncRNAs in regulating male reproductive development and emphasizes the contribution of ncRNAs, especially miRNAs and lncRNAs, in responses to abiotic stresses during this unique process in flowering plants.

Blocking miR528 function promotes tillering and regrowth in switchgrass

MiRNAs have been reported to be the key regulators involving a wide range of biological processes in diverse plant species, but their functions in switchgrass, an important biofuel and forage crop, are largely unknown. Here, we reported the novel function of miR528, which has expanded to four copies in switchgrass, in controlling biomass trait of tillering number and regrowth rate after mowing. Blocking miR528 activity by expressing short tandem target mimic (STTM) increased tiller number and regrowth rate after mowing. The quadruple pvmir528 mutant lines derived from genome editing also showed such improved traits. Degradome and RNA-seq analysis, combined with in situ hybridization assay revealed that up-regulation of two miR528 targets coding for Cu/Zn-SOD enzymes, might be responsible for the improved traits of tillering and regrowth in pvmir528 mutant. Additionally, natural variations in the miR528-SOD interaction exist in C3 and C4 monocot species, implying the distinct regulatory strength of the miR528-SOD module during monocot evolution. Overall, our data illuminated a novel role of miR528 in controlling biomass traits and provided a new target for genetic manipulation-mediated crop improvement.

Identification of cell-type specific alternative transcripts in the multicellular alga Volvox carteri

BACKGROUND

Cell type specialization is a hallmark of complex multicellular organisms and is usually established through implementation of cell-type-specific gene expression programs. The multicellular green alga Volvox carteri has just two cell types, germ and soma, that have previously been shown to have very different transcriptome compositions which match their specialized roles. Here we interrogated another potential mechanism for differentiation in V. carteri, cell type specific alternative transcript isoforms (CTSAI).

METHODS

We used pre-existing predictions of alternative transcripts and de novo transcript assembly with HISAT2 and Ballgown software to compile a list of loci with two or more transcript isoforms, identified a small subset that were candidates for CTSAI, and manually curated this subset of genes to remove false positives. We experimentally verified three candidates using semi-quantitative RT-PCR to assess relative isoform abundance in each cell type.

RESULTS

Of the 1978 loci with two or more predicted transcript isoforms 67 of these also showed cell type isoform expression biases. After curation 15 strong candidates for CTSAI were identified, three of which were experimentally verified, and their predicted gene product functions were evaluated in light of potential cell type specific roles. A comparison of genes with predicted alternative splicing from Chlamydomonas reinhardtii, a unicellular relative of V. carteri, identified little overlap between ortholog pairs with alternative splicing in both species. Finally, we interrogated cell type expression patterns of 126 V. carteri predicted RNA binding protein (RBP) encoding genes and found 40 that showed either somatic or germ cell expression bias. These RBPs are potential mediators of CTSAI in V. carteri and suggest possible pre-adaptation for cell type specific RNA processing and a potential path for generating CTSAI in the early ancestors of metazoans and plants.

CONCLUSIONS

We predicted numerous instances of alternative transcript isoforms in Volvox, only a small subset of which showed cell type specific isoform expression bias. However, the validated examples of CTSAI supported existing hypotheses about cell type specialization in V. carteri, and also suggested new hypotheses about mechanisms of functional specialization for their gene products. Our data imply that CTSAI operates as a minor but important component of V. carteri cellular differentiation and could be used as a model for how alternative isoforms emerge and co-evolve with cell type specialization.

An Integrated Analysis of microRNAs and the Transcriptome Reveals the Molecular Mechanisms Underlying the Regulation of Leaf Development in Xinyang Maojian Green Tea (Camellia sinensis)

Xinyang Maojian (XYMJ) tea is one of the world's most popular green teas; the development of new sprouts directly affects the yield and quality of tea products, especially for XYMJ, which has hairy tips. Here, we used transcriptome and small RNA sequencing to identify mRNAs and miRNAs, respectively, involved in regulating leaf development in different plant tissues (bud, leaf, and stem). We identified a total of 381 conserved miRNAs. Given that no genomic data for XYMJ green tea are available, we compared the sequencing data for XYMJ green tea with genomic data from a closely related species (Tieguanyin) and the Camellia sinensis var. sinensis database; we identified a total of 506 and 485 novel miRNAs, respectively. We also identified 11 sequence-identical novel miRNAs in the tissues of XYMJ tea plants. Correlation analyses revealed 97 miRNA-mRNA pairs involved in leaf growth and development; the csn-miR319-2/csnTCP2 and miR159-csnMYB modules were found to be involved in leaf development in XYMJ green tea. Quantitative real-time PCR was used to validate the expression levels of the miRNAs and mRNAs. The miRNAs and target genes identified in this study might shed new light on the molecular mechanisms underlying the regulation of leaf development in tea plants.

Unveiling the biosynthesis, mechanisms, and impacts of miRNAs in drought stress resilience in plants

Drought stress is one of the most serious threats to sustainable agriculture and is predicted to be further intensified in the coming decades. Therefore, understanding the mechanism of drought stress tolerance and the development of drought-resilient crops are the major goals at present. In recent years, noncoding microRNAs (miRNAs) have emerged as key regulators of gene expressions under drought stress conditions and are turning out to be the potential candidates that can be targeted to develop drought-resilient crops in the future. miRNAs are known to target and decrease the expression of various genes to govern the drought stress response in plants. In addition, emerging evidence also suggests a regulatory role of long non-coding RNAs (lncRNAs) in the regulation of miRNAs and the expression of their target genes by a process referred as miRNA sponging. In this review, we present the regulatory roles of miRNAs in the modulation of drought-responsive genes along with discussing their biosynthesis and action mechanisms. Additionally, the interactive roles of miRNAs with phytohormone signaling components have also been highlighted to present the global view of miRNA functioning under drought-stress conditions.

Small RNAs: Efficient and miraculous effectors that play key roles in plant-microbe interactions

Plants' response to pathogens is highly complex and involves changes at different levels, such as activation or repression of a vast array of genes. Recently, many studies have demonstrated that many RNAs, especially small RNAs (sRNAs), are involved in genetic expression and reprogramming affecting plant-pathogen interactions. The sRNAs, including short interfering RNAs and microRNAs, are noncoding RNA with 18-30 nucleotides, and are recognized as key genetic and epigenetic regulators. In this review, we summarize the new findings about defence-related sRNAs in the response to pathogens and our current understanding of their effects on plant-pathogen interactions. The main content of this review article includes the roles of sRNAs in plant-pathogen interactions, cross-kingdom sRNA trafficking between host and pathogen, and the application of RNA-based fungicides for plant disease control.

miR-18a expression correlates with ATM and p53 levels and poor prognosis in lymphomas

microRNAs (miRNAs) are non-coding, endogenous, small-molecule RNAs. They are involved in cell proliferation, differentiation, apoptosis, and metabolism. Additionally, they play an essential role in the development and progression of various malignancies. Recent research has revealed that miR-18a plays an important role in cancer development. However, its role in lymphoma is not yet fully understood. In this study, we investigated the clinicopathological characteristics and potential functional roles of miR-18a in lymphomas. First, we predicted the potential downstream genes of miR-18a using miRTarBase software and subjected these downstream genes to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses to determine the potential mechanisms of action of these genes. We found that these target genes were closely related to cellular senescence, the p53 signaling pathway, and other signaling pathways. From the predicted downstream target genes, ATM and p53 were selected as the target genes; their deletion in patients with lymphoma was detected using the fluorescence in situ hybridization technique. The results showed that some patients with lymphoma have a deletion of the ATM and p53 genes. In addition, the deletion rates of ATM and p53 were positively correlated with the expression of miR-18a. Next, the expression levels of miR-18a and the deletion rates of ATM and p53 were used for correlation and prognostic analyses with patient clinical information. The findings revealed a significant difference in disease-free survival (DFS) between patients with lymphoma with ATM deletion and those with a normal ATM gene expression (p < 0.001). Moreover, a significant difference in overall survival (OS) and DFS between patients with p53 deletion and those with normal p53 expression was observed (p < 0.001). The results indicate that the deletion of ATM and p53 downstream of miR-18a is closely associated with the development of lymphoma. Thus, these biomarkers may serve as key prognostic biomarkers for lymphomas.

Miniature Inverted-Repeat Transposable Elements: Small DNA Transposons That Have Contributed to Plant MICRORNA Gene Evolution

Angiosperms form the largest phylum within the Plantae kingdom and show remarkable genetic variation due to the considerable difference in the nuclear genome size of each species. Transposable elements (TEs), mobile DNA sequences that can amplify and change their chromosome position, account for much of the difference in nuclear genome size between individual angiosperm species. Considering the dramatic consequences of TE movement, including the complete loss of gene function, it is unsurprising that the angiosperms have developed elegant molecular strategies to control TE amplification and movement. Specifically, the RNA-directed DNA methylation (RdDM) pathway, directed by the repeat-associated small-interfering RNA (rasiRNA) class of small regulatory RNA, forms the primary line of defense to control TE activity in the angiosperms. However, the miniature inverted-repeat transposable element (MITE) species of TE has at times avoided the repressive effects imposed by the rasiRNA-directed RdDM pathway. MITE proliferation in angiosperm nuclear genomes is due to their preference to transpose within gene-rich regions, a pattern of transposition that has enabled MITEs to gain further transcriptional activity. The sequence-based properties of a MITE results in the synthesis of a noncoding RNA (ncRNA), which, after transcription, folds to form a structure that closely resembles those of the precursor transcripts of the microRNA (miRNA) class of small regulatory RNA. This shared folding structure results in a MITE-derived miRNA being processed from the MITE-transcribed ncRNA, and post-maturation, the MITE-derived miRNA can be used by the core protein machinery of the miRNA pathway to regulate the expression of protein-coding genes that harbor homologous MITE insertions. Here, we outline the considerable contribution that the MITE species of TE have made to expanding the miRNA repertoire of the angiosperms.

Overview of Repressive miRNA Regulation by Short Tandem Target Mimic (STTM): Applications and Impact on Plant Biology

The application of miRNA mimic technology for silencing mature miRNA began in 2007. This technique originated from the discovery of the INDUCED BY PHOSPHATE STARVATION 1 (IPS1) gene, which was found to be a competitive mimic that prevents the cleavage of the targeted mRNA by miRNA inhibition at the post-transcriptional level. To date, various studies have been conducted to understand the molecular mimic mechanism and to improve the efficiency of this technology. As a result, several mimic tools have been developed: target mimicry (TM), short tandem target mimic (STTM), and molecular sponges (SPs). STTM is the most-developed tool due to its stability and effectiveness in decoying miRNA. This review discusses the application of STTM technology on the loss-of-function studies of miRNA and members from diverse plant species. A modified STTM approach for studying the function of miRNA with spatial-temporal expression under the control of specific promoters is further explored. STTM technology will enhance our understanding of the miRNA activity in plant-tissue-specific development and stress responses for applications in improving plant traits via miRNA regulation.

Comparative analysis of physiological variations and genetic architecture for cold stress response in soybean germplasm

Soybean (Glycine max L.) is susceptible to low temperatures. Increasing lines of evidence indicate that abiotic stress-responsive genes are involved in plant low-temperature stress response. However, the involvement of photosynthesis, antioxidants and metabolites genes in low temperature response is largely unexplored in Soybean. In the current study, a genetic panel of diverse soybean varieties was analyzed for photosynthesis, chlorophyll fluorescence and leaf injury parameters under cold stress and control conditions. This helps us to identify cold tolerant (V100) and cold sensitive (V45) varieties. The V100 variety outperformed for antioxidant enzymes activities and relative expression of photosynthesis (Glyma.08G204800.1, Glyma.12G232000.1), GmSOD (GmSOD01, GmSOD08), GmPOD (GmPOD29, GmPOD47), trehalose (GmTPS01, GmTPS13) and cold marker genes (DREB1E, DREB1D, SCOF1) than V45 under cold stress. Upon cold stress, the V100 variety showed reduced accumulation of H₂O₂ and MDA levels and subsequently showed lower leaf injury compared to V45. Together, our results uncovered new avenues for identifying cold tolerant soybean varieties from a large panel. Additionally, we identified the role of antioxidants, osmo-protectants and their posttranscriptional regulators miRNAs such as miR319, miR394, miR397, and miR398 in Soybean cold stress tolerance.

Plant responses to drought stress: microRNAs in action

Drought is common in most regions of the world, and it has a significant impact on plant growth and development. Plants, on the other hand, have evolved their own defense systems to deal with the extreme weather. The reprogramming of gene expression by microRNAs (miRNAs) is one of these defense mechanisms. miRNAs are short noncoding RNAs that have emerged as key post-transcriptional gene regulators in a variety of species. Drought stress modulates the expression of certain miRNAs that are functionally conserved across plant species. These characteristics imply that miRNA-based genetic changes might improve drought resistance in plants. This study highlights current knowledge of plant miRNA biogenesis, regulatory mechanisms and their role in drought stress responses. miRNAs functions and their adaptations by plants during drought stress has also been explained that can be exploited to promote drought-resistance among economically important crops.

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.

The era of cultivating smart rice with high light efficiency and heat tolerance has come of age

How to improve the yield of crops has always been the focus of breeding research. Due to the population growth and global climate change, the demand for food has increased sharply, which has brought great challenges to agricultural production. In order to make up for the limitation of global cultivated land area, it is necessary to further improve the output of crops. Photosynthesis is the main source of plant assimilate accumulation, which has a profound impact on the formation of its yield. This review focuses on the cultivation of high light efficiency plants, introduces the main technical means and research progress in improving the photosynthetic efficiency of plants, and discusses the main problems and difficulties faced by the cultivation of high light efficiency plants. At the same time, in view of the frequent occurrence of high-temperature disasters caused by global warming, which seriously threatened plant normal production, we reviewed the response mechanism of plants to heat stress, introduced the methods and strategies of how to cultivate heat tolerant crops, especially rice, and briefly reviewed the progress of heat tolerant research at present. Given big progress in these area, the era of cultivating smart rice with high light efficiency and heat tolerance has come of age.

Roles of RNA silencing in viral and non-viral plant immunity and in the crosstalk between disease resistance systems

RNA silencing is a well-established antiviral immunity system in plants, in which small RNAs guide Argonaute proteins to targets in viral RNA or DNA, resulting in virus repression. Virus-encoded suppressors of silencing counteract this defence system. In this Review, we discuss recent findings about antiviral RNA silencing, including the movement of RNA through plasmodesmata and the differentiation between plant self and viral RNAs. We also discuss the emerging role of RNA silencing in plant immunity against non-viral pathogens. This immunity is mediated by transkingdom movement of RNA into and out of the infected plant cells in vesicles or as extracellular nucleoproteins and, like antiviral immunity, is influenced by the silencing suppressors encoded in the pathogens' genomes. Another effect of RNA silencing on general immunity involves host-encoded small RNAs, including microRNAs, that regulate NOD-like receptors and defence signalling pathways in the innate immunity system of plants. These RNA silencing pathways form a network of processes with both positive and negative effects on the immune systems of plants.

Trichoderma atroviride hyphal regeneration and conidiation depend on cell-signaling processes regulated by a microRNA-like RNA

The ability to respond to injury is essential for the survival of an organism and involves analogous mechanisms in animals and plants. Such mechanisms integrate coordinated genetic and metabolic reprogramming events requiring regulation by small RNAs for adequate healing of the wounded area. We have previously reported that the response to injury of the filamentous fungus Trichoderma atroviride involves molecular mechanisms closely resembling those of plants and animals that lead to the formation of new hyphae (regeneration) and the development of asexual reproduction structures (conidiophores). However, the involvement of microRNAs in this process has not been investigated in fungi. In this work, we explore the participation of microRNA-like RNAs (milRNAs) molecules by sequencing messenger and small RNAs during the injury response of the WT strain and RNAi mutants. We found that Dcr2 appears to play an important role in hyphal regeneration and is required to produce the majority of sRNAs in T. atroviride. We also determined that the three main milRNAs produced via Dcr2 are induced during the damage-triggered developmental process. Importantly, elimination of a single milRNA phenocopied the main defects observed in the dcr2 mutant. Our results demonstrate the essential role of milRNAs in hyphal regeneration and asexual development by post-transcriptionally regulating cellular signalling processes involving phosphorylation events. These observations allow us to conclude that fungi, like plants and animals, in response to damage activate fine-tuning regulatory mechanisms.

Integration of mRNA and microRNA analysis reveals the molecular mechanisms underlying drought stress tolerance in maize (Zea mays L.)

Drought is among the most serious environmental issue globally, and seriously affects the development, growth, and yield of crops. Maize (Zea mays L.), an important crop and industrial raw material, is planted on a large scale worldwide and drought can lead to large-scale reductions in maize corn production; however, few studies have focused on the maize root system mechanisms underlying drought resistance. In this study, miRNA-mRNA analysis was performed to deeply analyze the molecular mechanisms involved in drought response in the maize root system under drought stress. Furthermore, preliminary investigation of the biological function of miR408a in the maize root system was also conducted. The morphological, physiological, and transcriptomic changes in the maize variety "M8186" at the seedling stage under 12% PEG 6000 drought treatment (0, 7, and 24 h) were analyzed. With prolonged drought stress, seedlings gradually withered, the root system grew significantly, and abscisic acid, brassinolide, lignin, glutathione, and trehalose content in the root system gradually increased. Furthermore, peroxidase activity increased, while gibberellic acid and jasmonic acid gradually decreased. Moreover, 32 differentially expressed miRNAs (DEMIRs), namely, 25 known miRNAs and 7 new miRNAs, and 3,765 differentially expressed mRNAs (DEMRs), were identified in maize root under drought stress by miRNA-seq and mRNA-seq analysis, respectively. Through combined miRNA-mRNA analysis, 16 miRNA-target gene pairs, comprising 9 DEMIRs and 15 DEMRs, were obtained. In addition, four metabolic pathways, namely, "plant hormone signal transduction", "phenylpropane biosynthesis", "glutathione metabolism", and "starch and sucrose metabolism", were predicted to have important roles in the response of the maize root system to drought. MiRNA and mRNA expression results were verified by real-time quantitative PCR. Finally, miR408a was selected for functional analysis and demonstrated to be a negative regulator of drought response, mainly through regulation of reactive oxygen species accumulation in the maize root system. This study helps to elaborate the regulatory response mechanisms of the maize root system under drought stress and predicts the biological functions of candidate miRNAs and mRNAs, providing strategies for subsequent mining for, and biological breeding to select for, drought-responsive genes in the maize root system.

SlMIR164A regulates fruit ripening and quality by controlling SlNAM2 and SlNAM3 in tomato

MiRNAs are important posttranscriptional regulators of plant development. Many miRNAs, such as the conserved miR164 species, are encoded by families of MIRNA genes, but the specific roles of individual MIRNA genes are largely undefined. Here, we characterize the functions and regulatory mechanisms of SlMIR164A, one of the primary genes of Sly-miR164, in tomato. We show that SlMIR164A is preferentially expressed at late stages of fruit development and plays a vital role in controlling fruit ripening and quality. Loss of function of SlMIR164A by CRISPR/Cas9-mediated mutagenesis results in accelerated fruit ripening and enhanced chloroplast development, which leads to altered sugar and organic acid contents and affects the nutritional quality of fruits. We also show that SlMIR164A modulates fruit ripening and quality through specific target genes, SlNAM2 and SlNAM3, which control key regulators of chloroplast function and fruit ripening processes. MIR164 genes have been shown to play conserved roles in regulating organ ageing, such as leaf senescence and fruit ripening, in a variety of plants, but whether and how their family members in tomato exert the same function remain to be elucidated. Our results reveal a previously undiscovered role of SlMIR164A in ripening control, which will further our understanding of the actions of MIR164 family, as well as the mechanisms of fruit ripening and quality control in tomato. Moreover, as loss of SlMIR164A exhibits minor impacts on organ morphology, our results can be leveraged in tomato breeding for specific manipulation of fruit ripening and quality to facilitate tomato improvement in agriculture.

The evolutionary history of small RNAs in Solanaceae

The Solanaceae or "nightshade" family is an economically important group with remarkable diversity. To gain a better understanding of how the unique biology of the Solanaceae relates to the family's small RNA (sRNA) genomic landscape, we downloaded over 255 publicly available sRNA data sets that comprise over 2.6 billion reads of sequence data. We applied a suite of computational tools to predict and annotate two major sRNA classes: (1) microRNAs (miRNAs), typically 20- to 22-nucleotide (nt) RNAs generated from a hairpin precursor and functioning in gene silencing and (2) short interfering RNAs (siRNAs), including 24-nt heterochromatic siRNAs typically functioning to repress repetitive regions of the genome via RNA-directed DNA methylation, as well as secondary phased siRNAs and trans-acting siRNAs generated via miRNA-directed cleavage of a polymerase II-derived RNA precursor. Our analyses described thousands of sRNA loci, including poorly understood clusters of 22-nt siRNAs that accumulate during viral infection. The birth, death, expansion, and contraction of these sRNA loci are dynamic evolutionary processes that characterize the Solanaceae family. These analyses indicate that individuals within the same genus share similar sRNA landscapes, whereas comparisons between distinct genera within the Solanaceae reveal relatively few commonalities.

TRUEE; a bioinformatic pipeline to define the functional microRNA targetome of Arabidopsis

Central to plant microRNA (miRNA) biology is the identification of functional miRNA-target interactions (MTIs). However, the complementarity basis of bioinformatic target prediction results in mostly false positives, and the degree of complementarity does not equate with regulation. Here, we develop a bioinformatic workflow named TRUEE (Targets Ranked Using Experimental Evidence) that ranks MTIs on the extent to which they are subjected to miRNA-mediated cleavage. It sorts predicted targets into high (HE) and low evidence (LE) groupings based on the frequency and strength of miRNA-guided cleavage degradome signals across multiple degradome experiments. From this, each target is assigned a numerical value, termed a Category Score, ranking the extent to which it is subjected to miRNA-mediated cleavage. As a proof-of-concept, the 428 Arabidopsis miRNAs annotated in miRBase were processed through the TRUEE pipeline to determine the miRNA 'targetome'. The majority of high-ranking Category Score targets corresponded to highly conserved MTIs, validating the workflow. Very few Arabidopsis-specific, Brassicaceae-specific, or Conserved-passenger miRNAs had HE targets with high Category Scores. In total, only several hundred MTIs were found to have Category Scores characteristic of currently known physiologically significance MTIs. Although non-exhaustive, clearly the number of functional MTIs is much narrower than many studies claim. Therefore, using TRUEE to numerically rank targets  directly on experimental evidence has given insights into the scope of the functional miRNA targetome of Arabidopsis.

Dosage-sensitive miRNAs trigger modulation of gene expression during genomic imbalance in maize

The genomic imbalance caused by varying the dosage of individual chromosomes or chromosomal segments (aneuploidy) has more detrimental effects than altering the dosage of complete chromosome sets (ploidy). Previous analysis of maize (Zea mays) aneuploids revealed global modulation of gene expression both on the varied chromosome (cis) and the remainder of the genome (trans). However, little is known regarding the role of microRNAs (miRNAs) under genomic imbalance. Here, we report the impact of aneuploidy and polyploidy on the expression of miRNAs. In general, cis miRNAs in aneuploids present a predominant gene-dosage effect, whereas trans miRNAs trend toward the inverse level, although other types of responses including dosage compensation, increased effect, and decreased effect also occur. By contrast, polyploids show less differential miRNA expression than aneuploids. Significant correlations between expression levels of miRNAs and their targets are identified in aneuploids, indicating the regulatory role of miRNAs on gene expression triggered by genomic imbalance.

Deciphering the role of miRNA in reprogramming plant responses to drought stress

Drought is the most prevalent environmental stress that affects plants' growth, development, and crop productivity. However, plants have evolved adaptive mechanisms to respond to the harmful effects of drought. They reprogram their: transcriptome, proteome, and metabolome that alter their cellular and physiological processes and establish cellular homeostasis. One of the crucial regulatory processes that govern this reprogramming is post-transcriptional regulation by microRNAs (miRNAs). miRNAs are small non-coding RNAs, involved in the downregulation of the target mRNA via translation inhibition/mRNA degradation/miRNA-mediated mRNA decay/ribosome drop off/DNA methylation. Many drought-inducible miRNAs have been identified and characterized in plants. Their main targets are regulatory genes that influence growth, development, osmotic stress tolerance, antioxidant defense, phytohormone-mediated signaling, and delayed senescence during drought stress. Overexpression of drought-responsive miRNAs (Osa-miR535, miR160, miR408, Osa-miR393, Osa-miR319, and Gma-miR394) in certain plants has led to tolerance against drought stress indicating their vital role in stress mitigation. Similarly, knock down (miR166/miR398c) or deletion (miR169 and miR827) of miRNAs has also resulted in tolerance to drought stress. Likewise, engineered Arabidopsis plants with miR165, miR166 using short tandem target mimic strategy, exhibited drought tolerance. Since miRNAs regulate the expression of an array of drought-responsive genes, they can act as prospective targets for genetic manipulations to enhance drought tolerance in crops and achieve sustainable agriculture. Further investigations toward functional characterization of diverse miRNAs, and understanding stress-responses regulated by these miRNAs and their utilization in biotechnological applications is highly recommended.

Effect of aneuploidy of a non-essential chromosome on gene expression in maize

The non-essential supernumerary maize (Zea mays) B chromosome (B) has recently been shown to contain active genes and to be capable of impacting gene expression of the A chromosomes. However, the effect of the B chromosome on gene expression is still unclear. In addition, it is unknown whether the accumulation of the B chromosome has a cumulative effect on gene expression. To examine these questions, the global expression of genes, microRNAs (miRNAs), and transposable elements (TEs) of leaf tissue of maize W22 plants with 0-7 copies of the B chromosome was studied. All experimental genotypes with B chromosomes displayed a trend of upregulated gene expression for a subset of A-located genes compared to the control. Over 3000 A-located genes are significantly differentially expressed in all experimental genotypes with the B chromosome relative to the control. Modulations of these genes are largely determined by the presence rather than the copy number of the B chromosome. By contrast, the expression of most B-located genes is positively correlated with B copy number, showing a proportional gene dosage effect. The B chromosome also causes increased expression of A-located miRNAs. Differentially expressed miRNAs potentially regulate their targets in a cascade of effects. Furthermore, the varied copy number of the B chromosome leads to the differential expression of A-located and B-located TEs. The findings provide novel insights into the function and properties of the B chromosome.

A fungal effector suppresses the nuclear export of AGO1-miRNA complex to promote infection in plants

SignificanceIncreasing evidence demonstrates that small RNAs can serve as trafficking effectors to mediate bidirectional transkingdom RNA interference (RNAi) in interacting organisms, including plant-pathogenic fungi systems. Previous findings demonstrated that plants can send microRNAs (miRNAs) to fungal pathogen Verticillium dahliae to trigger antifungal RNAi. Here we report that V. dahliae is able to secret an effector to the plant nucleus to interfere with the nuclear export of AGO1-miRNA complexes, leading to an inhibition in antifungal RNAi and increased virulence in plants. Thus, we reveal an antagonistic mechanism that can be exploited by fungal pathogens to counteract antifungal RNAi immunity via manipulation of plant small RNA function.

miR160 Interacts in vivo With Pinus pinaster AUXIN RESPONSE FACTOR 18 Target Site and Negatively Regulates Its Expression During Conifer Somatic Embryo Development

MicroRNAs (miRNAs) are key regulators of several plant developmental processes including embryogenesis. Most miRNA families are conserved across major groups of plant species, but their regulatory roles have been studied mainly in model species like Arabidopsis and other angiosperms. In gymnosperms, miRNA-dependent regulation has been less studied since functional approaches in these species are often difficult to establish. Given the fundamental roles of auxin signaling in somatic embryogenesis (SE) induction and embryo development, we investigated a previously predicted interaction between miR160 and a putative target encoding AUXIN RESPONSE FACTOR 18 in Pinus pinaster (PpARF18) embryonic tissues. Phylogenetic analysis of AUXIN RESPONSE FACTOR 18 (ARF18) from Pinus pinaster and Picea abies, used here as a model system of conifer embryogenesis, showed their close relatedness to AUXIN RESPONSE FACTOR (ARF) genes known to be targeted by miR160 in other species, including Arabidopsis ARF10 and ARF16. By using a luciferase (LUC) reporter system for miRNA activity in Arabidopsis protoplasts, we have confirmed that P. pinaster miR160 (ppi-miR160) interacts in vivo with PpARF18 target site. When the primary miR160 from P. pinaster was overexpressed in protoplasts under non-limiting levels of ARGONAUTE1, a significant increase of miR160 target cleavage activity was observed. In contrast, co-expression of the primary miRNA and the target mimic MIM160 led to a decrease of miR160 activity. Our results further support that this interaction is functional during consecutive stages of SE in the conifer model P. abies. Expression analyses conducted in five stages of development, from proembryogenic masses (PEMs) to the mature embryo, show that conifer ARF18 is negatively regulated by miR160 toward the fully developed mature embryo when miR160 reached its highest expression level. This study reports the first in vivo validation of a predicted target site of a conifer miRNA supporting the conservation of miR160 interaction with ARF targets in gymnosperms. The approach used here should be useful for future characterization of miRNA functions in conifer embryogenesis.

MicroRNA482/2118, a miRNA superfamily essential for both disease resistance and plant development

MicroRNAs (miRNAs) are a class of 21-24 nucleotides (nt) noncoding small RNAs ubiquitously distributed across the plant kingdom. miR482/2118, one of the conserved miRNA superfamilies originating from gymnosperms, has divergent main functions in core-angiosperms. It mainly regulates NUCLEOTIDE BINDING SITE-LEUCINE-RICH REPEAT (NBS-LRR) genes in eudicots, functioning as an essential component in plant disease resistance; in contrast, it predominantly targets numerous long noncoding RNAs (lncRNAs) in monocot grasses, which are vital for plant reproduction. Usually, miR482/2118 is 22-nt in length, which can trigger the production of phased small interfering RNAs (phasiRNAs) after directed cleavage. PhasiRNAs instigated from target genes of miR482/2118 enhance their roles in corresponding biological processes by cis-regulation on cognate genes and expands their function to other pathways via trans activity on different genes. This review summarizes the origin, biogenesis, conservation, and evolutionary characteristics of the miR482/2118 superfamily and delineates its diverse functions in disease resistance, plant development, stress responses, etc.

Characterization of maize miRNAs responsive to maize Iranian mosaic virus infection

MicroRNAs (miRNAs) play key regulatory roles in the plant's response to biotic and abiotic stresses and have fundamental functions in plant-virus interactions. The study of changes in miRNAs in response to virus infection can provide molecular details for a better understanding of virus-host interactions. Maize Iranian mosaic virus (MIMV) infects maize and certain other poaceous plants but miRNA changes in response to MIMV infection are unknown. In the present study, we compared the miRNA profiles of MIMV-infected and uninfected maize and characterized their predicted roles in response to the virus. Small RNA sequencing of maize identified 257 conserved miRNAs of 26 conserved families in uninfected and MIMV-infected maize libraries. Among them, miR395, miR166 and miR156 family members were highly represented. Small RNA data were confirmed using RT-qPCR. In addition, 33 potential novel miRNAs were predicted. The data show that 13 miRNAs were up-regulated and 113 were down-regulated in response to MIMV infection. Several of those miRNAs are known to be important in the response to plant pathogens. To determine the potential roles of individual miRNAs in response to MIMV, miRNA targets, predicted interactions with circular RNAs and comparative transcriptome data were analyzed. The expression profiles of different miRNAs in response to MIMV provide novel insights into the roles of miRNAs in the interaction between MIMV and maize plants.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-022-03134-1.

Integrated analysis of the transcriptome, sRNAome, and degradome reveals the network regulating fruit skin coloration in sponge gourd (Luffa cylindrica)

Sponge gourd fruit skin color is an important quality-related trait because it substantially influences consumer preferences. However, little is known about the miRNAs and genes regulating sponge gourd fruit skin coloration. This study involved an integrated analysis of the transcriptome, sRNAome, and degradome of sponge gourd fruit skins with green skin (GS) and white skin (WS). A total of 4,331 genes were differentially expressed between the GS and WS, with 2,442 down-regulated and 1,889 up-regulated genes in WS. The crucial genes involved in chlorophyll metabolism, chloroplast development, and chloroplast protection were identified (e.g., HEMA, CHLM, CRD1, POR, CAO, CLH, SGR, CAB, BEL1-like, KNAT, ARF, and peroxidase genes). Additionally, 167 differentially expressed miRNAs were identified, with 70 up-regulated and 97 down-regulated miRNAs in WS. Degradome sequencing identified 125 differentially expressed miRNAs and their 521 differentially expressed target genes. The miR156, miR159, miR166, miR167, miR172, and miR393 targeted the genes involved in chlorophyll metabolism, chloroplast development, and chloroplast protection. Moreover, a flavonoid biosynthesis regulatory network was established involving miR159, miR166, miR169, miR319, miR390, miR396, and their targets CHS, 4CL, bHLH, and MYB. The qRT-PCR data for the differentially expressed genes were generally consistent with the transcriptome results. Subcellular localization analysis of selected proteins revealed their locations in different cellular compartments, including nucleus, cytoplasm and endoplasmic reticulum. The study findings revealed the important miRNAs, their target genes, and the regulatory network controlling fruit skin coloration in sponge gourd.

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).

A unique miR775-GALT9 module regulates leaf senescence in Arabidopsis during post-submergence recovery by modulating ethylene and the abscisic acid pathway

The submergence-induced hypoxic condition negatively affects the plant growth and development, and causes early onset of senescence. Hypoxia alters the expression of a number of microRNAs (miRNAs). However, the molecular function of submergence stress-induced miRNAs in physiological or developmental changes and recovery remains poorly understood. Here, we show that miR775 is an Arabidopsis thaliana-specific young and unique miRNA that possibly evolved non-canonically. miR775 post-transcriptionally regulates GALACTOSYLTRANSFERASE 9 (GALT9) and their expression is inversely affected at 24 h of complete submergence stress. The overexpression of miR775 (miR775-Oe) confers enhanced recovery from submergence stress and reduced accumulation of RBOHD and ROS, in contrast to wild-type and MIM775 Arabidopsis shoot. A similar recovery phenotype in the galt9 mutant indicates the role of the miR775-GALT9 module in post-submergence recovery. We predicted that Golgi-localized GALT9 is potentially involved in protein glycosylation. The altered expression of senescence-associated genes (SAG12, SAG29 and ORE1), ethylene signalling (EIN2 and EIN3) and abscisic acid (ABA) biosynthesis (NCED3) pathway genes occurs in miR775-Oe, galt9 and MIM775 plants. Thus, our results indicate the role for the miR775-GALT9 module in post-submergence recovery through a crosstalk between the ethylene signalling and ABA biosynthesis pathways.

Degradome comparison between wild and cultivated rice identifies differential targeting by miRNAs

Background

Small non-coding (s)RNAs are involved in the negative regulation of gene expression, playing critical roles in genome integrity, development and metabolic pathways. Targeting of RNAs by ribonucleoprotein complexes of sRNAs bound to Argonaute (AGO) proteins results in cleaved RNAs having precise and predictable 5` ends. While tools to study sliced bits of RNAs to confirm the efficiency of sRNA-mediated regulation are available, they are sub-optimal. In this study, we provide an improvised version of a tool with better efficiency to accurately validate sRNA targets.

Results

Here, we improvised the CleaveLand tool to identify additional micro (mi)RNA targets that belong to the same family and also other targets within a specified free energy cut-off. These additional targets were otherwise excluded during the default run. We employed these tools to understand the sRNA targeting efficiency in wild and cultivated rice, sequenced degradome from two rice lines, O. nivara and O. sativa indica Pusa Basmati-1 and analyzed variations in sRNA targeting. Our results indicate the existence of multiple miRNA-mediated targeting differences between domesticated and wild species. For example, Os5NG4 was targeted only in wild rice that might be responsible for the poor secondary wall formation when compared to cultivated rice. We also identified differential mRNA targets of secondary sRNAs that were generated after miRNA-mediated cleavage of primary targets.

Conclusions

We identified many differentially targeted mRNAs between wild and domesticated rice lines. In addition to providing a step-wise guide to generate and analyze degradome datasets, we showed how domestication altered sRNA-mediated cascade silencing during the evolution of indica rice.

Supplementary Information

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

Integrated transcriptome and small RNA sequencing in revealing miRNA-mediated regulatory network of floral bud break in Prunus mume

MicroRNAs is one class of small non-coding RNAs that play important roles in plant growth and development. Though miRNAs and their target genes have been widely studied in many plant species, their functional roles in floral bud break and dormancy release in woody perennials is still unclear. In this study, we applied transcriptome and small RNA sequencing together to systematically explore the transcriptional and post-transcriptional regulation of floral bud break in P. mume. Through expression profiling, we identified a few candidate genes and miRNAs during different developmental stage transitions. In total, we characterized 1,553 DEGs associated with endodormancy release and 2,084 DEGs associated with bud flush. Additionally, we identified 48 known miRNAs and 53 novel miRNAs targeting genes enriched in biological processes such as floral organ morphogenesis and hormone signaling transudation. We further validated the regulatory relationship between differentially expressed miRNAs and their target genes combining computational prediction, degradome sequencing, and expression pattern analysis. Finally, we integrated weighted gene co-expression analysis and constructed miRNA-mRNA regulatory networks mediating floral bud flushing competency. In general, our study revealed the miRNA-mediated networks in modulating floral bud break in P. mume. The findings will contribute to the comprehensive understanding of miRNA-mediated regulatory mechanism governing floral bud break and dormancy cycling in wood perennials.

MicroRNA840 (MIR840) accelerates leaf senescence by targeting the overlapping 3'UTRs of PPR and WHIRLY3 in Arabidopsis thaliana

MicroRNAs negatively regulate gene expression by promoting target mRNA cleavage and/or impairing its translation, thereby playing a crucial role in plant development and environmental stress responses. In Arabidopsis, the MIR840 gene is located within the overlapping 3'UTR of the PPR and WHIRLY3 (WHY3) genes, both being predicted targets of miR840* and miR840, the short maturation products of MIR840. Gain- and loss-of-function of MIR840 in Arabidopsis resulted in opposite senescence phenotypes. The highest expression levels of the MIR840 precursor transcript pre-miR840 were observed at senescence initiation, and pre-miR840 expression is significantly correlated with a reduction in PPR, but not WHY3, transcript levels. Although a reduction of transcript level of PPR, but not WHY3 transcript levels were not significantly affected by MIR840 overexpression, its protein levels were strongly reduced. Mutating the cleavage sites or replacing the target sequences abolishes the miR840*/miR840-mediated degradation of PPR transcripts and accumulation of WHY3 protein. In support for this, concurrent knockdown of both PPR and WHY3 in wild-type plants resulted in a senescence phenotype resembling that of the MIR840-overexpressing plant. This indicates that both PRR and WHY3 are targets in the MIR840-mediated senescence pathway. Moreover, single knockout mutants of PPR and WHY3 show a convergent upregulated subset of senescence-associated genes, which are also found among those induced by MIR840 overexpression. Our data provide evidence for a regulatory role of MIR840 in plant senescence.

Recent origination of circular RNAs in plants

Circular RNA (circRNA) is a kind of new regulatory RNA with diverse biological functions. Numerous circRNAs have been identified in many plant species; however, evolution of plant circRNAs remains largely unknown. In this study, we assembled full-length sequences of 6519 rice (Oryza sativa) circRNAs and analyzed their conservation in another 46 plant species based on comparison of sequences and expression patterns. We found that, at the genomic level, 8.7% of the 6519 circRNAs were conserved in dicotyledonous plants and 49.1% in Oryza genus. Meanwhile, 57.8% of parental protein-coding genes of the rice circRNAs originated recently after divergence of monocotyledonous plants, implying recent origin of the majority of rice circRNAs, a conclusion further supported by the results based on analysis of 4663 full-length circRNAs in Arabidopsis thaliana. Accordingly, we proposed three models to address the origination of different types of circRNAs. Taken together, the results obtained in this study provide new insights for the evolutionary dynamics of plant circRNAs and candidate circRNAs for further functional exploration.

Biogenesis, Trafficking, and Function of Small RNAs in Plants

Small RNAs (sRNAs) encoded by plant genomes have received widespread attention because they can affect multiple biological processes. Different sRNAs that are synthesized in plant cells can move throughout the plants, transport to plant pathogens via extracellular vesicles (EVs), and transfer to mammals via food. Small RNAs function at the target sites through DNA methylation, RNA interference, and translational repression. In this article, we reviewed the systematic processes of sRNA biogenesis, trafficking, and the underlying mechanisms of its functions.

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.

Tobacco curly shoot virus Down-Regulated the Expression of nbe-miR167b-3p to Facilitate Its Infection in Nicotiana benthamiana

Tobacco curly shoot virus (TbCSV) belongs to the genus Begomovirus of the family Geminiviridae, and causes leaf curling and curly shoot symptoms in tobacco and tomato crops. MicroRNAs (miRNAs) are pivotal modulators of plant development and host-virus interactions. However, the relationship between TbCSV infection and miRNAs accumulation has not been well investigated. The present study was conducted to analyze different expressions of miRNAs in Nicotiana benthamiana in response to the infection of TbCSV via small RNAs sequencing. The results showed that 15 up-regulated miRNAs and 12 down-regulated miRNAs were differentially expressed in TbCSV infected N. benthamiana, and nbe-miR167b-3p was down-regulated. To decipher the relationship between nbe-miR167b-3p expression and the accumulations of TbCSV DNA, pCVA mediation of miRNA overexpression and PVX based short tandem target mimic (STTM) were used in this study. It was found that overexpression of nbe-miR167b-3p attenuated leaf curling symptom of TbCSV and decreased viral DNA accumulation, but suppression of nbe-miR167b-3p expression enhanced the symptoms and accumulation of TbCSV. PRCP, the target gene of nbe-miR167b-3p, was silenced in plants using VIGS and this weakened the viral symptoms and DNA accumulation of TbCSV in the plants. Overall, this study clarified the effect of nbe-miR167b-3p on plant defense during TbCSV infection, and provided a framework to reveal the molecular mechanisms of miRNAs between plants and viruses.

Serrate-Associated Protein 1, a splicing-related protein, promotes miRNA biogenesis in Arabidopsis

MicroRNAs (miRNAs) are essential regulators of gene expression in metazoans and plants. In plants, most miRNAs are generated from primary miRNA transcripts (pri-miRNAs), which are processed by the Dicer-like 1 (DCL1) complex along with accessory proteins. Serrate-Associated Protein 1 (SEAP1), a conserved splicing-related protein, has been studied in human and yeast. However, the functions of SEAP1 in plants remain elusive. Lack of SEAP1 results in embryo lethality and knockdown of SEAP1 by an artificial miRNA (amiR^SEAP1 ) causes pleiotropic developmental defects and reduction in miRNA accumulation. SEAP1 associates with the DCL1 complex, and may promote the interaction of the DCL1 complexes with pri-miRNAs. SEAP1 also enhances pri-miRNA accumulation, but does not affect pri-miRNA transcription, suggesting it may indirectly or directly stabilize pri-miRNAs. In addition, SEAP1 affects the splicing of some pri-miRNAs and intron retention of messenger RNAs at global levels. Our findings uncover both conserved and novel functions of SEAP1 in plants. Besides the role as a splicing factor, SEPA1 may promote miRNA biogenesis by positively modulating pri-miRNA splicing, processing and/or stability.

Bm-miR172c-5p Regulates Lignin Biosynthesis and Secondary Xylem Thickness by Altering the Ferulate 5 Hydroxylase Gene in Bacopa monnieri

MicroRNAs (miRNAs) are small non-coding, endogenous RNAs containing 20-24 nucleotides that regulate the expression of target genes involved in various plant processes. A total of 1,429 conserved miRNAs belonging to 95 conserved miRNA families and 12 novel miRNAs were identified from Bacopa monnieri using small RNA sequencing. The Bm-miRNA target transcripts related to the secondary metabolism were further selected for validation. The Bm-miRNA expression in shoot and root tissues was negatively correlated with their target transcripts. The Bm-miRNA cleavage sites were mapped within the coding or untranslated region as depicted by the modified RLM-RACE. In the present study, we validate three miRNA targets, including asparagine synthetase, cycloartenol synthase and ferulate 5 hydroxylase (F5H) and elucidate the regulatory role of Bm-miR172c-5p, which cleaves the F5H gene involved in the lignin biosynthesis. Overexpression (OE) of Bm-miR172c-5p precursor in B. monnieri suppresses F5H gene, leading to reduced lignification and secondary xylem thickness under control and drought stress. By contrast, OE of endogenous target mimics (eTMs) showed enhanced lignification and secondary xylem thickness leading to better physiological response under drought stress. Taken together, we suggest that Bm-miRNA172c-5p might be a key player in maintaining the native phenotype of B. monnieri under control and different environmental conditions.

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.

Evolution and Phylogeny of MicroRNAs - Protocols, Pitfalls, and Problems

MicroRNAs are important regulators in many eukaryotic lineages. Typical miRNAs have a length of about 22nt and are processed from precursors that form a characteristic hairpin structure. Once they appear in a genome, miRNAs are among the best-conserved elements in both animal and plant genomes. Functionally, they play an important role in particular in development. In contrast to protein-coding genes, miRNAs frequently emerge de novo. The genomes of animals and plants harbor hundreds of mutually unrelated families of homologous miRNAs that tend to be persistent throughout evolution. The evolution of their genomic miRNA complement closely correlates with important morphological innovation. In addition, miRNAs have been used as valuable characters in phylogenetic studies. An accurate and comprehensive annotation of miRNAs is required as a basis to understand their impact on phenotypic evolution. Since experimental data on miRNA expression are limited to relatively few species and are subject to unavoidable ascertainment biases, it is inevitable to complement miRNA sequencing by homology based annotation methods. This chapter reviews the state of the art workflows for homology based miRNA annotation, with an emphasis on their limitations and open problems.

Genome-Wide Identification of TCP Transcription Factors Family in Sweet Potato Reveals Significant Roles of miR319-Targeted TCPs in Leaf Anatomical Morphology

Plant-specific TCP transcription factors play vital roles in the controlling of growth, development, and the stress response processes. Extensive researches have been carried out in numerous species, however, there hasn't been any information available about TCP genes in sweet potato (Ipomoea batatas L.). In this study, a genome-wide analysis of TCP genes was carried out to explore the evolution and function in sweet potato. Altogether, 18 IbTCPs were identified and cloned. The expression profiles of the IbTCPs differed dramatically in different organs or different stages of leaf development. Furthermore, four CIN-clade IbTCP genes contained miR319-binding sites. Blocking IbmiR319 significantly increased the expression level of IbTCP11/17 and resulted in a decreased photosynthetic rate due to the change in leaf submicroscopic structure, indicating the significance of IbmiR319-targeted IbTCPs in leaf anatomical morphology. A systematic analyzation on the characterization of the IbTCPs together with the primary functions in leaf anatomical morphology were conducted to afford a basis for further study of the IbmiR319/IbTCP module in association with leaf anatomical morphology in sweet potato.